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Improved technique has enabled more precise measure- 
ments to be made of the effective cross section for scattering 
of hydrogen, water, and ammonia for low velocity protons 
and of hydrogen for low velocity H2* and H;*. From the 
measurements of proton scattering in hydrogen and 
ammonia and H;* scattering in hydrogen, and the use of 
certain assumptions, the average force laws can be evalu- 
ated and the energy of combination estimated in the first 
two cases. The proton affinities thus calculated for hydro- 
gen and ammonia agree with previous estimates. For H2* 
in H2 an evaluation of force law and an estimate of energy 


of combination cannot be made, but the large interaction 
is a phenomenon of some interest. Further increase in 
precision, separation of scattering from neutralization, and 
extension of measurements to lower velocities of ions is 
necessary also to evaluate the force law and proton 
affinity of water. As measurements of this kind can be 
made using many kinds of ions, each of which can be 
scattered in a great variety of gases, much information of 
value can be obtained such as force laws, proton affinity 
or absolute basicity of substances, energies of solvation of 
ions, etc. 





F scattering experiments with protons can be 

performed and interpreted in the manner 
suggested for apparent electron affinities by 
Simons and Seward,! an experimental method of 
general applicability will be available for the 
determination of proton affinities. These are of 
particular interest because they are, in effect, 
absolute basicities. In addition, a variety of other 
ions could most probably be used and scattered 
in a variety of gases. By the use of certain as- 
sumptions average force laws could be evaluated 
and energies of combination estimated. This 
would give a means of estimating the energies of 
solvation of ions. If precise enough measurements 
could be made, an estimate of the variation of 
force law with distance might be possible. 

In the electron scattering work mentioned 
above classical laws were used to estimate the 
force law. It is recognized that this is a doubtful 
procedure for electrons and that it gives only a 


. ‘Simons and Seward, J. Chem. Phys. 6, 790 (1938); 
Seward and Simons, ibid. 7, 2 (1939). 


first approximation. For a more rigorous treat- 
ment the change of scattering with angle should 
be determined and wave equations used in the 
correlation. However, for protons and heavier 
ions the method can be used with much less 
uncertainty. 

Previous experiments on the scattering of low 
velocity hydrogen ions by hydrogen as a function 
of voltage were made by Holzer.? Independently 
and at nearly the same time Ramsauer, Kollath, 
and Lilienthal* reported similar measurements. 
The results of these investigations agree only in 
order of magnitude. In the extreme low velocity 
end of the measurements the precision is not high. 


APPARATUS 


The apparatus used is shown in Fig. 1. It 
consisted of a source of low velocity hydrogen 
ions, elements for focusing the ions into a beam, a 

2 Holzer, Phys. Rev. 36, 1204 (1930). 


3 Ramsauer, Kollath and Lilienthal, Ann. d. Physik 8, 
709 (1931). 
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. Diagram of apparatus. 


magnet, defining holes for selecting a beam of 
the desired mass and velocity, and a scattering 
chamber and collecting Faraday cage for meas- 
uring the fraction of the ions deflected from the 
beam at a known pressure of gas in the scattering 
chamber. 

The source of hydrogen ions was an arc in 
hydrogen similar to the design of Lamar and 
Luhr.‘ This consisted of a filament, anode and 
water-cooled nickel cathode. Ions formed in the 
anode arc were accelerated toward the cathode 
and some escaped into the focusing region 
through a small hole (0.529-mm diameter) in the 
face of the cathode. This hole was sufficiently 
small so that a high vacuum could be maintained 
in the region next to the arc. The spiral filament 
was centered around the hole in the cathode and 
was oxide coated. It was mounted on glass- 
covered tungsten rods, which passed through 
holes in the face of the anode, and was sealed into 
the system through a ground glass joint. The 
female of this joint extended inside the cathode 
and insulated it from the anode. The anode was 
made to fit snugly inside its sleeve to lessen the 
loss by diffusion of activated particles from the 
arc. Electrical contact was made to the cathode 
by a wire soldered to the copper cooling jacket. 
The cathode was sealed to the ground glass joint 
and to the vacuum system by De Khotinsky 
cement. 

Just beyond the cathode face a series of coaxial 
nickel cylinders of alternate potential served to 
focus the ions escaping from the arc. From these a 
beam was selected by two holes located between 
the focusing units and the magnet chamber. This 


4 Lamar and Luhr, Phys. Rev. 46, 87 (1934). 


region was evacuated by a mercury pump and 
liquid-air trap with a speed for air of 20 liters per 
second. The construction of this focusing and 
beam defining unit is shown in the insert (a) of 
Fig. 1. The focusing units were made of {-inch 
nickel tubing fastened by nickel strips to two 
glass rods. These rods fit into nickel sleeves 
soldered into the can supporting the beam- 
defining holes. The unit was held in position, 
with respect to the hole in the cathode, by the 
glass rods fitting into sockets in the cathode face 
and by a De Khotinsky seal between the sup- 
porting can and surrounding glass vacuum sys- 
tem. The electrical leads to these elements were 
made through cement seals. 

The magnet chamber was made of sheet copper 
and was connected to the glass parts with cement 
seals. Care was taken that the ion beam was 
shielded by metal from the surface charge of the 
glass connecting sleeves. This chamber contained 
windows for observation of the alignment of 
source chamber, focusing, and scattering units. It 
was evacuated by a separate vacuum system 
similar to the one described. The U-shaped mag- 
net was made of 2000 turns of No. 18 copper wire 
wound on a soft iron core. The pole pieces were 
4 cm in diameter and 0.8 cm apart. The position 
of the pole pieces is indicated as “Magnet” in 
Fig. 1. 

The arrangement of the scattering unit is 
diagrammed in insert (b) of Fig. 1. The chambers 
were made from }-inch nickel tubing to which 
were soldered nickel end plates containing holes 
of the desired sizes. The regions between cham- 
bers were protected from stray potentials by 
nickel shields held at the potential of the 
scattering chambers. These chambers and shields 
were supported on glass-covered tungsten rods. 
The unit was held in position by resting in a 
socket in the surrounding glass envelope and by 
the rods being fastened to a ring which slipped 
around the magnet chamber. The glass envelope 
was sealed by cement to the magnet chamber. 
The scattering gas entered through a tube sealed 
into the glass envelope. A heavy iron magnetic 
shield surrounded the scattering unit. 

The three separate units of this scattering 
portion of the apparatus are labeled D, S and C 
for defining, scattering, and collecting chambers. 
D was 4.0 cm long and had holes 1.95 mm in 
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diameter. S was 3.63 cm long, the first hole was 
2.85 cm in diameter and the second 5.75 mm in 
diameter. C was 10 cm long and had an entrance 
hole 6 mm in diameter. The distance between 
adjacent units was 1 mm. 

The wiring diagram of electrical connections is 
shown in Fig. 2. Constancy of potential sources 
was obtained by the use of batteries, held under 
continuous charge, where the current drain was 
appreciable. Leads No. 1, 2, 3 and 4 indicate 
contacts to focusing units. The switch S; made 
possible the use of the galvanometer for current 
measurements for units preceding the magnet 
and also in conjunction with the amplifier for 
measurements to the scattering can S, and the 
collecting can C. The switch S,, allowed currents 
to S and C to be measured independently or 
simultaneously. 

The currents to S and C were measured by 
allowing them to flow over the high known 
resistance Rig and by balancing the potential 
drop thus established with the potentiometer. 
The position of balance was found by alternately 
applying the potential of the chamber being 
measured and the potential beyond the potenti- 
ometer to the grid of the amplifier tube. No 
deflection of the galvanometer indicated that the 
potentiometer measured the potential drop over 
R,s. As the circuit had a sensitivity of 1 mv per 
cm galvanometer deflection and the highest 
usable value of Rig was 100,000 megohms, a 
current sensitivity of 1X10-' amp. per cm was 
obtained. As several minutes were required for 
the system to reach equilibrium when this very 
high resistance was used, resistors of lower value 
were used on all measurements except those of 
the smallest currents. It was found that sufficient 
linearity existed between the potential drop over 
tke resistor and the galvanometer deflection so 
chat the latter could be used as a direct measure 
of current. 

Special shielded switches were constructed of 
pools of mercury in sealing wax. The contacts 
Were controlled with rods of sealing wax that 
extended beyond the shield. Sealing wax was 
folund to be a suitable insulator even on days of 
hijzh humidity. 

he hydrogen used was produced by electrolysis 
of a,KOH solution. It was passed over a glowing 
platitum spiral to remove oxygen and through 
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P.O; to remove water. The water was carefully 
degassed before use and the ammonia was from a 
commercial cylinder retained by ammonium 
nitrate. The pressures were measured by a 
McLeod gage and an uncertainty of perhaps 10 
percent exists in the pressure measurements of 
water and ammonia. This would affect the 
absolute value but not the form of the curves 
since the pressure was constant throughout the 
measurements. Pressures in the scattering cham- 


ber ranged from 1X10-‘ to 1X10-? mm Hg. 


OPERATION 


Used at hydrogen pressures of 0.2 to 0.5 mm, 
the arc was found to be stable with a filament 
current of 3 to 5 amperes, anode current of 0.9 to 
1.3 amperes and cathode current of 0.2 to 0.4 
ampere. The potential drop in the anode arc 
varied from 60 to 90 volts with change of 
filament current and hydrogen pressure. Neither 
the total beam current nor the ion composition in 
the beam was sensitive to cathode potential in 
the range from 40 to 200 volts. The beam current 
increased with anode current but was limited to 
1.3 amperes by overheating of the anode plate at 
higher currents. The maximum current obtained 
in the beam was approximately 5 microamperes. 


Fic. 2. Wiring diagram. 
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The beam was found to consist principally of 
H+, H.* and H;* ions. The relative proportions 
of the ions for a given set of conditions changed 
with seasoning of the arc, the proton percentage 
increasing. Conditions were not adjusted to give 
the high proton percentages obtainable with this 
arc. As used for the proton scattering experi- 
ments, the output consisted of approximately 30 
percent Ht, 40 percent H;+ and 30 percent Het. 

A pressure ratio of 500 to 1000 was measured 
between the arc and the pump evacuating the 


TABLE I. Scattering cross section of 162-volt H3* ions 
in hydrogen. 








Arc 


Pressure 0.3 mm; filament 4.0 amperes; anode 1.0 
ampere; 88-volt arc drop; cathode 0.30 ampere; 80- 
volt arc drop. Total beam current 5 microamperes. 


Focusing 


Potentials of focusing units: No. 1, 41 volts, Nos. 2 and 
4, —40 volts, Nos. 3 and 5, 94 volts. 

Potential of defining holes and magnet chamber: 0 volts. 

Current to magnet chamber: 0.5 microampere. 

Pressure in focusing chamber: 6X 10-* mm. 


Ion percentage 
Currents measured to collecting and scattering chambers 
at peaks in magnetic spectrum: 
Ion H+ H,* H;+ 
Current (to C+S) X10" amp. 13.7 13.2 16.5 
H;* velocity 
Pressure in scattering unit 2.7X10-* mm; 20 volts 
deceleration to defining can. Currents measured to 
collecting and scattering chambers as a function of 


retarding voltage applied between these cans and the 
defining can. 











RETARDING CURRENT TO RETARDING CURRENT TO 
POTENTIAL C+S POTENTIAL C+S 
Ovolts 16.0X10-" amp.| 142 volts 10.210- amp. 

120 14.5 141 13.5 

155 —1.2 140 14.0 

150 —1.1 135 15.0 

145 —1.0 130 15.3 

144 —1.0 0 17.0 

143 33 











Hence the H;?* ions are principally 162 volts. 
H;* beam deflection by hydrogen 
Current measured to collecting, collecting and scattering, 
and scattering chambers as a function of pressure in 
the scattering unit. Ratio R is C/(C+S). 











PRESSURE Cc C+8S Ss R 
2.7X10-4 mm | 23.5X10- amp. | 23.8X10-! amp. | 0.3X10- amp. | 0.990 
21.6 21.6 0.3 0.990 
7.7X103 18.2 e | 22.9 4.8 0.793 
7.0 17.2 22.3 4.6 
7.1 17.3 22.2 4.6 0.790 
10.7X10-3 14.7 21.6 6.8 0.685 
11.0 15.1 22.1 6.9 0.684 
11.1 15.1 22.1 7.0 0.684 
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region of the focusing cans. Undoubtedly in the 
region immediately in the path of the beam the 
pressure was considerably higher and caused 
scattering of the beam by collision and loss by 
neutralization. Hydrogen diffusing from the arc 
resulted in a pressure of 10~* mm in the scattering 
cans. This gave appreciable scattering at the 
lowest velocities of the ions. 

When the ions were passed into the magnet 
chamber at 150 volts, fifteen percent of the total 
beam current could be focused into the magnet 
chamber. At 50 volts, two percent of the beam 


TABLE II. Hydrogen cross section for scattering protons. 











PROTON PRESSURE 
VELOCITY MM Hg X 108 
VoLts at 0°C RATIO Loc 1/R a 
0.4 0.348 0.458 58.2 
6 3.1 .204 .600 
45 Ie i .142 
9 2.7 50 301 47.4 
40 .685 .164 
10 3.1 .390 409 62.3 
34 .710 .148 
11 5.3 .510 .509 51.7 
A5 78 .108 
14 23 52 .284 65.7 
23 825 .084 
14.5 6.1 .316 .508 49.2 
39 .760 .120 
iS 5.4 .346 461 46.3 
39 .840 .076 
a7 .615 .210 
21 $3 423 373 40.4 
39 .825 .086 
2.6 .663 178 
29 3. 492 .308 28.4 
= .906 .043 
40 14.0 365 439 19.9 
30 .950 .022 
7.1 .640 .194 
50 12.0 .500 301 16.2 
.28 .940 .027 
7.0 .620 .208 
50 11.8 483 316 16.8 
.28 .940 .027 
6.7 .695 .158 
66 11.4 575 .240 13.1 
.27 .955 .020 
6.7 .768 115 
80 11.4 .640 .194 11.3 
.27 .965 .016 
6.7 815 089 ¥ 
100 11.4 .705 151 8.78 | 
a .980 .009 ’ 
135 7.2 .885 .054 5.0 
ae .968 014 6.0%3 
6.7 .868 .061 
144 11.4 .780 .108 
ao .974 .012 
152 14.2 .753 .123 95 
Fe .980 .009 
200 19.0 .730 .136 SB .03 
.23 .978 .010 
237 19.0 755 .122 4.85 
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current could be made to enter the magnet 
chamber. Focusing potentials were adjusted by 
experiment to give maximum current into the 
final collecting chamber. Potential differences 
between adjacent focusing units had to be set 
below 200 volts to prevent arc formation. 

The ion velocity distribution from the arc was 
sufficient that the H+, H.*+ and H;* currents 
overlapped in the magnetic spectrum if the beam 
was passed between the magnet poles at less than 
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Fic. 3. Variation of log 1/R with pressure for protons in Hoe. 


50 volts. To obtain lower voltage particles in the 
scattering unit, further deceleration was applied 
between the magnet chamber and the defining 
chamber of the scattering unit. 

The velocity distribution of the ions into the 
scattering and collecting chambers was deter- 
mined by observing the current as a function of 
retarding potential to these chambers. With the 
chambers free+from gas, increased opposing 
potentials cause first a slight decrease in current 
and then a sharp decrease to a small negative 
current. The sharp decrease shows the true 
velocity range. The other effects probably result 
from secondary electrons and penetration of 
potentials to spread the beam inside the defining 
can of the scattering unit. Both of these effects 
become of negligible importance with the de- 
fining, scattering and collecting chambers at the 
same potential as they are during measurements 
of scattering cross sections. 

The energy distribution of the ion beam as set 
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Fic. 4. Hydrogen cross section for scattering protons. 
S. Holzer; VY Ramsauer, Rollath and Lilienthal; O 
RFS. 


by the magnet was in accord with the law 
AE/E=Ar/r=0.03, 


where AE is the distribution of energy, E is the 
energy of the particle, 7 is the radius of curvature 
of the beam passing through the magnet, and Ar 
is the change in radius permitted by the geometry 
of the slit system. 

A series of typical data is illustrated in Table I. 
From these data the scattering cross section of 
162-volt H; ions in hydrogen can be calculated. 


RESULTS 


In this apparatus scattering measurements 
were made for H+ in Hz, H2O and NHs;, and for 
H.*+ and H;* in He. The usual method was 
employed of determining the ratio of current from 
the collecting can to that from both collecting 
and scattering cans. From the variation of this 
ratio with pressure the effective cross-sectional 
area, a, wascalculated. When no gas was admitted 
to the scattering chamber, this ratio was about 98 
percent for voltages of about 100. It fell off 
slightly as the velocity decreased. 

In Table II is given this ratio of current for 
various hydrogen pressures used in the scattering 
of protons. Representative data from this table 
are plotted in Fig. 3 and show that values of a 
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Fic. 5. Hydrogen cross section for scattering Hs. A Holzer; 
O RFS. 


calculated from these data are reasonably inde- 
pendent of pressure in the range of pressures 
used. 

The results of the experiments are given in 
Table III, except H+ in He and are shown 
graphically in Figs. 4, 5 and 6. @ is given in 
cm?/cm? at 1 mm Hg and 0°C. Results obtained 
by other experimenters are also shown in the 
figures. In Fig. 4 the full line is drawn to the 
equation aV=850 which is an inverse second- 


TABLE III. Scattering gas. 








HYDROGEN WATER AMMONIA HYDROGEN 
Ht Ht H* H;3* 
VELOCITY VELOCITY VELOCITY VELOCITY 
VoLTs VoLTs a VoLtTs a VoLts 
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270 
235 
217 
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power potential law of attraction. If we accept 
the calculations of Hirschfelder® the equilibrium 
configuration of H;* is approximately that of an 
equilateral triangle, and the position of the H+ is 
1.5A from the center of the hydrogen molecule. 
Assuming that the above force law holds to the 
equilibrium position, a value of 3.5 electron volts 
is calculated for the proton affinity. Hirschfelder's 
calculations show this to be over 3.3 electron 
volts. In making a calculation of this kind a 
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Fic. 6. Water and ammonia cross sections for scattering 
protons. Upper curve—water; lower curve—ammonia. 
O first run; 0 check run. 


method similar to that used by Simons and 
Seward is employed. After the scattering law was 
found and the potential law thus established, the 
average scattering angle for the scattering cham- 
ber had to be calculated. This was done with 
classical mechanics, using the potential law and 
the dimensions of the scattering chamber. In Fig. 
5 the full line is drawn to the equation a V*/7 = 180. 
This would correspond to scattering from a 
spherically symmetrical potential reacting with 
the ion according to an average 3.5 power po- 
tential law. In Fig. 6 the upper curve represents 
the scattering of H+ in H.O with the full line 
drawn through the experimental points, the 
lower curve is for proton scattering in ammonia 
with the full line drawn to the equation at 
a V!= 300. This corresponds to an inverse fourth- 


5 Hirschfelder, J. Chem. Phys. 6, 795 (1938). 
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power potential law of attraction from spherically 
symmetrical particles. Assuming the H+ distance 
from the nitrogen in NH,* to be the same as the 
N-H distance in ammonia, i.e., 1.0A, then the 
proton affinity or basicity of ammonia is calcu- 
lated to be 9.4 electron volts. The assumption is 
used that the law of force holds down to the 
equilibrium position. The Born-Haber cycle has 
been used to calculate this value.* Values ob- 
tained from ammonium iodide, bromide, chloride, 
and fluoride are 202.7, 208.6, 209.0 and 221.0 kg 
cal., respectively. This agreement is much better 
than is justified by the measurements; and 
Bichowsky and Rossini’ have calculated the heat 
of formation of NH,* in the gas by using the 
values of lattice energies for NH,F, NH,Cl, 
NH,Br and NH,I given by Sherman, and their 
own values for the other thermal quantities 
needed, and have obtained this quantity in kg cal. 


6 Sherman, Chem. Rev. IF, 150 (1932). 
7 Bichowsky and Rossini, Thermochemistry of Chemical 
Substances (Reinhold Publishing Corp., 1936), p. 269. 
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of 0.0, —17.3, —21.2 and —42.2 respectively. 
No curve is drawn for the scattering of H2* in 
Ho, as the experimental values were dependent 
upon the amount of deceleration of the ions, and 
opposing potentials showed negative ions or 
electrons present. This shows that H,* has a 
remarkable ability to ionize hydrogen. This is in 
agreement with the observations of other 
investigators. 

In all these experiments there is probably 
some neutralization of the incoming ion stream, 
and the measurements have minor errors due to 
ionization, secondary electrons from the walls, 
and apparatus defects. If higher precision can be 
obtained, neutralization and ionization separated 
from scattering, and the measurements extended 
to lower velocities of ions, the force laws could be 
established with more certainty, particularly for 
substances such as water and ammonia. It is also 
conceivable that higher precision will enable the 
variation of force law with distance to be 
estimated. 
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Hindered Rotation in Halogenated Ethanes 


GEORGE GLOCKLER AND CHARLES SAGE* 
Department of Chemistry and Chemical Engineering, State University of Iowa, Iowa City, Iowa 


(Received February 24, 1941) 


If it is assumed on the basis of available evidence that the “‘staggered’’ form of the molecule 
is more stable than the “‘eclipsed”’ one, then the eight haloethanes investigated by means of 
the Raman effect show the following number of forms: CCl;CF; and CCl;CHFCI one form 
(13 and 18 lines); CF:CICF.Cl, CF.BrCF2Br, CF:CICF2H and CF.CICFCl. two forms 
(23, 22, 21 and 24 lines) and CF2CICHFCI and CF,BrCFCIBr three forms (29 and 28 lines). 
These forms are considered to be distinct rotational isomers. 


HE problem of hindered rotation of the 
carbon-carbon single bond in ethane and 
substituted ethanes has been studied by various 
investigators by means of the Raman effect,! 


*This article is based upon a thesis presented to the 
faculty of the Graduate School of the University of 
Minnesota by Charles Sage in partial fulfilment of the 
requirements for the degree of Doctor of Philosophy. 

'Kohlrausch and Kahovec Ber. 73, 159 (1940) and 
others; Mizushima and co-workers, Sci. Pap. Inst. Phys. 
Chem. Research (Tokyo) 29, 111 (1936) and others; 
Langseth and Bernstein, J. Chem. Phys. 8, 410 (1940). 


electron diffraction,? dipole moments’ and ther- 
modynamic considerations.‘ The present work 
consists of a study of the Raman spectra of some 
halogenated ethanes and gives further evidence 


2? Wierl, Ann. d. Physik 13, 453 (1932); Beach and 
co-workers, J. Am. Chem. Soc. 61, 3127 (1939) and others. 

3Smythe and co-workers, J. Am. Chem. Soc. 57, 979 
(1935) and others; Mizushima and co-workers, Sci. Pap. 
Inst. Phys. Chem. Research (Tokyo) 25, 159 (1934) and 
others; Beach and Stevenson, J. Chem. Phys. 6, 635 (1938). 

‘For references see: Pitzer, Symposium on Thermo- 
dynamic Properties of Hydrocarbons, Cincinnati Meeting 
of American Chemical Society (1940). 
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for the existence of ‘‘rotational isomerism’’ due 
to hindered rotation. 


EXPERIMENTAL RESULTS 


The liquid state Raman spectra of the com- 
pounds shown in Table I have been obtained 
with our usual equipment.® All of the substances 
except CCIl;CHFCI were obtained from Dr. 
Benning of the Jackson Laboratory of the E. I. 
du Pont de Nemours Company. The sample of 
CCl;CHFCI was obtained from Professor Henne 
of Ohio State University. We wish to thank 
Dr. Benning and Professor Henne for the loan 
of these substances. 

The two low boiling liquids, CCIF,CCIF»2 ‘and 
CCIF,CHF:, were kept below the boiling point 
by the use of the low temperature apparatus 
previously described. The temperatures main- 
tained during the exposure were —18 and 
—25°C, respectively. The temperature of the 
sample of CCIF2,CHFCI, which boils at 28°, was 
maintained at about 30° by means of a glass 
jacket through which cold air was circulated. 
The other substances were run at the regular 
temperature of the light furnace, which is about 
40°C. 

The Raman shifts of the substances are given 
in Table II. 


DISCUSSION OF EXPERIMENTAL DATA 


The substituted ethanes have eight atoms and 
consequently will have eighteen fundamental 
vibrational frequencies. The number of active 


TABLE I. 








CoMPOUND FREON NUMBER BOILING Point°C Purity* 


CCIF,CCI.F F-113 47.6 
CCIF:CCIF F-114 3.6 
CCIF,CHF, F-124 —10 
CCIF,CHFCIf F-123 28 
CF.BrCF.Br 46-47 
CF.BrCFCIBr 93-94 
CF;CCls F-113a 45.8 
CCI;CHFC1 116 











* The purity was specified as follows by Dr. Benning: C is fractionated 
in the laboratory in an efficient column to give a boiling range of about 
0.2°C. D is laboratory fractionated to give a boiling range of 0.5°C. 
J is refrigerator grade plant product, probably 99 percent pure or 

tter. 

t Structure uncertain. 


5 Glockler and Davis, J. Chem. Phys. 2, 881 (1934); 
Glockler and Bachmann, Phys. Rev. 54, 970 (1938). 
6 Glockler and Renfrew, Rev. Sci. Inst. 9, 306 (1938). 
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vibrations appearing in the Raman effect will 
depend, in any specific case, on the selection 
rules as determined from the symmetry of the 
molecular configuration. The symmetry of the 
molecular configuration will depend on _ the 
rotational orientation of the two halves of the 
molecule. There are two types of rotational 
isomers which must be considered, the so-called 
“staggered” and ‘“eclipsed’’ forms. However, 
from the evidence of other investigators®* on 
the electron diffraction, dipole moments, and 


TABLE II. Raman shifts of haloethanes. 








Av (cm-!; estimated relative intensity; exciting lines: 

a=4358; b=4046; c=4077). 

CF3;CCl;: 181.1(8,ab) ; 264.1(9,abc) ; 329.3(00,a) ; 
366.7(7,ab); 428.0(10,abc); 561.4(5,a); 657.2(0,a); 
711.7(8,abc); 791.5(0,a); 849.6(5,ab); 906.9(1,q); 
1214.1(2,a); 1250.4(1,ad). 

CCI;CHFCI: 160.1(3,a); 185.7(5,a) ; 235.7(6,ab) ; 
265.1(5,ab) ; 290.2(5,ab) ; 353.0(8,abc) ; 385.8(6,abc) ; 
452.3(10,abc) ; 613.9(7,ab) ; 767.7(0,a) ; 803.3(0,ad) ; 
832.8(2,ab) ; 842.7(7,ab) ; 1030.9(2,ab) ; 1096.7 (2,ab) ; 
1261.5(2,ab) ; 1321.5(2,a); 2983.8(5,ad). 

CF.CICF2CI: 171.4(4,a) ; 207.4(0,a) ; 253.1(7,ab) ; 
309.5(6,ab); 324.7(5,ab); 360.0(2,a); 399.9(1,a); 
443.9(10 broad,abc) ; 497.7(4,ab) ; 503.4(1,a); 
542.8(1,ab); 556.5(00,ab) ; 677.0(8,abc) ; 699.7 (8 abc); 
723.4(2,ab) ; 734.1(2,ab); 914.7(00,a); 1043.5(1,a); 
1055.9(1,a) ; 1063.7 (4,ab) ; 1111.6(0,a) ; 1172.5(2,ad); 
1269.6(3,ad). 

CF.CICFCl.: 167.4(3,a); 180.1(0,a) ; 202.8(3,a) ; 
241.6(6,ab); 251.6(1,a); 288.3(6,abc); 307.0(2,a); 
316.0(3,ab); 351.3(3,ab); 371.8(1,a); 393.6(3,ad); 
436.2(1,a); 442.5(9,ab); 459.9(9,ab); 507.2(5,a0); 
531.8(5,ab) ; 632.2(4,ab) ; 653.7(10,abc) ; 903.1(2,ad) ; 
1045.7(5,ab); 1089.9(0,a); 1111.6(0,a); 1169.2(0,a); 
1210.5(3,ab). 

CF.CICF:H: 93.0(?)(2?,b); 195.0(4,a) ; 254.1(1,a); 
322.3(5,ab); 340.7(3,a); 362.1(3,a); 434.1(9,ad); 
459.8(6,ab); 516.5(2,a); 560.6(3,a); 639.6(7,abc) ; 
673.9(10,abc); 701.1(2,ab); 828.7(6,ac); 910.2(0,a); 
997.4(2,ab) ; 1078.3(2,a) ; 1124.9(4,ab) ; 1214.9(3?,d); 
1352.5(4,ab) ; 2999.8(7,ab). 

CF.CICHFCI: 175.0(4,a) ; 206.7(1,a) ; 247.2(1,a); 
266.3(6,ab); 311.5(6,ab); 338.4(4,ab); 393.9(5,ad); 
427.8(10,ab) ; 439.1(10,ab) ; 457.5(8,ab) ; 477.6(5,ad) ; 
495.5(2,a); 598.9(5,ab); 622.9(6,ab); 630.0(1,ad); 
658.0(10,ab); 785.7(1,ab); 795.5(7,ab) ; 805.4(7,ad) ; 
847.9(3,ab); 862.2(0,a); 993.3(7,ab); 1057.2(0,c); 
1090.1(0,a); 1186.3(1,@); 1225.6(4,6); 1277.0(4,ad); 
1346.5(4,a) ; 2994.2(8,ab). 

CF.BrCF;Br: 63.4(4,ab); 179.0(10,ab) ; 272.0(1,ad) ; 
282.9(4,ab) ; 309.5(1,ab) ; 320.2(10,abc) ; 334.2(2,a8) ; 
342.5(6,abc); 367.3(7,abc); 465.3(2,a); 525.0(2,c); 
648.3(5,abc); 679.1(6,abc); 687.8(4,ab); 733.2(0,a); 
871.0(2,a) ; 964.4(0,a) ; 1003.7(3,ab) ; 1021.2(10,abc) ; 
1098.5(2,ab) ; 1162.4(2,ab); 1249.4(2,ab). 

CF.BrCFCIBr : 68.0(6,a); 173.1(10,ab) ; 197.4(8,a) ; 
208.2(2,a); 235.1(9,abc); 259.4(7,ab); 286.8(2,ab); 
296.4(7,abc); 310.5(8,ab); 331.0(6,a); 344.1(8,ad); 
409.3(5,ab); 441.0(4,ab); 468.5(3,ab); 510.7(4,abc) ; 
609.9(2,ab); 630.5(3,a); 637.7(7,ab); 728.8(1,ad); 
845.0(0,a); 859.4(4,a); 872.6(4,a); 889.6(5,ad); 
1013.6(10,abc) ; 1077.1(1,a); 1104.1(1,a); 
1161.1(0,a); 1197.9(2,a). 
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low temperature Raman spectra of substituted 
ethanes, it is reasonable to consider only staggered 
forms as stable configurations. With this re- 
striction, it is possible to divide the substances 
investigated in this work into three groups: 
(1) those which may have only one molecular 
configuration: CF;CCl; and CCl;CHFCI; (2) 
those which may have two rotational iso- 
mers: CF2CICF.Cl, CF,CICFCl:,, CF:CICF:H, 
CF.BrCF.Br; (3) those which may have three 
rotational isomers: CF2BrCFCIBr and CF;Cl- 
CHFCl. The configurations of the various 
isomers are shown in Fig. 1 together with the 
symmetry classification. Using Kohlrausch’s’ 
tables of selection rules, the number of Raman 
active vibrations for each isomeric configuration 
is obtained and, since the substance exists as a 
mixture of rotational isomers, we may tabulate 
the maximum possible number of Raman active 


_' Kohlrausch, Der Smekal-Raman Effekt, Erginzungsband 
(Springer, Berlin, 1938). 
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fundamental vibrations as shown in Table III. 
Comparing the number of lines observed in the 
spectrum with the number of active fundamental 
vibrations, it is seen that the number of observed 
lines increases as the number of possible isomeric 
forms increases. That there is not complete 
correspondence between the theoretical number 
of active vibrations and the number of observed 
lines may be due to three factors: (1) some 
observed lines are overtones or combination 
frequencies; (2) the abundance of some of the 


TABLE III. Selection rules for halogenated ethanes. 








LINES 
OBSERVED 


ACTIVE 


COMPOUND ISOMERS VIBRATIONS 


CCIsCF s 
CCI;CHFCI 
CF.CICF.Cl 





one form 1 13 
one form 18 18 


23 


trans 


9) 
ig {27 





CF .BrCF Br trans Cor 9 | 27 22 
C. | 18f 


CF.CICF.H trans Cin 18 \ 36 21 


Ci | 18f 


CF CICFCl, Cu} 18 | 
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ml 1g; 36 


24 
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CF.CICHFCI 
CF .BrCFCIBr 


3X 18=54 
3X18 =54 


three forms all C; 


three forms all C; 














isomers may be so small that only the most 
intense lines are obtained; (3) some of the 
frequencies may be nearly the same for all the 
isomers, thereby giving just a broadening of the 
lines. The broad and diffuse nature of the high 
wave number valency shifts, as contrasted with 
the sharpness of the low wave number deforma- 
tion shifts, indicated that the valency vibrations 
are not appreciably affected by the rotational 
orientation of the two halves of the molecule. 
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Raman Spectra of Acetylenes 


IV. Carbon Isotope Effect in Acetylenes 


Forrest F. CLEVELAND, Department of Physics 


AND 


M. J. Murray, Department of Chemistry, Illinois Institute of Technology, Chicago, Illinois 
(Received February 14, 1941) 


The carbon isotope shift in the frequency near 2200 cm™ has been calculated for H—C* 
=C”—H, CH;—C®=C8—H, CH;—C8 =C"®—H and CH;—C#=C”— CH; from the valence 
force equations, assuming a four-mass, linear oscillator. The results were —27, —26, —50 and 
— 33.5 cm“, respectively. Weak lines which may correspond to isotope shifts were observed 
for the monosubstituted acetylenes, 1-heptyne (—21 and —53 cm™), 4-methoxy-1-butyne 
(—24 and —53) and phenylethyne (—24 and —54), and for the disubstituted acetylenes, 
1-phenyl-1-butyn-3-one (—37), 1-phenyl-1-butynol-3 (—33), 1-phenyl-1-butyne (—29), 3- 
hexyne (—26), 3-octyne (—36) and 5-decyne (—38 cm™). Shifts of —27 and —33 cm“ for 
acetylene and dimethylacetylene, respectively, have been reported by Glockler and Renfrew. 





INTRODUCTION 


HE weak line at 1934 cm appearing with 

the strong 1961 line in the Raman spectrum 

of liquid acetylene has been explained by 

Glockler and Renfrew! as the triple-bond fre- 

quency of the isotopic molecule, HC®=C¥"H. 

Using the secular equation of Glockler and 

Wall,? they calculated the shift in frequency to 

be —32 cm-, as compared with the experimental 
value of —27. 

A similar explanation was given by Glockler 
and Renfrew* to account for the weak lines 2201 
and 2280 which appeared with the strong lines 
2235 and 2313 in the Raman spectrum of liquid 
dimethylacetylene. The shift, taken as the 
difference of the means of the (supposedly) 
resonance doublets was —33 cm~!. The shift in 
frequency calculated upon the assumption that 
the molecule constituted a two-body oscillator 
was — 20. Crawford‘ has recalculated the isotopic 
shift from a normal coordinate treatment, as 
follows: The shift in the fundamental near 


* Presented at the Philadelphia menting of the American 


Physical Society, December 1940. The first three papers 
of the series may be found in J. Am. Chem. Soc. 60, 2664 
(1938) ; 61, 3546 (1939); 62, 3185 (1940). 

1G. Glockler and M. M. Renfrew, J. Chem. Phys. 6, 
340 (1938). 

937) Glockler and F. T. Wall, J. Chem. Phys. 5, 813 
(1937). 

3G. Glockler and M. M. Renfrew, J. Chem. Phys. 6, 
408 (1938). 

4B. L. Crawford, Jr., J. Chem. Phys. 7, 555 (1939). 


2270 cm~! was calculated to be 46.4. The shift 
in the 1126 frequency was calculated to be 12.8, 
resulting in a 25.6 shift for the overtone fre- 
quency. The mean of the two shifts, 36.0 cm~', 
was then the expected difference in the means 
of the two doublets. 

More recently, Herzberg® has shown that the 
weak 2062 line accompanying the strong line at 
2094 in the Raman spectrum of HCN can be 
quantitatively explained as due to the isotopic 
molecule HC®N. 

Since the triple-bond line in the Raman 
spectra of substituted acetylenes is very intense, 
it occurred to the authors that the line due to 
the isotopic molecule might be observable on 
spectrograms of substituted acetylenes other 
than dimethylacetylene. 


CALCULATED ISOTOPE SHIFTS 


In order to have comparable values of the 
expected shifts, it seemed desirable to calculate 
them by use of the same potential function. 
Accordingly, calculations were made from the 
valence force equations,® for the molecules 
H—C¥=C®—-H, CH;—C®=C*®—H and 
CH;—C"=C”—H, assuming a four-mass, linear 
oscillator. Observed frequencies were used to 

5 G. Herzberg, J. Chem. Phys. 8, 847 (1940). 


6 Cf. K. W. F. Kohlrausch, Der Smekal-Raman Effett, 
Ergdénzungsband (Julius Springer, Berlin, 1938), p. 67. 
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calculate the force constants, which in turn were 
used to calculate the frequencies of the isotopic 
molecule. The shifts so obtained were —27, —26 
and —50 cm, respectively. For CH;—C® 
=C"—CH+;, the shift in the 2270 fundamental 
was calculated to be —43, the shift in the 
overtone of the 1126 frequency was calculated 
to be —24, thus giving an expected difference 
between the means of the doublets of —33.5 
cm~!. Comparison between the values calculated 
by this method, previously calculated values and 
observed values is made in Table I. 

In many cases substitution of other groups for 
the methyl group in monosubstituted acetylenes 
causes little change in the triple-bond frequency. 
For example, the frequency in methylacetylene 
is 2125 and in amylacetylene it is 2119 cm™. 
It therefore seems probable that the effective 
mass of such groups is not greatly different from 
that of the methyl group. This is in agreement 
with the conclusion reached by Harkins and 
Bowers,’ from a study of a series of aliphatic 
bromides, that the effective mass of the portion 
of the hydrocarbon chain (of two or more 
carbon atoms) which vibrates with respect to 
the bromine atom is 17. Using this value for the 
mass of the end group in R—C”®=C”—H, the 
calculated frequency was found to be only 3 cm™! 
lower than for a mass of 15. Moreover, the 
frequency of the isotopic molecule would also 
be lowered (probably by about the same amount) 
by this increase in mass of the end group and 
it thus seems reasonable to expect that the 
magnitude of the isotope shifts in other mono- 
substituted acetylenes would be nearly the same 
as those calculated above for methylacetylene. 

Similar arguments apply to disubstituted 
acetylenes. Here, however, the situation is 
complicated by the resonance interaction. The 
shift would be expected to be —43 cm~ if there 
were no resonance interaction and somewhat 
less than this value for cases in which there was 
such interaction. 


EXPERIMENTAL 
Monosubstituted acetylenes 


For the monosubstituted acetylenes, 1-hep- 
tyne, 4-methoxy-1-butyne and phenylethyne 


7™W. D. Harkins and H. E. Bowers, Phys. Rev. 38, 
1845 (1931). 
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weak lines were observed on the Raman spectro- 
grams which would correspond to the isotope 
shifts listed in the upper part of Table II. The 
lines corresponding to the smaller shifts were 
close to the strong triple-bond line, thus making 
their measurement difficult and the results 
somewhat uncertain. However, the half-width 
of the heavy line was only 9 cm~ so that the two 
lines should be just resolved if their separation 
is that calculated above. 


TABLE I. Calculated and observed carbon isotope shifts for 
acetylene, methylacetylene and dimethylacetylene. 





OBSERVED 
(cm~!) 


CALCULATED 
(cm~'!) 





PREs- 
ENT 
RE- 

SULTS 


—27 | -27 


COMPOUND RAMAN 








Acetylene 


—26 | — 
-50 | — 


Methylacetylene 














Dimethylacetylene 


—20 | —36.0 





@ Glockler and Renfrew. 
>’ Crawford. 


The lines corresponding to the larger shifts 
fell very close to the position (—57 cm-') of the 
triple-bond line excited by Hg 4347A and are 
thus open to question. Rough intensity measure- 
ments indicated that the intensity of Hg 4347A 
is 1/100 to 1/200 that of the Hg 4358A line, 
which was used for excitation of the spectra.’ 
Thus the intensity of the triple-bond line excited 
by Hg 4347A would be of the same order of 
magnitude as that of the line due to the isotopic 
molecule. It seems probable that the observed 
lines represent unresolved doublets of the two 
lines in question. If the two lines have the same 
influence upon the setting of the cross hair, the 
corrected isotope shifts would be —49, —49 
and —51i cm~, respectively. An attempt was 
made to use Hg 4047A for excitation in order to 
clear up this question. This, however, gave so 
much continuous background that it was im- 
possible to observe any lines of small intensity.** 

** Note added in proof.—While the manuscript was in 
the hands of the Editor, Mr. Robert D. Sieg, working in 
this laboratory, obtained the Raman spectrum of 1-hexyne 
(prepared by Mr. R. E. Dineen) with excitation by both 


Hg 4358A and Hg 4047A. The continuum was weak and 
the isotope shifts could easily be observed by both exciting 





392 F. F. CLEVELAND 
In the infra-red absorption spectrum of gaseous 
methylacetylene, Crawford® observed a strong 
band at 2151 cm~ and a weak band at 2101 cm 
which he interpreted as the combination fre- 
quency 642+1444. It seems that a reasonable 
alternative explanation would be the assignment 
of this band to the isotopic molecule CH;—C® 
=C"—H. The corresponding isotope shift would 
be —50 cm~—!. This value is included in Table I. 


Disubstituted acetylenes 


Weak lines which may correspond to fre- 
quencies of isotopic molecules were observed for 
the disubstituted acetylenes listed in Table III, 
which gives the displacements and relative 
intensities of the lines observed in the 2200-cm-! 
region for these compounds. The weak lines in 
question are the ones indicated by the intensity 


TABLE II. Carbon isotope shifts corresponding to observed 

















weak lines. 
SHIFT 
CoMPOUND MONOSUBSTITUTED | DISUBSTITUTED 
1-Heptyne —21 —53 
4-Methoxy-1-butyne —24 —53 
Phenylethyne —24 —54 
1-Phenyl-1-butyn-3-one —37 
1-Phenyl-1-butynol-3 — 33 
1-Phenyl-i-butyne —29 
3-Hexyne —26 
3-Octyne —36 
5-Decyne —3% 

















lines. The shifts obtained were —21 and —52 with 4358A 
and —25 and —47 cm™ with 4047A. This greatly increases 
the evidence in favor of the assignment of the weak lines 
to frequencies of the isotopic molecules. It is interesting to 
note that if the —52 shift (by 4358A) is corrected as above, 
one obtains the same value (—47 cm~) as with excitation 
by 4047A. 
8 B. L. Crawford, Jr., J. Chem. Phys. 8, 526 (1940). 


AND M. J. 
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TABLE III. Frequencies observed in the 2200-cm™ region for 
the disubstituted acetylenes. 








COMPOUND OBSERVED FREQUENCIES 





1-Phenyl-1-butyn-3-one 
1-Phenyl-1-butynol-3 
1-Phenyl-1-butyne 


2099(w), 2136(m), 2167(w), 2205(s) _ 
2200(w), 2233(s) 
2209(w), 2238(s) 


3-Hexyne® 2203(w), 2229(s), 2248(s), 2299(m) 
3-Octyne 2199(w), 2235(s), 2296(m) 
5-Decyne 2195(w), 2233(s), 2245(m), 2296(fs) 











@ w=weak, m =medium, fs =fairly strong, s =strong. 
+’ This compound was synthesized by Mr. H. J. Taufen. 


symbol w. In the case of 1-phenyl-1-butyn-3-one,’ 
it seems probable that 2205, 2136 is a resonance 
doublet and that 2167, 2099 is the corresponding 
resonance doublet for the isotopic molecule. If 
this is true, the corresponding isotope shift 
would be —37 cm~. Semiquantitative measure- 
ment of the relative intensities was made by 
taking 100-, 2-, 1- and }-minute exposures on 
the same film. Visual comparison of the strong 
line in the short exposures with the weak line 
in the long exposure indicated that the intensity 
ratio was between 100 : 1 and 50 : 1. The ratio 
expected from the relative abundance of the 
isotopes was 50:1, since approximately two 
percent of the molecules have the C™ isotope in 
the triple-bond position. 

In the other five cases, a weak line with a 
smaller Raman displacement was observed near 
the strong line. Failure to observe a correspond- 
ing weak line near the medium intensity line may 
be due to the smaller intensity. The isotope shifts 
listed in Table II for these five compounds repre- 
sent the differences in frequency between the 
strong and weak lines in Table III and are thus 
somewhat arbitrary. 


9 This compound was synthesized by Mr. R. E. Dineen. 
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A Partition Function for Normal Liquids 


JoHN WALTER AND HENRY EyRING 
Frick Chemical Laboratory, Princeton University, Princeton, New Jersey 


(Received January 29, 1941) 


A theory of liquids designed to account for rate as well as thermodynamic properties is 
given. By considering the molecules of the liquid to be distributed at random among solid-like 
and gas-like equilibrium positions, where the introduction of a new equilibrium position 
requires a volume increase of about 4 the molecular volume, it is possible to derive an explicit 
partition function for liquids. The model embodies the features previously used in the theory 
of viscosity and diffusion. The thermodynamic properties of argon, nitrogen and benzene are 


calculated and compared with experiment. 





INTRODUCTION 


E are interested here in constructing par- 
tition functions for liquids which embody 
the known relationship of the liquid to the solid 
and vapor states and which will give an adequate 
picture of the rate processes as well as the 
thermodynamic properties of liquids. This prob- 
lem has been previously examined from the 
reaction rate point of view by Roseveare, Powell 
and Eyring.'! These authors have formulate the 
essential features that a model of the liquid state 
must possess if it is to reproduce the observed 
rate processes correctly, and we have been able 
to retain these features in a partition function 
which gives satisfactory results for the thermo- 
dynamic properties. 

For the detailed requirements imposed by the 
rate processes we refer the reader to their paper, 
but we shall here summarize their results. The 
idea that empty equilibrium positions or ‘‘holes’’ 
are necessary for viscous flow shows, when 
applied to liquids, that the formation of such 
holes plus any additional activation energy re- 
quires about one-third the heat of vaporization 
for ordinary liquids and about one-tenth the heat 
of vaporization for metallic liquids. Since the 
entire heat of vaporization is required to form a 
hole the size of a molecule, it is clear that a 
volume considerably smaller than that normally 
occupied by a molecule is sufficient to form a new 
equilibrium position. The actual volume required 
for the formation of a new equilibrium position 
can be calculated by applying the theory of 
absolute reaction rates to the effect of pressure on 


'W. E. Roseveare, R. E. Powell and H. Eyring, J. App. 
Phys. (in press). 


393 


viscosity, since the derivative of the free energy 
of activation with respect to pressure is the 
increase in volume required to provide a new 
equilibrium position. Frisch, Kincaid and Eyring,’ 
using Bridgman’s* high pressure data, found in 
this way that a new equilibrium position for 
normal nonmetallic liquids occupied about 1/7 
the volume of a molecule and for metallic liquids 
about 1/23 of the molecular volume. 

These facts, together with the fact that a 
typical nonmetal expands about 10 percent 
on melting and a metal about 3 percent led 
Roseveare, Powell and Eyring! to suggest that 
substances acquire upon melting about 0.7 mole 
of new equilibrium positions per mole of mole- 
cules. If the mole of molecules were randomly 
distributed among the 1.7 moles of equilibrium 
positions, the entropy of melting would be 
approximately the two entropy units observed 
for many substances. 

The observation of Batschinski‘ that compen- 
sating changes of temperature and pressure which 
leave the volume unchanged do not greatly affect 
the viscosity substantiates the view that the 
chief activation free energy for viscosity is used 
in providing the equilibrium positions. A conse- 
quence of the above views is that liquids should 
have about the same viscosity at their melting 
points.! This is a rather good rule, and the 
exceptions seem to be cases where some special 
structural factor such as hydrogen bonds would 
be expected to complicate the simple picture. 


2 D. Frisch, J. Kincaid and H. Eyring, J. App. Phys. 11, 
75 (1940). 

3P, M. Bridgman, The Physics of High Pressure (Mac- 
millan, 1931). 

4 A. Batschinski, Zeits. f. physik. Chemie 84, 643 (1913). 
































































394 J. WALTER 
With these facts in mind, we can now formulate 
our problem more specifically. 


THE PARTITION FUNCTION 


Following the above discussion, we assume 
that the liquid has N equilibrium positions of the 
type present in the solid, and , new equilibrium 
positions of volume V,/n, with m, given by 
n,=nN(V—V,)/V;, where V, and V are the 
volumes per molecule in the solid and liquid, 
respectively. For a molecule in one of the solid 
type equilibrium positions we use the standard 
harmonic oscillator partition function 


e-0/2T \3 
f. exp (E,/RT) = (— en, exp (E,/RT), 


where 0=2@pebye, and E, is the potential energy 
of the solid. For a molecule in one of the gas-like 
equilibrium positions we use the partition func- 
tion for a molecule in a free volume V,/n, which 
is 
(2rmkT)? =| 

h® n 


Xexp (—6E,/RT), 


f, exp (—8E,/RT) = | 


where 6 is a function of volume to be discussed 
later. 

If the N molecules are distributed at random 
among the N+, positions, there will be a 
combinatory term in the complete partition func- 


TABLE I. Argon. 








VOLUME (cc) VAPOR PRESSURE (ATMOS.) 





T°K CALCULATED OBSERVED! CALCULATED OBSERVED! 
83.85 28.27 28.03 0.660 0.674 
87.44 28.81 28.69 1.000 1.000 
89.95 29.17 29.07 1.310 1.323 
97.71 30.42 30.15 2.735 2.671 
111.87 32.9 32.63 8.02 7.40 
122.34 35.9 35.08 15.04 13.58 
137.59 43 41.02 33.14 28.29 
147.93 52 51.68 ~48 43.19 


Melting point constants Critical constants 


CALc. Oss. CALc. Oss.1° 
Tm  82.9° 83.85° 7, 142° 150.66° 
AVm 3.14cc 3.05cc Ve 78.7 cc 75.26 cc 
ASm 3.40 3.35 P, 594atmos. 48.00 atmos. 
n= 8.35 
é= 60.08 
E,= 2067.4 cal. 
V.s= 24.98 cc® 
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tion equal to (V+,)!/N'!n,!. This assumption of 
random distribution of the molecules over the 
equilibrium positions leads to the correct par- 
tition function for the two limiting cases of solid 
and gas. It is not possible to say how accurate 
this expression is in the liquid range, although the 
fact that the thermodynamic properties of the 
liquid can be reproduced with a reasonable 
choice of the function 6 leads one to believe that 
it cannot be seriously incorrect. 

We have assumed that the most stable con- 
figuration will be the one in which the molecules 
are distributed among the two kinds of equi- 
librium positions in a ratio equal to the volume 
fractions of the two kinds, i.e., a fraction V,/ V in 
the solid type and a fraction (V—V,)/V in the 
gaseous type. Using this assumption, the com- 
plete partition function becomes 


f={(f. exp (E./RT))""!" 
(N +m)! 


x ( be — 6F, RT (V—V3)/V (\N____ 
fa exp ( /RT)) }! Wal 


From the activation energy for viscous flow, 
we know that the formation of a hole requires a 
greater fraction of the heat of vaporization than 
would be the case if the energy required were 
directly proportional to the volume of the hole. 
This indicates that 6 must be a positive quantity 
in the liquid range. If the partition function is to 
extrapolate smoothly to that of a perfect gas, 6 
must vanish for large volumes. The exact form of 
5 must be determined from experimental data. 

If we introduce the parameter y, defined by 
n,=~7N, we can write the Helmholtz free energy 
A in the simple form 


A= -Rr{ (in j+—) 


¥ 
n+~¥ 





|. 








+—(n h——) +i (147)+7 In 


-¥ 


For y=0, we have the equation for the solid ; for 
y—«, we have the equation for a perfect gas, 
provided that 6-0 as yo. The various 
thermodynamic functions are given explicitly in 
terms of A when a suitable form has been chosen 
for 6 and n has been given a numerical value. 








to | 
an 
foll 
clos 
for 
rare 


obs 
mel 
app 
duc 
the 
moc 
has 
nun 
X-ra 
Hor 
sum 
equi 
mol 
pict 
the 
U 
pres 
assu 
der 
xX 1¢ 
for 
with 
b= ( 
1/(n 
rath 
plete 
limit 
holes 
requ 
For 





en SP — e 





PARTITION FUNCTION 


APPLICATION TO LIQUID ARGON 


This type of partition function is first applied 
to liquid argon, as this is the simplest example of 
a normal liquid. Since the heavy rare gas liquids 
follow the law of corresponding states very 
closely ,5 we expect the partition function derived 
for argon to apply without change to the heavier 
rare gas liquids. 

The constant » was determined, from the 
observed entropy® and volume change® on 
melting, to be 8.35, indicating that on this picture 
approximately one mole of holes has been intro- 
duced at the melting point. Rice’ has discussed 
the entropy of melting of argon by means of a 
model in which 6 varies with volume and in which 
1 mole of larger size holes which later subdivide 
has been introduced, giving a’ coordination 
number 10 in the liquid in agreement with the 
x-ray diffraction results of Miller and Lark- 
Horowitz.* Lennard-Jones and Devonshire® as- 
sume that in the melting process a number of new 
equilibrium positions equal to the number of 
molecules has been introduced. On the above 
picture the number of holes is a linear function of 
the volume and may increase without limit. 

Using 6=60,° E, was calculated from the vapor 
pressure” of the solid at the melting point 
assuming the vapor to be represented by a van 
der Waals equation with constants a=1.346 
X10® atmos./cc, and 6=32.2 cc. The simplest 
form of 6 which would give fair agreement 
with the experimental data was found to be 
6=(V,/V)?(2y+1/(n—1))—. The small constant 
1/(n—1) in the denominator of this expression is 
rather arbitrary. If this value is used, the com- 
plete expression for the potential energy in the 
limit of small y’s is E= E,(1—vy), that is, the first 
holes to be introduced cost the amount of energy 
required to vaporize a molecule from the solid. 
For 

1 


2y>——-, 
n—1 


°K. Pitzer, J. Chem. Phys. 7, 583 (1939). 

°K. Clusius, Zeits. f. physik. Chemie B31, 459 (1936). 

70. K. Rice, J. Chem. Phys. 9, 121 (1941). 
(1940) P. Miller and K. Lark-Horowitz, Phys. Rev. 58, 207 
_*J. E. Lennard-Jones and A. F. Devonshire, Proc. Roy. 
Soc. A169, 317 (1939). 

Int. Crit. Tab. 
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Fic. 1. Calculated and observed specific heats of argon. 


which is nearly true in the normal liquid range, 
the potential energy is 


vy ts7° 
z=2,—-—(—) |. 
Vo 2nX\V 

The volume of the liquid under the observed 
vapor pressure was calculated at a given temper- 
ature by determining the volume for which 
Azr+pV_ was a minimum. From the value of Az, 
at this volume the vapor pressure was then 
calculated, using the relation (A+PV)riquia 
=(A+pV)cas. The vapor was assumed to be a 
van der Waals gas with the same values of the 
constants as those used in calculating E,. The 
melting point was calculated by determining the 
temperature at which the vapor pressures of the 
solid and liquid were equal, and the entropy of 
melting calculated from the calculated volume 
change on melting at this temperature. The 
results are given in Table I. 

The equation of state, given by the relation 
p=-—(0A/dV)z, is, for large liquid volumes 

RT en ie fi—In fs) 


p=—jn In 
V. y V2 


E.psV.\? 3 sV.\' 

gahG? “shed 
This is a van der Waals type equation of state, 
and the critical constants were calculated by 
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Fic. 2. The virial coefficient of argon. 


determining the temperature and volume at 
which the p— V isotherm has a point of inflection. 
The values so obtained are listed in Table I. 

The specific heat for this type of partition func- 
tion, to the approximation that the specific heat 
of the solid is classical, is C,=(3R/2)(1+V;/V). 
C, is given by the thermodynamic relation 


cemcver( 2) (22) mca 


The observed" and calculated values are plotted 
in Fig. 1. 

Experimental values were not available for the 
compressibility 8, but since the calculated coeff- 
cient of thermal expansion a is approximately 
correct, as may be seen from the calculated 
volumes, the percentage error in 6 is of the same 
order as that in C,—C,. 

_The above partition function extrapolates 
smoothly to that of a perfect gas. To see how well 
the equation of state represents the gaseous 
phase, the second virial coefficient, given by 


B(T)= 


y 2n—1 ——r E, 
4 ——+In f,—In f.-— 
2n : RT 


was plotted as a function of temperature in Fig. 


1! Tables Annuelles de Constantes 8, 101 (1931). 
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2, along with the observed” values and the 
second virial coefficient of the van der Waals 
equation. The agreement is comparable to that of 
the van der Waals equation of state. By plotting 
A as a function of volume, and determining the 
two points at which the curve had a common 
tangent, it would have been possible to determine 
the volumes of the liquid and gaseous phases at 
a given temperature, as well as the equilibrium 
pressure, directly from the partition function. 
The use of the van der Waals equation to 
represent the gas simplified the calculations 
considerably. 

In the above treatment the constant and the 
function 6 have been determined from the 
properties of the liquid. The value n=8.35 is 
approximately the value expected from previous 
studies of viscosity. 6 cannot be determined very 
uniquely; the form taken above gives a simple 
and reasonable expression for the potential 
energy of the liquid. The constants could be 
chosen in such a way that for y—>0 the correct 
value for the energy required to make a hole in 
the solid would be obtained. From the energy of 
activation for diffusion in the solid, this value is 
found to be of the order of one-half the energy of 
vaporization." 


Liguip NITROGEN 


Since nitrogen has a transition point in the 
solid and a normal entropy of melting, the 


TABLE II. Nitrogen. 








VAPOR PRESSURE 
VOLUME “- (ATMOS.) 
T°K CALC. Oss."° CALC. 


63.14 31.84 31.95 0.1229 
68.41 32.61 33.09 0.2918 
77.32 33.87 34.70 0.968 
90.63 36.05 37.70 3.434 
99.52 38.44 40.51 6.71 
111.79 42.1 46.08 14.63 


Melting point constants Critical constants 
CALC. Oss. CALC. Oss.” 
T, 62.86° 63.14°%5 7, 141.0° 126.0° 
AVm 249cc 2.64cc4# V, 90.3 cc 90.1 cc 
ASm 2.64 2.7315 P, 484atmos. 33.49 atmos. 
n= 8.35 
@=51' 
E,= 1660 cal. 
V,=29.31 cc'4 











R. H. Fowler and E. A. Guggenheim, Statistical 
Thermodynamics (Cambridge, 1939). 

13 A. E. Stearn and H. Eyring, J. Phys. Chem. 44, 955 
(1940). 
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PARTITION FUNCTION 


assumption was made that the molecules were 
rotating freely in both solid and liquid. The 
rotational degrees of freedom could thus be 
factored out of the partition function, as was the 
vibrational degree of freedom. Using these as- 
sumptions, the value of » calculated from the 


entropy of melting was 8.35; identical with that | 


found for argon. The same form of 6 was used as 
in the case of argon with constants chosen to give 
the melting temperature correctly ; the result was 
6=(V,/V)?(2.707+0.18)—. The liquid volumes, 
vapor pressures, and critical constants were 
calculated in the same manner as were those of 
argon; the results are given in Table II.'*1® 


BENZENE 


In treating benzene the rotational degrees of 
freedom must of course be included. Lord, 
Ahlberg and Andrews'* found that the specific 
heat of solid benzene could be accounted for by 
assuming that each benzene molecule had three 
torsional degrees of freedom with the same 
frequency as the lattice vibrations. The molecules 
were assumed to vibrate in phase, so that the 
specific heat could be represented by six Debye 
terms, with 6) =150. The specific heat of liquid 
benzene indicates that rotation is highly re- 
stricted in the liquid phase. We have, therefore, 
assumed that the molecules in the solid type 
equilibrium positions were undergoing torsional 
vibrations, while those in the gas type positions 
were rotating freely. The internal frequencies 
were assumed equal in solid, liquid and gas. For 
the solid, a value of @ was chosen so that the 
entropy from six Einstein terms at the melting 
point was equal to the entropy of six Debye 
terms with @)>=150, plus the entropy of lattice 
expansion.!® ” was taken to be 8.35, and a new 
value for 6 for the solid type positions in the 
liquid was chosen to give the entropy of melting. 
These values for the 6’s are 6=80.0, 6’ =55.7. E, 
was calculated from the vapor pressure of the 
solid at the melting point, and the constants in 6 
taken so as to give the melting temperature cor- 


*F. Simon, M. Ruhemann and W. Edwards, Zeits. f. 
physik. Chemie B6, 331 (1930). 

“ W. F. Giauque and J. O. Clayton, J. Am. Chem. Soc. 
55, 4882 (1933). 

*°R. C. Lord, J. R. Ahlberg and D. H. Andrews, J. Chem. 
Phys. 5, 649 (1937). 


FOR NORMAL LIQUIDS 


TABLE III. Benzene. 








T°K VOLUME (cc) VAPOR PRESSURE (MM) 
Lc. Os Oss.!° 


S.10 CALc. 


36.0 

95.2 

325.7 
753 
1319 
4400 
11500 





36.0 

94.4 

326.3 
760 
1335 
4335 
10640 


87.28 
89.36 
92.80 
95.96 
98.46 


87.24 
89.27 
93.32 
96.92 
100.3 


278.6 
298.1 
328.1 
353.3 
373.1 
423.1 110.0 106.8 
473.1 125.0 118.2 


T°K Cy(Lig.)-Cy(Gas)'? Cp-—Cy'8 

CALc. Oss. Cac. OBs. 
298.1 11.1 110 74 98 
328.1 109 10.8 7.9 9.7 


Melting point constants 
CALc. Oss.” 
Tm 278.6 278.6 
AVm 10.24cc 10.28cc V 
ASn 8.42 8.46 P 
n= 8.35 
6’=55.7 
E,= 11,334 cal. 
V.=77.00 cc!” 
Moments of inertia 
146,146,292 x 10-* g (cm)? 


B (Atmos, ~!)!8 
CALc. Oss 
155 10-6 97x 10-6 
211x10-§ 1221076 


Critical constants 
CALc. Opss.!0 
T 590 562 
289 cc 256 cc 
61.7 atmos. 47.9 atmos. 








rectly. For benzene, 6=(V,/V)?(0.87y+0.03)—. 
The partition function for the liquid does not 
join on smoothly to that of the solid in this case, 
but could be made to do so by having @ go 
continuously over into 6’. The calculated and 
observed properties are given in Table III.’ "8 

Preliminary calculations indicate that the 
same model may be successfully applied to the 
calculation of the thermodynamic properties of 
liquid metals. A value of ” considerably greater 
than the 8.35 used for the normal liquids is 
required, in agreement with the ideas previously 
developed from a study of the viscosity of liquid 
metals. 

From these results it appears that it is possible 
to base a thermodynamic treatment of liquids 
upon the same model which accounts for the 
observed rate processes. The volumes of the new. 
equilibrium positions introduced into the solid 
are of the same order as those previously de- 
termined from the study of viscosity and diffu- 
sion in liquids. 

In conclusion, the authors wish to express their 
appreciation to Mr. R. B. Simpson for assistance 
with certain of the numerical calculations. 


- 55 Kincaid and H. Eyring, J. Chem. Phys. 6, 620 
(1938). 

18 R, E. Gibson and J. Kincaid, J. Am. Chem. Soc. 60, 
511 (1938). 
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Second Virial Coefficients of Polar Gases* 
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An equation is given for the classical second virial coefficient of a polar gas in terms of the 
parameters appearing in an intermolecular potential energy which includes London and dipole 
attraction and inverse-power repulsion. The equation is successfully fitted to the data for HO 
and NHs with allowance for the small quantum correction, and the derived values of the 
London constant are in good accord with theoretical estimates. The defects of the model and 
its relation to the concept of the hydrogen bond are briefly discussed. 





INTRODUCTION 


HE characteristic physical properties of 
polar substances indicate that the electro- 
static interactions between the permanently 
asymmetric molecules comprise a very significant 
portion of the total van der Waals forces. 
Nevertheless, this indication has not yet been 
developed into a completely satisfactory quanti- 
tative correlation between theoretical and experi- 
mental values of the second virial coefficient. 
The attainment of such a correlation would not 
only be valuable in the theory of intermolecular 
forces, but might well be of assistance in the 
development of more accurate methods for calcu- 
lating the thermodynamic properties of gas 
mixtures containing polar components. The prac- 
tical importance of many high pressure equilibria 
involving polar gases makes the development of 
such methods highly desirable. 

The first attempts at a theory of polar gases! 
were distinctly unsatisfactory because at that 
time the quantum-mechanical theory of dis- 
persion forces? had not been developed. Thus 
electrostatic and induction forces alone were 
found insufficiently large to account for all of 
the van der Waals attraction unless unreasonable 
molecular parameters were assumed. In later 
treatments* ‘ dispersion forces were included and 
the electrostatic contribution was represented 
solely by the leading term (dipole-dipole inter- 


* Contribution No. 473 from the Research Labora- 
tory of Physical Chemistry, Massachusetts Institute of 
Technology. 

1W. H. Keesom, Physik. Zeits. 22, 129 (1921); H. 
Falkenhagen, Physik. Zeits. 23, 87 (1922). 

2 F. London, Zeits. f. physik. Chemie B11, 222 (1930). 

3H. Margenau, Phys. Rev. 36, 1782 (1930). 

4G. Briegleb, Zeits. f. physik. Chemie B23, 105 (1933). 


action), but a considerable disagreement with 
experiment was still apparent, especially in the 
case of steam. It was suggested* that the dis- 
crepancy was due to a rather large quadrupole 
moment in the water molecule, which was esti- 
mated to be about 5 X10-*6 e.s.u. 

In recent years it has become apparent that 
the properties of the condensed phases of many 
polar substances cannot be explained adequately 
if the molecular charge distribution is represented 
by a point dipole, and the concept of the hydro- 
gen bond has become very useful in correlating 
many such properties.’ It is therefore not un- 
reasonable to expect that dipole-dipole inter- 
action alone would also be insufficient to explain 
the properties of some gases. The association of 
hydrogen cyanide, for example, has recently been 
discussed from the hydrogen bond point of 
view.® The estimated quadrupole moment of the 
H.O molecule, 5X10-*° e.s.u., is nevertheless 
somewhat higher than values derived from simple 
reasonable models of the charge distribution. 
Furthermore, a fundamental error in the treat- 
ments of Margenau and Briegleb lies in their 
tacit assumption that electrostatic energy and 
dispersion energy contribute additively to the 
second virial coefficient, whereas in fact the inter- 
molecular potential energy appears in the phase 
integral as an exponential. A correct derivation 
from the simple dipole-dipole model leads to a 
much more satisfactory result, since the terms 
omitted through the assumption of additivity are 
of appreciable magnitude. 


5L. Pauling, The Nature of the Chemical Bond (Cornell 
University Press, 1940), second edition, Chapter IX. 

6 W. F. Giauque and R. A. Ruehrwein, J. Am, Chem. 
Soc. 61, 2626 (1939). 
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DEVELOPMENT OF EQUATIONS 


It is apparent that some economy must be 
exercised in the choice of a function to represent 
the intermolecular energy, especially for such 
anisotropic molecules as those to be considered 
here. All terms in the potential energy between 
two such molecules will depend on their orienta- 
tion, but some much more so than others. In 
order to avoid great complexity, it is convenient 
to treat all terms save the electrostatic energy 
between the permanent dipoles as independent 
of orientation. 

Accordingly, the energy of repulsion is taken 
to vary simply as an inverse power of the inter- 
molecular separation higher than six. The dis- 
persion energy is combined with the induction 
energy and written as a single spherically sym- 
metrical term varying as the inverse sixth power 
of the separation. Actually there are also higher 
terms involving eighth, tenth and higher inverse 
powers of the separation, but these are much 
smaller and it is possible to estimate their effect 
on the result. Although the dispersion and 
induction energies depend on orientation, they 
never change sign, so that it appears permissible 
in the dilute gas to use merely their values 
averaged over all orientations. Finally, the 
permanent dipole interaction, which depends 
most intimately on orientation, is given its usual 
explicit form as a function of the angular 
coordinates. Higher electrostatic terms are pur- 
posely neglected. 

The complete potential energy between a pair 
of similar polar molecules is thus written as 


U=)r-*—cr-*§— rg, (s>6) (1) 


where g=2 cos 6; cos 62:—sin 6; sin 62 cos ¢. 
In this expression » is the permanent dipole 
moment, 0; and 62 are the inclinations of the 
two dipole axes to the intermolecular axis, and ¢ 
is the azimuthal angle between them. 

The second virial coefficient in either of the 
equations of state 


pV/RT=1+B/V+::: (2) 
(3) 


may be found, for classical statistics, from the 


or pV=RT+Bp+::- 


COEFFICIENTS OF 


POLAR GASES 
well-known formula 


Beiass = 4N f f (1—e-U!*T)r°*drdQ, — (4) 
0 Q 


where 


fao=f f f sin 0; sin 02 d0,d0.dg=8r. 
2 0 0 0 


N is Avogadro’s number and & is Boltzmann’s 
constant. The dimensions of B thus calculated 
are those of volume per mole. 

The integration of Eq. (4) with the U defined 
by Eq. (1) is straightforward, and resembles the 
method employed by Lennard-Jones’ for non- 
polar gases. After part of the exponential is 
expanded and the definition of the gamma- 
function is utilized, the result of the integration 
becomes 


s—3 
; => (eNot/3)y4 r (=) 


o kSn/2 
-> Canp2*yins—6n—64) /2 , 
uk 
n=1 7 


0 


n G.3 (6n—6k-—3 
Cum (a) (= ‘ 
n! s 


5 


where 


The various new symbols appearing in Eq. (5) 
are defined in the following equations: 


a= /c, (6) 
y=c/kTo’, (7) 
p=n04/c, (8) 


1 1 m=k k as 
Gi=— | gao-—— F ()—— (0) 
Q 


Sr 2k+1"=° m+ 


and (sa) and (;.) are binomial coefficients. The 
series of Eq. (5) may be shown to converge for 
all finite values of y. When the dipole moment 
(and hence p) vanishes, Eq. (5) reduces to the 
result given by Lennard-Jones.’ If the exponent 
s is allowed to approach infinity while o is kept 
finite (which corresponds to representing the 
repulsion by impenetrable spheres of diameter ¢), 
Eq. (5) reduces to 


o kSn/2 
Beisss = (27 No*/3)[1—-Y YS Anp**y"], (10) 
n=1 k=0 


where 
G;. 1 


n! 2n—2k—1 
7 J. E. Lennard-Jones, Proc. Roy. Soc. A106, 463 (1924). 


Anz = (2x) 
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so that only integral powers of 1/T now appear. 
Equation (10) will be used in some of the com- 
parisons with experiment to follow. Further 
checks on the correctness of the result are 
afforded by letting either » or y become zero, 
in which cases Eq. (10) reduces to two expressions 
previously given by Keesom."® 

It is necessary to consider a correction to the 
classical result which is due to the quantization 
of rotational energy. This effect might be ex- 
pected to be fairly large for molecules with large 
permanent moments, but Margenau® and Kirk- 
wood’® have shown that it is rather small, though 
just large enough to receive consideration. Ex- 
plicit equations are given by Kirkwood for the 
simplest case of two identical diatomic polar 
molecules, and these may be used with sufficient 
confidence for the accuracy here required. The 
first several terms of the result given by use of 
Kirkwood’s equations (24a) and (31) were evalu- 
ated, and it was concluded that the quantum 
correction, AB,, to be added to the classical 
value of B given by Eqs. (5) or (10) could be 
written as 


AB, = (24No?/3)h2u'/36n°Jo%(RT)3+--+, (11) 


where J is taken as the geometric-mean moment 
of inertia about the two axes perpendicular to 
the dipole axis, and h is Planck’s constant. In no 
experimental case to be considered does AB, ever 
become as much as 1 percent of the total B. 
Since the uncertainty in the data is probably at 
least 1 percent, the quantum correction could 
really be omitted, although it will be included in 
the calculations to give better ideas of its 
magnitude. 

Before proceeding to fit Eq. (5) to experiment, 
TABLE I. Second virial coefficients of NH; in cc/mole, hard 


sphere model, Eq. (10). o=3.18A, c= 1.057 X 10-*8 
erg-cm®, p= 1.47 X 10-8 e.s.u. 








Obs.-CALc. 


—3.5 
+0.6 
+1.4 
+0.8 
—0.8 
—2.0 


T°(K) 


250 
285.7 
333.3 
400 
500 
600 


ABg Beale Bons 


—500.2 —503.7 
—316.1 —315.5 
—198.1 —196.7 
—119.9 —119.1 
— 65.9 — 66.7 
—38.6 —40.6 


Belass 


— 502.0 
—317.3 
— 198.9 
—120.4 
— 66.1 
— 38.7 











8 W. H. Keesom, Comm. Phys. Lab. Leiden, Suppl. 24B, 
32 (1912). 

®H. Margenau, Zeits. f. Physik 64, 584 (1930). 

10 J. G. Kirkwood, J. Chem. Phys. 1, 597 (1933). 
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it is necessary to consider further the exponent s 
appearing in the repulsive energy. Since the 
observations are limited to regions well below 
the Boyle temperature, the choice of s is very 
uncertain, but can have very little effect on the 
accuracy with which the data can be represented. 
Once a value of s has been chosen, however, the 
values of o and ¢ are determined from the data 
with considerable precision. It is, therefore, 
possible to exercise a little judgment in the final 
selection of s by deciding whether the corre- 
sponding values of ¢ and ¢ are in sufficient accord 
with other information about molecular size or 
attraction. Since one of the objects of the present 
work is to compare the experimentally derived 
values of c with those given by the London 
theory, the value of o alone may be used as a 
criterion of the adequacy of a chosen value of s. 

It should be mentioned that the application of 
the equations given above involves considerable 
tedium. For this reason, comparisons with experi- 
ment are restricted to the hard sphere model, 
Eq. (10), and to Eq. (5) with s=24. Values of s 
near 12 are customarily used with the Lennard- 
Jones equation for the light nonpolar gases, but 
Wohl" has shown that somewhat higher powers 
for substances with larger dispersion forces will 
bring the derived values of ¢ in better accord 
with London’s theory. In view of these considera- 
tions, the value s= 24 was adopted as reasonable 
for highly polar substances, a multiple of 6 being 
chosen for mathematical simplicity. For both 
Eqs. (5) and (10) it was found that powers of 1/T 
up to about ten were required in order to sum 
the series with sufficient accuracy at the lowest 
temperatures. Since Eq. (10) is somewhat easier 
to handle, it was investigated first. A systematic 
objective procedure of deriving best values of 
and c could not be developed, and it was neces- 
sary to re-evaluate the entire series for each set 
of trial values of o and c until satisfactory agree- 
ment with the observed values was attained. 
Consequently, the given values of « and c cannot 
be claimed to be the absolutely best values, 
although they give an excellent fit of the data. 
The experimental dipole moments, obtained from 
measurements of the dielectric constants of the 
gases, are employed throughout. 


1K, Wohl, Zeits. f. physik. Chemie B14, 36 (1931). 





VIRIAL 


TABLE II. Second virial coefficients of H2O, in cc/mole, hard 
sphere model, Eq. (10). c=3.16A, c=11.03X 10-*° 
erg-cm®, p= 1.83 X10 e.s.u. 








T°(K) 


400 
500 


Belass 


— 348.6 
— 168.0 
—99.8 


Bobs OBS.-CALC. 
— 346.9 
— 166.4 


—990 +0.1 
—-654 —2.1 





+0.5 





AMMONIA 


Keyes” has correlated the observed!* second 
virial coefficients of ammonia with the empirical 
equation 


Bnu;(cc/gram) = 2.4—(2316/T) 


exp (7.744X104/T2), (12) 


which covers, within the experimental error, the 
entire range of observation from 238° to 598°K. 
This expression has been used to find the given 
“observed”’ values of B for ammonia, the mo- 
lecular weight being taken as 17.031. In Table I 
is presented a comparison of several observed 
values with those calculated by Eq. (10) with 
o@=3.18A, c=1.057 X10-*8 erg-cm®, and 4» = 1.47 
X10-!8 e.s.u. (mean experimental value). The 
agreement is within the experimental error 
throughout, and the value of 3.18A for o corre- 
sponds exactly to setting the constant term in 
Eq. (12) equal to (27No*/3)/17.031. Since x-ray 
measurements!! reveal the N-—N distance in 
solid ammonia to be 3.38A, the hard sphere 
diameter here employed may be considered 
reasonable enough. Final calculations on am- 
monia for Eq. (5) with s = 24 have not been made, 
but several trials have indicated that for this 
model ¢ and ¢ remain practically unchanged. 


STEAM 


Measurements of the second virial coefficient 
of steam!® extend from 311° to 733°K, but com- 
parisons are here restricted to temperatures 
above 400°K, to avoid prohibitive labor in sum- 
ming the series of Eqs. (5) or (10). Keyes, Smith 


2 F. G. Keyes, J. Am. Chem. Soc. 60, 1761 (1938). 

8 C. H. Meyers and R. S. Jessup, Refrig. Eng. 11, 345 
(1925); J. A. Beattie and C. K. Lawrence, J. Am. Chem. 
Soc. 52, 6 (1930). 

4H. Markand E. Pohland, Zeits. f. Krist. 61, 532 (1925). 

FG. Keyes, L. B. Smith and H. T. Gerry, Proc. Am. 
Acad. Arts Sci. 70, 319 (1936); S. C. Collins and F. G. 
Keyes, ibid. 72, 283 (1938). 
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and Gerry give the correlative empirical equation 


Bu.0(cc/gram) = 1.89 — (2641.6/T) 


Xexp (1.858 10°/T?), (13) 


which was used, together with the molecular 
weight 18.016, to find the given ‘‘observed”’ 
values. The measurements by Collins and Keyes'® 
of low pressure Joule-Thomson coefficients from 
38° to 125°C led them to correct Eq. (13) 
slightly, but above 400°K this correction is 
entirely negligible. 

Computations for the hard sphere model of 
Eq. (10) were first made. Equation (13) corre- 
sponds to a diameter o of 2.99A, which is con- 
siderably greater than the value 2.76A derived 
for the O—O distance in ice from x-ray measure- 
ments,!® or the diameter 2.53A determined from 
viscosity measurements with the aid of the 
Sutherland formula.'? With Eq. (10), however, 
it was not found possible to represent the data 
adequately with such low values of o, the slope 
of B against 1/7 thus obtained being always 
too great. Satisfactory representation of the data, 
as shown in Table II, was secured with Eq. (10), 
by taking o=3.16A and c= 11.03 K10-® erg-cm*. 
The application of Eq. (5) with s=24 leads to 
smaller diameters, more in accord with the 
x-ray and viscosity data. In Table III it is 


TABLE III. Second virial coefficients of HO, in cc/mole, 
inverse-power model, Eq. (5), with s=24. o=2.76A, 
c=7.04X 10-*9 erg-cm,® p= 1.83 X 10“ e.s.u. 








Belass AB 
— 352.6 
— 232.7 
— 167.0 
— 126.5 
—99.2 
—74.6 
— 54.3 


T°(K) 


Beale 


— 349.5 
— 230.5 
— 165.4 
—125.4 
—98.3 
—74.0 
—53.9 


Bobs Opss.-CALc. 


346.9 +26 
312 07 
-166.4 —1.0 
—1%i —O7 
-990 —-0.7 
-745 —05 
—543 —-04 


450 
500 
550 
600 
6665 
750 


| 
++t++4+4+4| 
Sees te Ge 
PRHORAW | 








shown that the experimental second virial coeffi- 
cients are adequately represented with the values 
o=2.76A, c=7.04X10-® erg-cm®, and p»=1.83 
X10-"* e.s.u., the latter being the figure selected 
by Dorsey’ as the most reliable. 


16 W. H. Bragg, Proc. Phys. Soc. London 34, 98 (1922); 
W. H. Barnes, Proc. Roy. Soc. A125, 670 (1929). 

17 Unpublished computations by Professor F. G. Keyes. 

18N. E. Dorsey, The Properties of Ordinary Water- 
Substance (Reinhold Publishing Corporation, 1940), p. 48. 
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There are, unfortunately, no other polar gases 
for which sufficiently extensive and precise data 
are available. 


DISCUSSION 


The term —cr~® in the potential energy of 
Eq. (1) was taken to represent the entire dis- 
persion and induction energy, for which theory 
indicates an expression of the form 


u= — (Cy7~® + Cor 8 +6379 + - + +) 
— (2u?ar-6-+- ++). 


The last term is the induction energy averaged 
over all orientations, » being the permanent 
dipole moment and & the mean static polar- 
izability, and the quantities ¢1, C2, cz; are likewise 
averaged over all orientations. The terms in C2 
and cs; are small, so that their effect can be 
adequately estimated by considering only the 
first term in the result. Accordingly we write 


is) «@ 
-- f r—‘dr= [ ur*dr, 
oC fo 


where uw is the quantity defined by Eq. 
This integration leads to the result 


6201(14+3¢2/5¢,0?+3¢3/7¢104) + 2y2a. 


(14) 


(15) 


(14). 


(16) 


Since Margenau" has given equations by means 
of which it is possible to estimate ¢2/c, and ¢3/¢1, 
it is easy to transform the derived values of c 
into derived values of c;, which can then be 
compared directly to theoretical estimates of 
this quantity. 

The equations given by Margenau may be 
written as 


C2/¢1=5(&ao/f)}, (17) 


where do is the radius of the first Bohr orbit in 
hydrogen and f is the so-called oscillator strength 
of the single characteristic frequency which best 
reproduces the dispersion of the refractive index 
of the gas. For ammonia we take &= 2.24 10-*4 
cm’, f=2.72, which are the values selected by 
Margenau.” Then the constant c= 1.057 X 10-*8, 
which was found to represent the data with 
Eq. (10), leads to the value 7.8X10-® erg-cm® 
for c;. A discussion of the theory of c; will not be 
undertaken here, since this quantity has already 


€3/Cy = 105 Gao/8f, 


19H, Margenau, Rev. Mod. Phys. 11, 1 (1939). 
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been considered for these molecules. Margenau’s 
most recent estimate” is 7.0X10-® erg-cm'®, 
while London” gives 9.3 X 10-® erg-cm®. 

For steam, Dorsey”! recommends &=1.59 
X10-** cm’, and the value f=2.41 is given by 
Wolf and Herzfeld.” From the value c=7.04 
X 10-* erg-cm® for the preferred model of Eq. (5) 
with s=24, the equivalent value of c; is found 
to be 4.7X10-* erg-cm® while for the hard 
sphere model of Eq. (10) it is 7.9X10-*. The 
theoretical value given by London” is 4.7 KX 10-®. 

The good agreement between the theoretical 
and derived values of c; suggests that, for a 
calculation of the second virial coefficient, the 
electrical anisotropy of many molecules may be 
sufficiently well expressed entirely in terms of the 
dipole moment. Since, however, the B values are 
very sensitive to the molecular diameter, which 
was here chosen to give the best fit of the data, 
it may be that the derived potential energies 
are somewhat fictitious. For the particular case 
of steam, it is certain that the actual hindrance 
of rotation is far greater than that derived from 
the point-dipole model, so that the quantum 
correction given by Eq. (11) may be much too 
small. It might be possible to shed some light 
on this question by going to the other extreme 
and treating the second virial coefficient of 
steam as the association constant for a rather 
rigid hydrogen-bonded dimer. 

More generally, the point-dipole model should 
be least successful for such substances as the 
alcohols and the carboxylic acids, in which the 
dipole moment is large and not centrally situated. 
For a treatment of dense gases or condensed 
phases it is also certain to be inadequate. There 
is, for example, no property in the potential of 
Eq. (1) which can lead to the tetrahedrally co- 
ordinated structure of ice. On the other hand, 
this model has been shown to give a much better 
approximation for steam than was generally 
believed, and it should certainly prove to be 
useful for the semi-empirical treatment of gaseous 
mixtures containing polar components. 

The author is indebted to Professor George 
Scatchard for helpful criticism. 


* F. London, Trans. Faraday Soc. 33, 8 (1937). 

aN. E. Dorsey, reference 18, p. 49. 

2K. L. Wolf and K. F. Herzfeld, Handbuch der Physik, 
Vol. 20, p. 493. 
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A critical study of the equilibrium and molecular data for phosphorus and some of its com- 
pounds has been made and is described. The data are used to calculate the entropies and free 


energies of the substances considered. 





INTRODUCTION 


HE importance of phosphorus compounds 
in many branches of chemistry makes it 
desirable to have reliable values for their thermo- 
dynamic constants over a range of temperatures. 
Only a few reversible reactions are known whose 
equilibrium states have been determined, but 
there are a number of phosphorus compounds for 
which the molecular structures and energy levels 
have been established with sufficient certainty to 
permit their use for the calculation of thermo- 
dynamic quantities. It is the purpose of the 
present paper to present and discuss the free 
energies, entropies and heat contents (enthalpies) 
calculated from equilibrium, molecular and ther- 
mal data for a number of phosphorus compounds. 
As a result of correlating the commonly used 
thermodynamic quantities with the analogous 
quantities arising in the theories of statistical 
mechanics,! the following relations, which apply 
to one mole of a perfect, non-degenerate gas, 
are obtained; these relations form the basis for 
the derived formulas that were used in making 
the calculations presented below. 


E=E,+RT°0 log Q/aT, 
S=(E-—E,)/T+R log Q—R log No+R, (1) 
F=E,—RT(log Q—log Nv), 


where Ep is the arbitrarily chosen energy of the 
system at 0°K, F is the Gibbs-Lewis free energy, 
S is the entropy, No is Avogadro’s number, and 


Q=)>° exp (—e,/RT), 
n=1 


é, being the energy of a single molecule in mth 
quantum state; the summation is taken over all 


Tolman, Statistical Mechanics (Oxford University 


Press, 1938). 


accessible quantum states including translation. 
From these relations it follows at once that 


H°—H=RT(Ta log Q/aT +1), 


where the zero superscript refers to the value of 
H when the gas is in the standard state chosen. 
The free energy equation in terms of Q is also 
readily derived from Eqs. (1) and is 


(F°—H,°)/T = R(log No—log Q). (2) 


For an isothermal reaction between substances 
maintained in their standard states one obtains, 


AH, =AF°— RTA(log No—log Q) 
= —RT log K—RTA(log No—log Q), (3) 


and AH,’ may be evaluated if the equilibrium 
constants K and the Q’s are known. AH,’ is 
independent of 7, and its constancy for a series 
of temperatures is a test of the consistency of 
the experimental K and Q values; the test is not 
always sufficiently sensitive, but it is quite 
useful in general. 

The equilibrium and molecular data for phos- 
phorus compounds which are available for use in 
connection with Eq. (3) are those for phosphorus 
vapor, phosphine and phosphorus pentachloride. 
Molecular data alone are available for a number 
of volatile phosphorus compounds, and may be 
used to calculate (F°—H,°)/T for a moderate 
range of temperatures. It is unfortunate that 
data are not available for the calculation of the 
thermodynamic constants for the three well- 
known oxides of phosphorus, namely, P.O«, 
P3046, and P.O 40. 

The fundamental constants used in the calcu- 
lations were taken from the International Critical 
Tables; this was done in order to assure con- 
sistency with the thermodynamic quantities pub- 
lished by other investigators. 
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TABLE I, AH values calculated from equilibrium 
and molecular data. 


P2(g) = 2P(g) 








1473 
3.859 
50,500 


1073 
10.027 
49,400 | 


| 
1173 | 1273 | 1373 
8.061 | 6.452 | 5.065 


49,600 | 49,900 | 50,200 


ae Katmos , 
— AH)’, cal. | 








PHOSPHORUS VAPOR 


Preuner and Brockméller? investigated the 
equilibria in phosphorus vapor in the tempera- 
ture range 800 to 1200° and concluded that the 
molecules P;,, P2 and P were present together at 
measurable pressures. Later, Stock, Gibson and 
Stamm’ made a similar study, and their results 
showed that the partial pressures of monatomic 
phosphorus were much too small to be observed 
even at 1200°. It will be shown by the use of 
Eq. (3) that the conclusions of Stock, Gibson 
and Stamm are correct, and that some interfering 
substances must have been present in the experi- 
ments of Preuner and Brockmédller. 

The emission band spectrum of P:2 has been 
carefully analyzed by Herzberg and by Ashley, 
and from their results we have derived the 
expression w= (n+4)780.76—(n+3)?2.98 cm 
for the vibrational energy levels; the moment of 
inertia of Pz is 90.47X10-* g cm®. The normal 
state of the molecule is '2 while that of the atoms 
into which it dissociates is 4S, the dissociation 
energy being 5.008 electron volts (115,450 cal./ 
mole). In Table I are presented the AH,° values 
calculated by means of Eq. (3) from the data of 
Preuner and Brockméiller. It is evident that AH,” 
is neither satisfactorily constant nor is it in 
agreement with the spectroscopically determined 
value, 115,450 cal./mole. 

Table II contains the values of Katmos calcu- 
lated from spectroscopic data alone; it is be- 
lieved that these equilibrium constants are much 
more reliable than those shown in Table I. It is 
evident that the dissociation of diatomic phos- 
phorus is extremely slight even at 1200°. 

Further evidence for the superiority of Stock, 
Gibson and Stamm’s equilibrium data appears 
when the reaction P,(g)=2P2(g) is considered. 

2 Preuner and Brockmiller, Zeits. f. physik. Chemie 81, 
129 (1913). 

3 Stock, Gibson and Stamm, Ber. 45, 3527 (1913). 


4 Herzberg, Ann. d. Physik 15, 677 (1932). 
5 Ashley, Phys. Rev. 44, 919 (1933). 
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The P, molecule has a tetrahedral structure as 
was shown by the electron diffraction experi- 
ments of Maxwell, Hendricks and Mosely ;*° the 
P—P distance is 2.21A. Venkateswaran’ photo- 
graphed the Raman spectrum of liquid phos- 
phorus and found three lines corresponding to 
the three allowed fundamentals for a tetratomic, 
tetrahedral molecule. The frequencies and, in 
parentheses, their degeneracies are 362(2), 463(3) 
and 607(1) cm~!. In Table III are presented the 


TABLE II. Values of Katmos calculated from 
spectroscopic data. 


P2(g) = 2P(g), Katmos = Pp*/Pp, 
AH>° = 115,450 cal. 








| 1273 
5400 


37.47 


1373 1473 
17 X 104) 36 10° 
33.99 | 30.96 


1073 | 1173 
0.74 | 91 
| 41.54 | 


T°K 
Katmos x 107° 
—In Bi otros 46.35 


| | 








experimental values of Katmos obtained by both 
Preuner and Brockmdller (PB), and Stock, 
Gibson and Stamm (SGS) together with the 
corresponding AH,° values for the reaction 
P,(g) =2P.(g). The constancy of the AH,” values 
obtained using the experimental results of Stock 
and co-workers is very satisfactory and is to be 
contrasted with the large and unidirectional 
variations shown by the AH,° calculated from 
Preuner and Brockmédller’s data. 

The results of our calculations for both of the 
reactions considered in this section show quite 
definitely that the measurements of Stock, 
Gibson and Stamm are free from any major 
errors, and further that dissociation of diatomic 
phosphorus is very slight even at 1200°. 

From existing data it is now possible to cal- 
culate a fairly reliable value for the entropy 
of white phosphorus, P(w). Macrae and Van 
Voorhis® studied the vapor pressures of both 
white and liquid phosphorus, and from their 
results we may write 


4P(w) = P,(g), A F°o93 = 5700 cal., 
AF ogg = 13,200 cal. ’ 


no correction has been made for gas imper- 
fections. The following entropies for elementary 


6 Maxwell, Hendricks and Mosley, J. Chem. Phys. 3, 
699 (1935). 

7 Venkateswaran, Proc. Ind. Acad. 3, 260 (1935). 
any and Van Voorhis, J. Am. Chem. Soc. 43, 547 
1921 
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phosphorus were obtained from molecular data 
and the vapor pressure result. 

P(w) 
10.1 


P2(g) 
52.15 


Ps(g) 
66.92 


P(g) 


S298 cal./deg. 38.99 


Entropy determinations of the four more com- 
mon modifications of solid phosphorus, namely 
P(w, a), P(w, 8), P(red) and P(black),® would 
add considerably to our understanding of this 
important element. 


PHOSPHINE 


The only reliable measurements of the equi- 
librium constants for the reaction 


P,(g) +6H:=4PH,(g) 


are those of Ipatiew and Frost." They worked at 
total pressures of from 18 to 227 atmospheres 
and over the temperature range 627 to 771°K. 
The equilibrium concentrations of phosphine 
were determined by chemical analysis, and the 
partial pressures of P.(g) were assumed to be 
those in equilibrium with the red (violet) phos- 
phorus present as a solid phase. They used the 
vapor pressure measurements of Smits and 
Bokhorst" in this connection. Since the vapor 
pressures of red (violet) phosphorus do not 
appear to be accurately reproducible, the equi- 
librium measurements of Ipatiew and Frost may 
be subject to errors on this account. 


TABLE III. The values of AHo® from molecular and 
equilibrium data. 
P,(g) = 2P2(g), Katmos= Pp,?/Pp, 


tae 4 
Kstmos (PB) 
AH)®, cal. | 
Katmos (SGS) — {0.00645 .0375 .1671 

AH; - | 53,540 | 53,590 $3,670 | 








1173! 1273 | 1373 | 
0329! .1053| .2763 | 


1473 
.6579 


| 1073 | 
0.00855 | 
48,460 | 49,730| 51,000 | 52,310 | 53,440 


6118 
53,660 








The infra-red and Raman spectrum of phos- 
phine has been the subject of a number of 
investigations. We have chosen the frequency 
assignment of Sutherland” for our calculations. 
Only one moment of inertia, that about the axis 


* Jacobs, J. Chem. Phys. 5, 945 (1937). 

10 Ipatiew and Frost, Ber. 63, 1104 (1930). 
“sien and Bokhorst, Zeits. f. physik. Chemie 91, 249 

See Hibben, The Raman Effect (Reinhold Publishing 
Corporation, 1939); Yost and Anderson, J. Chem. Phys. 
2, 624 (1934). 
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perpendicular to the symmetry axis of the 
pyramidal molecule, is obtainable from the infra- 
red bands. Stevenson" has shown recently that 
the bond angle H—P—H is probably 93°. Table 
IV presents the calculated values of A/7,°. All of 
Ipatiew and Frost’s data are not presented, but 
only enough to show the representative varia- 
tions in AH,”. 

The variations in A/J,° are considerable. 
AH,®=1530 cal. has been chosen as the most 
acceptable average, and it has been used to 
calculate the values of K and In K shown in 
Table V. Greater precision in the measurements 
of this equilibrium is without doubt possible, 
and no other method seems as adequate for the 
estimation of AH since a direct determination 
does not appear feasible ; indirect methods would 
unavoidably involve the large heats of formation 
of phosphorus and hydrogen halides or oxides. 
The values in Table V are consistent with the 
entropy measurements made by Giauque and 
Stephenson“ on phosphine. 


PHOSPHORUS PENTACHLORIDE 


Measurements of the dissociation constants 
for the reaction 


PCI5(g) = PCls(g) +Clo(g) 


have been made by Holland and Nernst," and 
by Fischer and Jiibermann.'® In neither case 


were the measurements made with great pre- 

TABLE IV. The dissociation of phosphine. w,(1)||=2320, 
wo(1)|!=970, w3(2) 1 =2430, w(1) 1 =1118 cm. =], 
=6.22X 10~ g cm?, Is=7.00X 10-* g cm’. 


4 6 
Katmos = Ppu,/Pp,Pu. 








T°K 
627 


—In Katmos Ho’cal. = T°K 
25.08 
25.60 
25.31 
25.51 
25.90 
26.30 
26.00 
27.24 
26.85 


—In Katmos 
28.41 


28.16 
28.83 





1,110 771 
1,756 
1,398 
860 
1,381 
1,923 
1,516 
2,398 
1,835 


677 








18 Stevenson, J. Chem. Phys. 8, 285 (1940). 
ma and Giauque, J. Chem. Phys. 5, 149 
1 . 
be _ and Nernst, Zeits. f. Electrochemie 18, 234 
1912). 

Fischer and Jiibermann, Zeits. f. Anorg. u. angew. 
Chemie 235, 340 (1938). 
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TABLE V. The dissociation of phosphine. AH, = 1530 cal.; 
A Fog = — 18,370 cal. 


P.(g) +6H2(g) = 4PH3(g) 


AND D. M. YOST 


TABLE VI. The dissociation of phosphorus 
pentachloride vapor. 


PCl3(g) +Clo(g) = PCls(g) 








T°K 627 677 731 


Katmos X 10” 10 5.6 te 
—In Katmo 25.3 25.9 26.7 
P(w) -+3/2H2(g) = PH:;(g), A Fog, = —3,140 cal. 











cision, but the results are in good agreement with 
each other. Electron diffraction experiments by 
Schomaker and Sargent!’ show the PCI; molecule 
to be a trigonal bipyramid with an axial P—Cl 
distance of 2.12A and an equatorial P—Cl dis- 
tance of 2.03A. Moureau, Magat and Witroff'® 
have investigated the Raman spectrum and find 
for w1, w2(2), w3(2), ws, w5(2) and we(2) the values 
100, 190, 271, 392, 449, and 495 cm™, respec- 
tively. w; and ws are not known; they are per- 
mitted in the infra-red but not in the Raman 
spectrum. 

In principle it is possible to evaluate w; and ws 
from the equilibrium data. To do this con- 
veniently the function 7B(T) for the reaction 
was defined by the relation 


F°— Hy’ 


=AH),"—RT In Q7Qs, 


where Q; and Qs are partition functions of the 
form (1—e~'*/+7)—! for the frequencies w; and ws, 
respectively. Since 


—} (—) 
A{ ——— } =( ——— 
( z T PCl; 


, F°—H,° (—=) 
—RIn mill ainieemene a pee 
Or0s ( i ) PCI, T Cle 


and the right-hand side may be evaluated from 
the known energy states of PCl;, PCls and Cla, it 
follows that TB(T) may be evaluated from 
A((F°—H,°)/T)’ and the equilibrium data. The 
following four sets of values for w;, ws were 
assumed 


@7 200 250 450 600 cm— 
Ws 500 450 600 600 cm=! 
17 Schomaker and Sargent, private communication. 


18 Moureau, Magat and Witroff, Comptes rendus 205, 
276 (1937). 


Katmos —AH0° T°K Katmos —AH0 
122.63 21.44 | 502.2 1.752 21.35 
88.82 21.33 | 439 37.18 21.25 
67.96 21.35 | 443 30.33 21.27 
45.22 21.61 | 462 15.80 21.62 
28.53 21.34 | 485 4.09 21.42 
24.09 21.24 | 495 3.843 21.34 
15.52 21.22 | 534 0.4767 21.35 
13.42 21.26 | 556 0.2048 21.31 
8.88 21.20 | 574 0.1074 21.25 
8.00 21.40 | 613 0.0242 20.92 
5.816 21.20 | 623 0.0234 21.21 
2.175 21.54 | 631 0.0186 21.21 
Mean 21.32 














TABLE VII. Internuclear distances in phosphorus 
compounds, in A. 








MOLECULE P-O, P-S ZX-P-X 


P, 

P, — Tetrahedral 
PH; ‘ — 93° 
PF; 7 — 104 








PCI; 101 
PBrz 100 
Cl;PO f 5 106 
C1;PS . 9. 107 
Cl;PCl, 1 Trigonal 
bipyramid 


CONAUWP PwWre | 2 
| 








1 Ashley, reference 5 of text. 
2 Maxwell, Hendricks and Mosley, reference 6 of text. 
3 Stevenson, reference 13 of text. 
4 Brockway and Wall, J. Am. Chem. Soc. 56, 2373 (1934). 
a9 eee Hampson, Sutton and Jones, Trans. Faraday Soc. 33, 852 
6 Brockway and Beach, J. Am. Chem. Soc. 60, = (1938). 
7 Beach and Stevenson, J. Chem. Phys. 6, 75 (1938). 
8 Schomaker and Sargent, reference 17 of text. 


and for each set, 7B(T)+RT In Q;,Q3;=AH,° was 
calculated at the temperatures of the equilibrium 
experiments. It was found that with w;=450 cm“! 
and ws=600 cm, AH, did not drift with 
temperature as Table VI shows. 


FREE ENERGY TABLES 


No reliable equilibrium data are available for 
the calculation of the free energies and heat 
contents of phosphorus compounds other than 
those described above. Fortunately there exist 
adequate molecular data and heats of formation 
for several of the important halogen-containing 
compounds. In Tables VI—-X are assembled both 
the molecular and energy values for these 
compounds. 
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THERMODYNAMIC PROPERTIES OF PHOSPHORUS 


TABLE VIII. Frequencies, in cm™, degeneracies and symmetry numbers of phosphorus compounds. 








SYMMETRY 
|| 61|| ol 2| NUMBER 


770 — — — 2 

372(2) 463(3) 607 — 12 

2320 970 2430 1118 

890 531 840 486 

PCls $11 257 480 ‘190 

PBrs 400 162 380 116 

Cl;PO 1290 267 581 193 
Cl;PS 748 245 539 

Degeneracy 1 1 2 2 





ty ¢ 
NNn 








! Ashley, reference 5 of text. 

2 Venkateswaran, reference 7 of text. 

3 Hibben and Yost and Anderson, reference 12 of text. 

4 Cabannes and Rousset, Ann. de physique 19, 229 (1933). 
5 Simon and Schulze, Naturwiss. 25, 669 (1937). 


TABLE IX. Heat content differences, H—Hp°, in kcal./mole for phosphorus compounds, 








Ps PH; PF; PCls PCls PBrs 


3.37 2.42 3.84 5.22 4.25 
4.23 2.90 3.65 4.74 6.62 5.21 
5.09 3.39 4.39 5.65 8.03 6.16 
6.89 4.45 6.00 7.50 10.91 8.06 
8.74 5.61 7.71 9.40 9.98 

10.62 6.88 9.50 11.32 11.93 

12.53 8.22 13.26 13.89 

14.46 9.66 15.21 

16.41 11.16 17.16 

18.37 12.72 

20.32 14.36 

22.26 15.97 

24.24 17.67 

26.21 19.40 
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TABLE X. Free energies, —((F°—Ho°)/T), and entropies of gaseous phosphorus compounds, in cal./mole deg. 








P Pe Ps PH; PF; PCls PBrs PCI; Cl:PO 


42.10 54.37 61.66 68.89 66.76 63.45 
43.42 55.99 63.75 71.28 69.53 65.77 
44.54 57.41 65.61 73.31 71.97 67.85 
46.47 59.98 68.87 76.82 76.30 71.57 
48.13 62.24 71.68 79.72 80.08 74.74 
49.61 64.29 74.09 - 82.32 83.41 77.69 
50.95 66.13 76.30 84.64 86.36 80.32 
52.19 67.84 78.31 89.14 82.60 
53.34 69.37 80.07 91.56 84.76 
54.43 
55.45 
56.41 
57.34 
58.24 


50.20 











The values of (F—H,°)/T and H-—H, in’ connection with known thermochemical'® and 
Tables IX and X may be used in connection with published entropy values to calculate the A F293 
similar tables for other substances when the free of formation shown in Table XI. These AF%a 
energy of a given reaction is desired. Table X values refer to the elements in their commonly 
contains also the standard entropies at 25° of ~ 19 Bichowsky and Rossini, Thermochemistry (Reinhold 
the compounds listed, and these were used in Publishing Corporation, 1936). 
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TABLE XI. Standard free energies of formation of some 
gaseous phosphorus compounds at 25°. 








PCls 
74,800 


POCI; 


127,300 


PH; 
3140 


PCls 
65,300 


CoMPOUND 


—AFo9g, cal./mole 











known standard states. The well-defined a white 
phosphorus is taken as the standard state for 
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this element notwithstanding the fact that the 
red and black modifications are the more stable 
at room temperature. 

The free energy and entropy values given in 
this paper are to replace those calculated several 
years ago by Anderson and Yost.” 


20 Anderson and Yost, J. Chem. Phys. 4, 529 (1936). 
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Electric Polarization of Carboxylic Acids in Dilute Solutions of Nonpolar Solvents 


I. The Relation of Electric Polarization to the Association of Carboxylic 
Acids in Hydrocarbon Solvents 


HERBERT A. PouL,* Marcus E. Hosss Anp PAuL M. Gross 
Department of Chemistry, Duke University, Durham, North Carolina 


(Received February 10, 1941) 


The association and polar character of various carboxylic acids have been studied in the 
solvents benzene and heptane by means of measurements of the electric polarization in very 
dilute solutions. The data have been analyzed by assuming an equilibrium between dimers and 
monomers of the acid molecules. A rather large anomalous contribution to the polarization of 
the dimer has been found. Possible origins of this polarization have been briefly considered. 


W§ may consider, as a simplified picture, 
that concentrated solutions of the car- 
boxylic acids in hydrocarbons consist of a mixture 
of large molecular aggregates of the acid sur- 
rounded by the solvent. With dilution these 
aggregates break up until a state is reached 
where the dimer type of molecule predominates. 
This species persists over a large concentration 
range, and only at relatively low concentrations 
does considerable dimer dissociation take place 
leading to the appearance of appreciable fractions 
of the single species. 

At present little is known of the extent and 
character of the dimer dissociations and still less 
is known of the polar character of the acid 
molecules which are present at low concentra- 
tions. Determinations of the dielectric constant 
and densities of these dilute solutions provide, 
through the polarizations derived from these 
values, a simple and effective means of studying 

* Part of thesis of Herbert A. Pohl submitted in partial 


fulfilment of the requirements for the Ph.D. degree in 
Chemistry at Duke University. 


the molecular structure of these acids.! This has 
been carried out for a number of carboxylic 
acids in benzene and heptane with a view to 
gaining further information about the following: 
(1) the factors governing the dimer-monomer 
dissociation equilibrium, (2) the effects produced 
by different groups attached to the carboxyl 
group and by changing solvents, (3) the values 
of the molar polarizations and apparent moments 
of the single and double species of acid molecules 
and (4) the evidence for the existence of an 
appreciable ‘‘atom polarization”’ in the dimer. 


EXPERIMENTAL 


Dielectric constant differences between the 
solvent and the solution were measured with the 
heterodyne beat, parallel capacitance method. 
The electrical circuit has been described else- 
where.” The condenser containing the solutions 
_ 1 The first use of polarization data in studying association 
in solution was by Hojendahl, ‘‘Studies of dipole-moment,”’ 
Thesis, Copenhagen (1928), p. 33. 


2 Hobbs, Jacokes and Gross, Rev. Sci. Inst. 11, 126 
(1940). 
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to be measured consists of three concentric 
monel cylinders, the second insulated from the 
others by quartz rods. The condenser has a 
capacitance in air of approximately 450 uuf and 
requires approximately 400 ml of liquid. Details 
of construction have been given previously.’ 
Densities were determined to five parts per 
million with a pycnometer similar to that used 
by Connell.* The temperature during the density 
and dielectric constant measurements was main- 
tained constant to better than +0.002°C by 
means of a water bath. All measurements were 
made at 30.0°C. Solutions were made up by 
weighing out stock solutions and adding, in the 
absence of moisture and oxygen, small amounts 
of the solution to the solvent. The stock solutions 
were made and the necessary transfers effected 
with the aid of a ‘‘dry box,’ a special weight 
burette, a transfer apparatus, and the use of 
dry Ne to transfer the stock solution to the 
solvent. 


PURIFICATION OF MATERIALS 


Benzene 


Kahlbaum benzene which had been dried over 
P.O; was refluxed over sodium for at least ten 
hours and then distilled through a 180-cm 
Dufton column. All except a small low boiling 
fraction distilled within 0.01° range. The benzene 
was further purified by fractional crystallization 
until a freezing range of 0.01°C or less was 
obtained. Before use in a measurement, the 
benzene was refluxed over sodium for at least 
two hours and then, after two hundred cc were 
distilled to rinse moisture from the distilling 
column, the remainder was distilled directly into 
the transfer flask. 

The boiling point found was 80.10°, d,*° 
= 0.86835. Timmermans’ gives b.p. = 80.20°, d,° 
= 0.86844. Wojciechowski® gives b.p.=80.093° 
+0.002° on benzene from four sources. Cohen 
and Buij’ give b.p.=80.2° and d,*° =0.86844. 


_* Pohl, Hobbs and Gross, ‘‘Dielectrics,” Ann. New 
York Acad. Sci. 40, 389 (1940). 
1933) Vosburgh and Butler, J. Chem. Soc. 933 


aodgytimmermans and F. Martin, J. Chim. Phys. 23, 733 
* Wojciechowski, Rocaniki Chem. 16, 534 (1936). 
7E. Cohen and J. S. Buij, Zeits. f. physik. Chemie B35, 
270 (1937). 


ACIDS 


Heptane 


The heptane used had been previously washed 
with sulphuric acid, sodium hydroxide, and then 
fractionally distilled from sodium. This material 
was twice redistilled from sodium before use in 
the present work. It boiled over a 0.04°C range. 
This was refluxed for six hours over sodium 
before distilling it into the transfer flask for the 
measurement, d,°°= 0.67475. 


Acetic acid 


Two liters of C.P. acid were found to freeze at 
16.0°. This was fractionally crystallized until the 
freezing point was 16.4°, and was then frac- 
tionally distilled in a Widmer column. The 
freezing point of 16.63°+0.05° was unchanged 
by further crystallization. Timmermans® gives 
16.55°C as the freezing point. 


Formic acid 


A sample of Kahlbaum formic acid was 
fractionally crystallized once and then distilled 
in an all-glass Widmer still. B.p. 100.2°, uncor- 
rected, range 0.15°C. This acid was fractionally 
crystallized six times and was found to melt at 
8.5°C+0.1 over a range of 0.07°C. Timmermans® 
gives 8.40°C. The critical solution temperature 
in anhydrous benzene according to the method 
of Ewins® was determined and found to be 
73.4°+0.3°C thus indicating that the present 
sample probably contained less than 0.02 percent 
water. 


Propionic acid 


Fractional distillation is regarded by Timmer- 
mans® as a satisfactory method of purification. 
A sample of Kahlbaum acid was fractionally 
distilled in a Widmer column four times. The 
final sample of 45 grams boiled over a range of 
0.07°C. B.p.=141.5+0.2°; Timmermans® gives 
141.35°C. A rough determination of the freezing 
point gave a value of —22°C+1°C. Timmermans 
gives —20.8°C. 


Trimethyl acetic acid 


A sample of Eastman Kodak material was 
fractionally distilled, b.p. 164.15+0.20°C. After 


8M. J. Timmermans and Mme. Hennaut-Roland, J. 
Chim. Phys. 27, 401 (1930). 
° A. J. Ewins, J. Chem. Soc. 105, 530 (1914). 
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TABLE I. Chlorobenzene in benzene (30°C). 











fe X10 AE X108 P2 | f2X108 AE X108 P2 
0.0683 0.236 81.6] 1.089 3.829 82.5 
0.2408 0.836 81.8] 2.140 7.517 82.4 
0.3408 1.189 82.1} 2.642 9.280 82.5 
0.3730 1.284 81.4| 5.810 20.470 82.6 
0.5192 1.810 82.1| 8.300 28.870 82.0 
0.6606 2.31 82.4 ; 
0.7756 2.710 ewes go oe 
0.8690 3.020 81.9 | fo PRE= 











TABLE II. Formic acid in benzene (30°C). 











fo X108 AE X108 P2 | feX108 AE X108 Ps 
0.1121 0.423 68.8 | 2.286 4.045 39.0 
0.1601 0.479 57.2} 6.392 8.376 32.2 


0.4902 1.280 51.5 | 8.464 10.294 30.7 
0.6981 1.673 48.4 | 12.965 14.129 28.8 





1.634 3.161 41.5 | 15.103 15.931 28.3 
Ad 
=- =0.0946 
fe 








TABLE III. Formic acid in heptane (30°C). 








t2X108 AE X108 P2 f2X108 AE X108 P2 


0.2877 0.091 20.9 | 3.368 0.677 17.5 
0.3235 0.123 22.8 | 3.459 0.639 17.1 
0.8806 0.240 19.6 | 9.162 1.572 16.7 
0.9900 0.218 18.1 | 15.515 2.812 17.0 
2.233 0.454 17.6 
Ad 
— =0.0 
hi 820 














distillation it was fractionally crystallized three 
times, redistilled and the crystals dried over 
P.O, for twenty-four hours. 


Monochloroacetic acid 


An Eastman Kodak sample was fractionally 
crystallized three times and was then recrystal- 
lized six times from benzene. The crystals after 
each precipitation were washed with cold benzene 
and dried over 99-percent H2SO, in a vacuum 
desiccator in accordance with the method sug- 
gested by P. M. Gross.” M.p.=62.3+0.5°, 
range =0.2°. 


Benzoic acid 


A dry sample of standard benzoic acid from 
the National Bureau of Standards was used 
without further treatment. 


oar” M. Gross, Ph.D. dissertation, Columbia University, 
1919, 
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TABLE IV. Acetic acid in benzene (30°C). The numbers 
in parentheses at the left of this table indicate the four series 
of measurements made on this compound. This 1s included 
to show the agreement obtained between the different series. 














fe X108 AE X108 P2 S2X108 ASE X108 P» 


(1) 0.0548 0.150 58.7 | (2) 1.557 1.661 33.8 
(2) 0.1351 0.299 50.9| (3) 2.498 2.315 31.7 
(1) 0.2506 0.475 46.0| (1) 2.826 2.454 30.9 
(2) 0.2885 0.531 45.3 | (4) 2.931 2.502 30.6 
(2) 0.5894 0.865 39.9 | (3) 4.940 3.709 29.1 
(3) 0.7380 1.017 38.5/ (1) 7.946 5.146 27.6 
(1) 1.039 1.249 35.8) (1) 12.37 7.236 26.6 
(4) 1.541 1.634 33.7 | (3) 14.18 8.323 26.5 


Ad _ 0.0817 


Sa 














TABLE V. Acetic acid in heptane (30°C). 














fo X108 AE X108 P2 | f2X108 AE X108 P: 

0.1543 0.072 28.2 | 3.3053 0.911 22.8 

0.3771 0.136 25.2 | 5.544 1.497 22.6 

0.6091 0.198 24.2 6.218 1.650 225 

1.3903 0.418 23.5 | 11.302 2.917 re Be 

1.7120 0.499 25.2 | 11.531 2.998 22.3 
“ =0. 1164 


Butyric acid 

Kahlbaum acid was fractionally distilled three 
times in a Widmer column and then fractionally 
crystallized twice. The product of the second 
crystallization gave a melting point of —5.7°C 
with a range of 0.1°C. 


RESULTS 


Explanation of tables 


Tables I-X contain the data obtained by 
measurement of the dielectric constants, densi- 
ties and mole fractions of the various solutions. 
fz=mole fraction of the solute; AE=increment 
in the dielectric constant of the solution over 
that of the solvent; P2,=average polarization 
of the solute per molecular weight Mz of solute 
stated ; Ad/fo=ratio of the increment of. density 
to mole fraction. 

For convenience of reference the constants 
for the solvents which were used in calculating 
the polarizations are listed here: Mpenzene = 78.05: 
Mreptane= 100.19; d, benzene=0.86835; d,* 
heptane =0.67475; Erenzene (30°) = 2.2627; 
Eucptene (30°) = 1.910. 
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TABLE VI. Propionic acid in benzene (30°C). 








AE X103 P2 f2 X10 ASE X108 
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TABLE VIII. Trimethyl acetic acid in benzene (30°C). 


f2 X108 AE X10° 


P2 | fe X103 AE X108 





0.237 
0.249 
0.660 
0.694 
1.333 


56.7 
56.1 
45.9 
44.4 
38.3 





1.302 
3.474 
3.883 
6.908 
7.695 


Ad _ 0.0781 


Se 


1.371 
2.616 
2.631 
4.245 
4.281 








TABLE VII. Butyric acid in benzene (30°C). 








AE X108 


Pe 


t2X108 


AE X108 





0.188 
0.193 
0.331 
0.508 
0.856 
0.923 
1.221 


76.4 
61.2 
60.2 
51.8 
50.8 
47.4 
47.6 


Ad 





1.063 
1.517 
3.223 
5.627 
8.510 
10.330 
10.565 


—- =0.0906 


Se 


1.198 
1.588 
2.534 
3.977 
5.000 
5.798 
6.428 








The last value is an average of the value given 


by Smyth" and that of deBruyne.” 


EXPERIMENTAL ERROR 


0.0376 
0.1082 
0.4162 
0.6252 
1.118 


S2X108 


0.090 
0.205 
0.482 
0.604 
0.899 


TABLE IX. Benzoic acid in benzene (30°C). 


AE X108 


72.4 
61.8 


| 2.102 


te X108 


1.411 
1.746 
2.027 
3.195 


AE X108 





0.1154 
0.1627 
0.3516 
0.5825 
1.342 


0.361 
0.465 
0.856 
1.266 
2.473 


2.368 
| 4.357 


| 5.482 


8.096 


9.727 


3357 


3.919 
6.528 
7.967 
11.134 


12.150 








fo X108 


AE X108 


P2 | feX103 


AE X108 


TABLE X. Monochloroacetic acid in benzene (30°C). 


P2 





0.1051 
0.1358 
0.6193 
0.6328 


0.665 
0.841 
3.120 
3.189 


114.3 | 
112.2 | 
95.1 
95.2 


1.884 
4.729 
5.248 
10.066 


7.970 
17.516 
19.167 
33.984 


83.1 
75.1 
74.4 
70.1 


The validity and precision of the results given 
are affected by the possibility that, in these very 
dilute solutions, errors due to adsorption, im- 
purity, and normal experimental errors are 
larger than those usually encountered. An 
examination of the scattering of the experimental 
points justifies confidence in the data with regard 
to normal random error. The magnitude of 
constant errors such as those due to conductivity 
of the solutions, evaporation losses from the 
solution, purity of the chemicals used, inductance 
corrections, nonlinearity of the condenser with 
dielectric constant, introduction of moisture in 
handling, etc., have been carefully minimized. 
The magnitude of the effect of adsorption is very 
difficult to estimate, but it is felt that adsorption 
does not play an important role at concentrations 
above 5 X 10-5 mole fraction, since measurements 
made on chlorobenzene in benzene show no 
significant deviations above this concentration. 
A complete analysis of the errors of the dielectric 


1 Smyth and Walls, J. Am. Chem. Soc. 54, 1859 (1932). 
® deBruyne, Davis and Gross, J. Am. Chem. Soc. 55, 
3936 (1933). 


1.862 7.901 83.3 | 11.095 37.196 


Ad 
=0.3846 
fr 








constant-density method has been detailed 


elsewhere.® 
DISCUSSION 


There can be conceived several possible 
explanations for the curvature of the carboxylic 
acid P2 vs. fz curves. The most likely one appears 
to be the following. In dilute solutions there 
exists a mixture of highly polar, single acid 
molecules and slightly polar, or nonpolar double 
molecules. As the total acid concentration is 
increased, association of single to double mole- 
cules produces a less polar mixture of acid 
molecules. The fall of the P»2 vs. fe curves is 
interpretable on this basis. In order to analyze 
the results quantitatively on the basis of this 
postulated mechanism, we have assumed that 
in these very dilute solutions the activities could 
be replaced by concentrations and the molar 
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Fic. 1. Calculated curves for the acids. 


volume of the solution by that of the solvent. 
We have postulated an equilibrium of the form 
K=Cs?- V/Cp where Cs and Cp are, respectively, 
the moles of single and double molecules per 
unit volume, and V is the molar volume of the 
solvent. It has been shown elsewhere® that K 
can be expressed in terms of polarizations 
according to the following equations. 


" (P2—Pp/2)?-2fe 2K» 
(Ps—P.)(Ps—Pp/2) Ps—Pp/2 





(1) 


It is convenient to use the quantity Ko, 
which is defined by Eq. (1), instead of K in 
evaluating Ps and Pp as follows: 


P sKo=foP2?—feP2Pp+fe(Pp/2)*+P2Ko (2a) 





or 
Pp [Koe?+4f2Ko(Ps—Pp/2) ]?-—K 
PaeiiPt 0 2k o( L's D 0 (2b) 
2 2fs 
= Pp/2+Gs(Ps—Pp/2). ’ (2c) 


In these equations K is expressed as moles of 
solute per molar volume of solvent, P2 is the 
molar polarization of the solute based on the 
molecular weight of the singles, Ps is the molar 
polarization of the single molecules, Pp is the 
molar polarization of the double molecules, Gs 
is the fraction of the molecules in singles form, 
and f: is the mole fraction of solute on the basis 
of the singles molecular weight. 

In the previous paper* the constants K, Ps, 
and Pp of Eq. (1) were evaluated from the 
experimental data by the method of averages. 


However, because of the extreme sensitiveness 
of the values of K to slight errors in P, in the 
very dilute region, it has been found desirable 
to employ the following more sensitive and 
critical method of determining the three con- 
stants. The method of averages (previously 
employed) is first used to obtain the approximate 
values of Ps and Pp; then with these values of 
the constants, K is calculated for each concen- 
tration by substituting the experimental P»2 and 
f2 values in Eq. (1). If the values employed for 
these constants (Ps and Pp) are those which 
most closely represent the experimental data, 
the computed values of K should exhibit random 
variation around a mean value which is constant 
over the entire concentration range. In general, 
however, the value of the equilibrium constant 
showed some variation with concentration with 
this preliminary choice of Ps and Pp values. 
It is evident from Eq. (1) that the values of K 
are quite sensitive to the choice of P's in the very 
dilute regions where P2 approaches Ps, and to 
the choice of Pp/2 in the more concentrated 
regions where P2 is but little greater than Pp/2. 
By a method of successive approximation, the 
values of Ps and Pp/2 are adjusted to give the 
most constant value of K over the entire con- 
centration range. While the two methods of 
solution often give practically identical values 
of the constants, the values found by the new 
method gave closer agreement with the experi- 
mental points. A comparison of the polarizations 
calculated by the new method with the values 
determined experimentally is given in Table XI 
for the case of acetic acid in benzene. A com- 
parison of this table with Table 13 of our 
previous paper*® will show that the constants 
found by the new method fit the experimental 
data more closely than those determined by the 
method of averages. 


TABLE XI. Acetic acid in benzene. 











Fz X108 Peale Pexp f2 X108 Peale Pexp 
0.0548 58.5 58.7 1.557 33.8 33.8 
0.1351 51.3 50.9 | 2.498 31.6 31.7 
0.2506 46.1 46.0 | 2.826 31.1 30.9 
0.2885 45.1 45.3 | 2.931 30.9 30.6 
0.5894 39.7 39.8 | 4.940 29.1 29.1 
0.7380 38.1 38.5 | 7.946 27.7 27.6 
1.039 36.0 35.8 | 12.37 26.6 26.6 
1.541 33.9 33.7 | 14.18 26.3 26.5 
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ACIDS 


TABLE XII. Constants of acids. 








Acid Ps Pp 


PS(B) Ps(A+0) P p(A+0) 





39.4 
31.8 
44.0 
43.2 
51.1 
62.2 
71.0 
90.8 
114.8 


Formic 

Formic* 

Acetic 

Acetic* 

Propionic 

Butyric 
Trimethylacetic 
Benzoic 
Monochloroacetic 


76 
(68)t 
71 
(69)+ 
76 
78 
86 


85 
126 


Om = WNOWES ss | RX 
Coto em WOROH- !| X 





8.3 

8.3 
13.0 
13.0 
17.6 
22.2 
26.8 
32.4 
17.7 


15.2 
18.0 
17.2 
15.9 
17.8 
17.4 
26.0 
79.4 








* Heptane measurements; all others in benzene. 


+ Equilibrium constant too small for accurate determination of Ps. The constants also are approximate. 


It is difficult to find an independent basis for 
appraisal of the real significance of this close 
agreement and its relation to our basic assump- 
tion as to the equilibrium in the solution. The 
equations may at least be considered as con- 
venient expressions which represent the observed 
data and which, at the same time, offer a 
reasonable quantitative interpretation of the 
process occurring in the dilute solutions. 

In addition to this analytical treatment of the 
experimental data, the results were also analyzed 
by several graphical methods with substantially 
the same result. 

Table XII gives the results of the analytical 
evaluation of the constants Ps, Pp and K 
(computed from Ko) for the various acids. The 
other data listed in this table are the electron 
polarizations Pz), the atomic plus orientation 
polarizations of the single molecules Ps:4+o), 
and of the double molecules Pp,4+0), as well as 


TABLE 


calculated dipole moments, ws and yup, of these 
molecules. 

A comparison of the values of the moments 
found with those reported by other investigators 
is given in Table XIII. As may be seen, we have, 
for purposes of this comparison, called Pp .4+o) 
a permanent dipole polarizability and have thus 
calculated up. A consideration of this point will 
be given later in this paper. 

If we consider first our ws values for acetic 
and propionic acid in relation to those of Zahn, 
it is seen that within the limits of legitimate 
comparison of solution and gas measurements 
the values are in good agreement. The value of 
us for formic acid seems to be anomalous in 
that it is larger than the gas value. That formic 
acid may be a special case is further indicated 
by the fact that Wenzke, using dioxane as the 
solvent, found the very high value of 2.05D. 

Our values for the apparent moment of the 
various doubles molecules differ considerably 


XIII. 











SuTTON*2 


(Gas) BRIEGLEB?® 


PRESENT WORK 


WENZKE* BENZENE HEPTANE 





Formic 
MS 
MD 

Acetic 


1.15 


1.02 

Benzoic 
Ls 

MD 0.79 


2.07 1.82 
1.06 
1.74 1.68 
0.94 
1.75 1.68 
0.88 


1.61 
1.13 


1.75 








+ These values of ug and up were calculated from Ps and Pp values given by Coop, Davidson and Sutton by using » =0.0127 [(P —Pp)T]}}?. 
1C. T. Zahn, Phys. Rev. 37, 1516 (1931). 

* Coop, Davidson and Sutton, J. Chem. Phys. 6, 905 (1938). 
*G. Briegleb, Zeits. f. physik. Chemie B10, 205 (1930). 

‘ Wilson and Wenzke, J. Chem. Phys. 2, 546 (1934). 

°K. L. Wolf, Physik. Zeits. 31, 227 (1930). 
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from those previously recorded, but are, however, 
rather consistent among themselves. All previous 
solution measurements have been made in such 
ranges of concentration as to make it very 
probable that sensible amounts of either mono- 
mers or of higher polymers than the doubles 
were present. Because of this, it is felt that the 
lack of agreement is in large measure due to the 
fact that no correction for these contributions 
to the total polarization has heretofore been 
made. It therefore seems justifiable to regard 
our values for the dimer polarizations as the 
most reliable ones at present available. 

A further examination of the data of Table XII 
reveals the following generalities. First, the 
values of the singles molecules, exclusive of 
chloroacetic, are essentially constant—-particu- 
larly in those cases in which the values are 
considered most accurate. This is in agreement 
with expectation since the carboxyl group is the 
origin of the moment. It may be noted that the 
values of the singles are in general considered 
not precise enough to show the small differences 
in moment that might exist when groups which 
are but slightly different are attached to the 
terminal carboxy] group. 

The existence of a rather large and sub- 
stantially constant value for Pp a+o) for all 
molecules examined (exclusive of chloroacetic) 
is also revealed by Table XII. This polarization 
might originate from any or all of the following: 
(a) a thermal bending of the two halves of the 
dimer so that a net permanent dipole polarization 
results; (b) a solvent effect such as suggested 
by Frank; (c) a permanent lack of electrical 


13F, C. Frank, Proc. Roy. Soc. (London) A152, 171 
(1935). 








symmetry of the acid dimer ring; (d) an atomic 
type of polarization; (e) the presence of polar 
aggregates not considered in the equilibrium 
which was assumed in the calculations; (f) a 
stepwise breaking of the dimer ring. 

The possibility that (a) is the origin of this 
polarization seems to be excluded by the fact 
that the net carboxyl group moment is nearly 
at right angles to the line joining the two 
carbons of the two carboxyl groups. Any bending 
perpendicular to the ring such as folding of the 
two molecules about a line joining the two 
hydrogen bonds would probably contribute very 
slightly, if at all, to the polarization. The 
suggestion in (b) is unlikely since the values of 
Mp found are very large for such solvent effects 
and because essentially the same values are found 
for acetic acid in benzene and heptane. The 
difference in values for formic acid in these two 
solvents is regarded as an anomaly due to the 
formic acid and is not thought to be a difference 
due to the usual type solvent effect. 

Of the remaining suggestions (d) appears to 
be the most likely explanation of such a large 
value of Pp(a,o). The basis for this belief will 
be considered in more detail in the following 
paper. 

Another observation that may be made on 
examining Table XII is the low dissociation 
constants for acetic and formic acid in heptane 
as compared to the values in benzene solutions 
(Kvenz/Knept).!° This indicates that the unsatu- 
rated and more highly anisotropic benzene 
molecule is a very much better dissociating agent 
than the heptane molecule. 


14 C, S. Brooks and M. E. Hobbs, J. Am. Chem. Soc. 
62, 285 (1940). 
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Electric Polarization of Carboxylic Acids 


II. Polarization of Heavy Acetic Acid and of the Three Fluorobenzoic Acids in Benzene and the 
Anomalous Polarization of the Dimer 


ARTHUR A. MAryott,* MArcus E. Hosss anp PauL M. Gross 
Department of Chemistry, Duke University, Durham, North Carolina 


(Received February 10, 1941) 


The electric polarizations of acetic acid-d and o- m- and p-fluorobenzoic acids have been 
determined in benzene solution at concentrations ranging from approximately 10-° to 10-2 
mole fraction of solute. On the assumption of a monomer-dimer equilibrium, the polarization 
of the monomer and of the dimer molecules and the equilibrium constant have been calcu- 
lated. The unexpectedly high polarization of the dimer is best interpreted on the basis of the 
assumption of an atomic polarization of about 20 cc. 


N the previous paper of this series! the polar- 
izations of a number of carboxylic acids were 
determined as a function of the concentration 
in dilute solutions of benzene and heptane. The 
variation of polarization with concentration was 
satisfactorily interpreted on the assumption that 
there existed only a monomer-dimer equilibrium 
in the concentration range studied. It was noted 
in particular that the acid dimer had a polariza- 
tion from 16 cc to 25 cc greater than one might 
expect from a molecule with the dimer structure.” 
This additional polarization might have for its 
origin a permanent electrical dissymmetry in all 
or a fraction of the polymerized molecules or an 
abnormally high atomic polarization. With a 
view to throwing further light on this question, 
measurements have been made on benzene solu- 
tions of heavy acetic acid and 0-, m-, and p-fluoro- 
benzoic acids. The concentration range covered 
was from approximately 10-° to 10 mole 
fraction of solute. 


MATERIALS 


Acetic acid-d was prepared in this laboratory 
by Mr. Henry Linschitz from acetyl chloride and 
99.6 percent D.O. The details of the preparation 
and physical constants will appear in a separate 
note. M.p. 15.66+.05°C (range 0.2°). 

0-Fluorobenzoic acid was made by oxidizing 
purified o-fluorotoluene in the manner described 


* Ethyl-Dow Chemical Company Post Doctorate Fellow, 
Duke University, 1940-1941. 

Pohl, Hobbs and Gross, J. Chem. Phys., this issue. 
a and Brockway, Proc. Nat. Acad. Sci. 20, 336 


in Organic Syntheses.’ The product was purified 
by recrystallization several times from benzene. 
F.p. = 123.4-123.6°C. 

m-Fluorobenzoic acid was obtained by oxidation 
of pure m-fluorotoluene and purified by recrystal- 
lization from benzene. F.p. = 122.5-—122.8°C. 

p-Fluorobenzoic acid was prepared from purified 
p-fluorotoluene in the same manner as were the 
ortho- and meta- acids. It was recrystallized 
twice from toluene. F.p. = 180.5-180.7°C. 


EXPERIMENTAL 


In the measurements which follow, the appa- 
ratus was modified in one respect; oil replaced 
the water in the thermostat. This was done 
primarily to reduce corrosion of the exposed 
metal parts of the condenser and its supports. 
A series of determinations of the polarization of 
acetic acid in benzene made after this change 
agreed, well within experimental error, with 
those obtained with the water thermostat. The 
experimental procedure, calculations, and method 
of analyzing the results have been described 
previously.! In the present case, because of the 
very limited solubility of the fluorobenzoic acids 
in benzene, the more concentrated solutions were 
made by adding the solid acid directly to the 
solution in the measuring condenser. To prevent 
the solid from falling down into the condenser 
where it would be difficult to determine when it 
had completely dissolved, a cloth plug was 


3 Organic Syntheses (John Wiley and Sons, 1930), Vol. X, 
20. 
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inserted in the lower constricted part of the 
receiving bulb until the acid had dissolved. 

The experimental data are given in Table I in 
which f2 is the mole fraction of solute, P2, its 
polarization calculated on the basis of the 
monomer molecular weight, AC, the capacity 
change of the condenser in arbitrary units after 
introduction of the given amount of solute. The 
vacuum capacity of the condenser was equal to 
912.5 of the arbitrary units or approximately 
450 uyf. The ratio, Ad/f2, of the change in density 
upon addition of the solute to the mole fraction 
of solute was assumed to be linear with con- 
centration’ throughout the range studied and 
only the densities of the pure solvent and of the 
most concentrated solution were determined. All 
measurements were carried out at 30.0°C. A meas- 
urement involved the determination of the di- 
electric constant of a solution of the given 
substance at four or five different concentrations. 
The o- and m-fluorobenzoic acids were measured 
twice and p-fluorobenzoic acid and the heavy 


TABLE JI, 








Je X108 AC P2 f2X108 AC 





Heavy acetic acid in m-Fluorobenzoic acid in 
benzene at 30°C benzene at 30°C 


0.0435 0.115 0.0487 0.217 1 
0.2667 0.483 .0 | 0.0828 0.353 1 
0.4883 0.718 5 | 0.2420 0.935 
0.7198 0.904 1 | 0.4359 1.548 
0.8960 1.117 .9| 1.007 3.307 

1.771 1.785 1} 1.369 4.288 
2.011 1.862 8} 2.851 8.678 

2.122 1.948 32.7 | 3.269 9.815 
5.795 3.888 .7| 4.179 12.50 

5.934 4.049 9 | 5.974 17.58 

9.239 5.434 4| 7.047 20.56 

9.521 5.554 Ad/f.=0.498 


14.67 7.799 oe eT 
-Fluorobenzoic acid in 
16.01 8.257 benzene at 30°C 


Ad/fs=0.0912 (av.) 0.0650 0.175 


o-Fluorobensoic acid in 0.0717 0.192 
benzene at 30°C 0.0793 0.214 
0.0902 0.598 139.7 | 0.2062 0.417 
0.2913 1.626 125.2 | 0.2920 0.534 
0.4654 2.470 118.5 | 0.2942 0.516 
0.6636 3.269 112.1 | 0.5110 0.837 
1.354 6.065 104.7 | 0.6042 0.913 
1.771 7.567 101.5 | 0.8782 1.258 
4.081 15.814 95.0} 1.164 1.483 
4.150 15.929 94.4 | 1.178 1.503 
7.361 26.68 91.0} 2.289 2.613 
9.118 $2.21 89.4 | 3.068 3.334 
12.64 43.28 87.6 Ad/f2=0.510 
Ad/f2=0.513 (av.) 
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* Schulz, Zeits. f. physik. Chemie B40, 151 (1938). 


O-FLUOROBENZOIC ACID 


Fic. 1. Calcu- 
lated curves 
(solid lines) and 
experimental 
points for the 
acids. 


M-FLUOROBENZOIC ACID 


P-FLUOROBENZOIC ACID 


© ACETIC ACID 
* ACETIC ACID-D 





4 10 


6 8 
Sf, X 10° 


acetic acid three times, the separate measure- 
ments agreeing within experimental error. 

Figure 1 shows the variation of polarization 
with mole fraction for the various acids. 


DISCUSSION OF RESULTS 


In Table II are listed some additional data cal- 
culated from the experimental results in Table I. 
The symbols have the same significance as 
previously indicated.! The dipole moment in 
Debye units has been calculated from the equa- 
tion, 4=0.221(P(4+0))?. 

The results for the fluorobenzoic acids seem to 
indicate that the carboxylic acid dimer structure 
is symmetrical and that the source of the 
apparent moment of the dimer in all probability 
lies in a rather large atomic polarization. The 
Pp:axo) term for p-fluorobenzoic acid is no 
greater (actually less) than that for benzoic acid. 
This indicates that the C—F bond moments in 


this dimer are antiparallel. If the Pp a+o) ° 


term were largely due to some electrically un- 
symmetrical configuration such as might arise 
from slow thermal bending of the halves of the 
dimer or association to unsymmetrical higher 
polymers, a larger value of Pp aio) would be 
expected for the p-fluorobenzoic acid dimer than 
for the benzoic acid dimer because of contribu- 

tions from the C—F bond moments. 
In the light of this interpretation, the data fcr 
ortho- and meta-fluorobenzoic acid dimers can 
) 
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ELECTRIC 


POLARIZATION OF CARBOXYLIC 


ACIDS 


TABLE Il . 








K X104 


PS(B) Ps(A+0) 





43.8 
44.0 
155.0 
150.4 
85.2 


Heavy acetic 3.0 68 
Acetic* , 71 
o-Fluorobenzoic : 158 
m-Fluorobenzoic 136 
p-Fluorobenzoic 106 


13.0 55 
13.0 58 
32.4 126 
32.4 104 
32.4 74 








* Data on ordinary acetic acid in benzene taken from the preceding paper for purposes of comparison. 


also be best accounted for by assuming a high 
atomic polarization together with an additional 
contribution to the polarization from the fluorine 
substituents in the ring. The fluorine contribution 
might originate either from free rotation of the 
phenyl groups or from the effect of an equal 
mixture of cis- and trans-forms.* The two forms 
for the ortho-acid are represented below. The 
two pheny! groups are considered to be coplanar. 
It seems quite probable that free rotation of the 


O---—-H—O F 
eee. pv . ees 
»—C CY S 
/ ~~ 4 
O—H---O 


68 2 08 


cis 
5  0---H—O 

oe ar 
\ i WA —_— A 


O—H ---O 
“trans” 


phenyl groups is inhibited by resonance although 
from the point of view of dipole calculations both 
assumptions are equivalent. If the C—F bond 
moment is taken as 1.46D, one calculates an 
orientation polarization contribution for the two 
C—F bonds of 66 cc when the fluorine atoms 
are in either the ortho- or meta-position. Sub- 
tracting the sum of this value and the electronic 
polarization of 65 cc from the total dimer 
polarization leaves a residual term of 24 cc for 
ortho-fluorobenzoic acid dimer and 19.5 cc for 
the dimer of meta-fluorobenzoic acid. These 
values are of about the same magnitude as the 
anomalous part of the polarization in the other 
acid dimers studied.! 

There is, however, an alternative explanation 
of the dimer polarization. In each case the 


* Brooks and Hobbs, J. Am. Chem. Soc. 62, 2851 (1940). 


polarization, Pp,:a,o), may be regarded as the 
sum of two terms: one a dipole contribution from 
the part of the molecule not associated with the 
carboxyl group ring and the other a relatively 
constant term associated with the ring. Two 
vectors at right angles would lead to such an 
additive relation; and, consequently, the possi- 
bility of a permanent moment of approximately 
1D associated with the carboxy] ring and directed 
parallel to the line joining the two pheny! groups 
is not excluded. However, it is difficult to see 
how an electrical dissymmetry of this sort could 
arise in the highly symmetrical ring of the dimer. 

A further possibility that has been suggested 
is that the rapid fall of the P2 values with in- 
creasing concentration is a consequence of hin- 
drance of rotation due to dipole interaction. 
Such interaction might conceivably account for 
the anomalous atomic polarizations we have 
postulated. However, as no general treatment of 
this interaction appears to be available in a form 
applicable to our data this possibility could not 
be tested as an explanation. 

Another possible explanation of this anomaly 
is that a fraction of the dimer molecules exists in 
a form in which one of the hydrogen bonds is 
broken. The data on p-fluorobenzoic acid indi- 
cate that the phenyl groups are in such positions, 
relative to each other, that the two C—F bonds 
are antiparallel. Because of the rather low energy 
of the hydrogen bond, as compared to an ordi- 
nary chemical bond one might reasonably expect 
considerable low frequency bending about this 
bond in the case that only one such bond con- 
nected the two halves of the dimer. If this 
bending did occur, it is highly probable that 
there would result a sensible increment in the 
value of the polarization associated with the 
dimer of the para-acid. This is not observed. 
A further piece of evidence against the partially 
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dissociated form is the practically constant 
value of the anomalous term. This constancy 
would require that the fraction of the dimer in 
the partially dissociated form remain essentially 
the same in all the aromatic and aliphatic cases 
studied. In view of the marked variation of the 
equilibrium constants this seems unlikely. Fur- 
thermore, the data on acetic acid previously re- 
ported show that the anomalous polarization is 
practically constant in going from benzene to 
heptane, although the tendency for the dimer to 
dissociate into single molecules is very much less 
in heptane than in benzene. 

In view of the above evidence it seems that 
only an atomic polarization term offers a reason- 
able interpretation of the anomaly. Large atomic 
polarizations have been observed by Sutton® for 
a number of compounds, particularly the metallic 
acetyl acetonates. In addition, Sutton,’ by deter- 
mining the polarization of formic acid vapor at 
various temperatures and using Coolidge’s® vapor 
density data to determine the polarization of the 
dimer, found that the dimer polarization is 
apparently independent of temperature, thus 
indicating that an orientation type of polariza- 
tion is not involved in the case of the dimer 
moment in that instance. 


6 Coop and Sutton, J. Chem. Soc. 1938, 1269. 

938) Davidson and Sutton, J. Chem. Phys. 6, 905 
(1938). 

8 Coolidge, J. Am. Chem. Soc. 50, 2166 (1928). 


HOBBS 


AND GROSS 

There exist no previous data on the moments of 
the fluorobenzoic acids. However, the moments 
for the monomers found here have the values to 
be expected in relation to the moments of some 
other halogenated benzoic acids in dioxane.® If 
the C—F and COOH group moments are taken 
as 1.45 and 1.64D, respectively, and the angle 
between the ¢—C axis and the COOH group 
moment is taken as 76°, one calculates 2.4, 2.2 
and 1.9D for the moments of 0-, m- and p-fluoro- 
benzoic acids, respectively. The experimental 
values 2.5, 2.2, and 1.9D agree well with those 
calculated. 

Comparison of the data on heavy and light 
acetic acids in Table II and in Fig. 1 shows that 
there is little difference in the Ps and Pp values 
of the two acids. The atomic polarization of the 
heavy acid dimer is 17.8 cc and not significantly 
different from the value 18.0 cc found for 
ordinary acetic acid. 

The equilibrium constant for o-fluorobenzoic 
acid is considerably higher than those for the 
other isomers. This is probably due to a stabiliza- 
tion of the o-fluorobenzoic acid monomer through 
hydrogen bond formation with the adjacent 
fluorine atom. A more complete consideration of 
the equilibrium constants will be presented at a 
later date, when additional data including meas- 
urements at different temperatures and constants 
for other related acids have been obtained. 
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Thermodynamic treatment shows that for chemical 
reactions whose equilibrium constants go through maxima 
(or minima) the relationship 


In K—In K,,=p(T—0)?+q(T—6)'+--- (a) 


holds in the neighborhood of 7,,=6. K,» is the extreme 


value of K. The constants and g are equal to 


(ACp) aACp\ 1 
so > =o / A ” _— 1/ ( ) —ry 
are? 2— — '(2/3)(ACpo— (1/6) \ 7 » Re 


Equation (a), aside from the cubic and higher terms, has 
been proposed by Harned and Embree to apply to the 
experimental data on acidic dissociation. It is shown that 
for acidic dissociation in aqueous solution the properties 


I. INTRODUCTION 


T has been known for some time that the 

dissociation constants for certain acids and 
bases pass through maxima near room tempera- 
ture.! This phenomenon and the general temper- 
ature dependence of acid and basic dissociation 
have been discussed several times.?-> In this 
paper we shall develop from thermodynamics 
the theoretical equation to give the temperature 
dependence of dissociation and then apply a 
model for the process to estimate certain con- 
stants. Our equation will be compared with the 
empirical formula which has been suggested. 

The thermodynamic treatment of Part II, 
although simple, has not, to our knowledge, 
been given before. For the mechanism of the 
dissociation process of Part III we follow pretty 
close to the theory which has been expressed, 
at least implicitly, a number of times, especially 
by Gurney* and Baughan.* However, we believe 
that our treatment is more detailed than theirs 
and may be more useful. 


II. THERMODYNAMIC TREATMENT 
The dissociation constant K for any acid can 


* Assistant Professor, Kyoto Imperial University, Japan. 
use and Embree, J. Am. Chem. Soc. 56, 1050 

* Walde, J. Phys. Chem. 39, 478 (1935). 

* Gurney, J. Chem. Phys. 6, 499 (1938). 

* Baughan, J. Chem. Phys. 7, 951 (1939). 

* Harned and Owen, Chem. Rev. 25, 31 (1940). 


of water almost completely determine the values of p 
and q. Estimates of these constants are made. The mecha- 
nism of the dissociation process is discussed, and it is 
shown that the equation for the dissociation of water can 
be represented by 


(m+n+2)H,O—H;tO(nH.O)+OH-(mH,0),  (b) 


where m+n+2=16. This is in accord with the result of 
Harned and Fallon for mixtures of water and dioxane 


logio Po/p=12.2N:, (c) 


in which p is the ion product for the mixture which con- 
tains a mole fraction N»2 of dioxane; and po is the ion 
product for pure water. 


be written 
AH AS 


(AlT), (AS) 
In K—In Ke= ——+ viet 
RT 


R Ré R 


(1) 


Here Kg is the dissociation constant for some 
definite temperature 6. We now make the 
following expansions, 


p 


OA 
AH= (AH (AC,)o(T-8) +1(— ) (T—6)?, 
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r AC, 
AS— (as) f —dT 
o I 


r| 


T—6 T—6\? T—06\* 
GY 
6 6 6 


Substitution of these into (1) gives 
(AH)¢ 
In K—In Ke= (T-—6) 
Re? 
(AC;) 
; ‘\er-0 
2Re? 
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Fic. 1. Schematic representation of the energy changes 
accompanying the dissociation of HA. Only the main terms 
of Eq. (6) are shown. The ionic atmosphere term and dipolar 
terms are both small. 


The result is actually an infinite series in 
(T—0@), but we have retained only the first 
three terms. This amounts to the assumption 
that 0?AC,/dT?=0, which is quite justifiable. 
Up to this point the development is perfectly 
general. Now we shall consider the case of a 
maximum value for K. At the maximum point 
we have 


dInK AH 
--( ) =0; (AH)n=0. 
oT RT"? « 








If we take 0 as Tmax, Eq. (2) now becomes 








AC,) 
ee a 
dAC,\ | (T—4)' 
-{3 acne (= r), | Re 
The empirical equation 
In K—In K,= —p(T-—6)* (4) 


has been proposed by Harned and Embree! to 
apply to the dissociation of acids. They give 
p=1.15 X10~ as a general constant for all acids. 
Walde? has re-examined the experimental data 
and concluded that they do not actually fit a 
strictly quadratic equation. At this point we 
cannot say how large the cubic term in (3) is, but 
we expect that, in general, it should not be 
vanishingly small, especially for large values of 
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(T— 86). We shall now proceed to estimate the 
constants in Eq. (3), and to do this we need a 
more detailed model of the dissociation of acids. 


III. A MopEL FoR THE DISSOCIATION 
Process oF ACIDS 


The dissociation constant for an acid which 
dissociates in the manner 


HA+H,0=—H;0*+A-— (5) 
we shall give as 
— RT |In K=AF=AF (dissociation into H+A) 
+ E(ionization H) — E(electron affinity A) 


ey 1 1 87e7c\ } 
Gre“ ae) 
2\r+ r- DRT 
+AF’+AF (hydration ions) 
(=~ = 
Ta~ THA 2D+1 
Here D is the macroscopic dielectric constant 
of the solvent, e the unit charge; wa- is the 
dipole moment of the A~ ion; ¢ is the concen- 
tration of the H;+O ion; r+ and r— are the radii 
of the ions taken as suitable spheres. The other 
symbols have their usual significance. 

Figure 1 diagrams our imaginary process by 
which the acid becomes dissociated. Here the 
work done in the process is the ordinate and the 
H-—A distance is the abscissa. The molecule is 
first dissociated into radicals (AF (dissociation 
into H+A)). This process is carried out in such 
a way that there is no change in the energy of 
hydration. An electron is now transferred from 
H to A (E (ionization potential H) —£ (electron 
affinity A)). The energy difference in the solution 
between these two quantum states for the 
electron concerned should be practically the 
same as in a vacuum, since the only screening 
which matters for an electron is caused by 
electrons more close to the nucleus about which 
it moves: The presence of the solvent molecules 
should not alter this term. There is energy 
stored up in the dielectric due to the creation 
of the electric field and this is given by the 
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Born term® 


esi 1 1 

— —+—)(—-1). 

2D\r+ r—-/\X\D 

The ion atmosphere interaction is given by the 
next term. The H* ion is now bound to a water 
molecule to form H;+O. This energy and such 
free energy changes as the difference in resonance’ 
energy between the acid and the ion are included 
in AF’. The ions are now hydrated, AF (hydration 
ions). The last term appears when there is a 
difference in the dipole moment of the acid and 
negative ion.* For this small effect we have 
written only the first term of a series. This 
corresponds to treating the molecules as spherical 
and having radii rs— and ra with dipoles at 
their centers. 

Some of the terms in (6) are ‘‘electrostatic,”’ 
i.e., are concerned with the storing up of energy 
in the dielectric medium, and some are non- 
electrostatic. The hydration energy will be 
partially electrostatic and partially nonelectro- 
static. We shall assume that by proper choice 
of the ion radii we can represent the electrostatic 
part of the free energy as 


e? 1 1 
ar=—(—+—). (7) 
2D\r+ r—- 


Thus we write the free energy of dissociation as 


esi 1 
AF=6+— —+—), (8) 


2D\r+ r—- 
where ® is the sum of the dissociation energy of 
HA (into H+A) and the nonelectrostatic part 
of the hydration energy. Such factors as differ- 
ence in resonance energy between the ion and 
acid’ will be included partially in one term and 
partially in the other. We are considering only 
infinitely dilute solutions (no ion atmosphere) 
and neglecting dipole effects. 

This is the expression for AF which Baughan‘ 
used, although his model was not as detailed. 


® Born, Zeits. f. Physik 1, 45 (1920). 
_* Pauling, The Nature of the Chemical Bond (Cornell 
University Press, Ithaca, 1940). 
* Kirkwood, J. Chem. Phys. 2, 351 (1934). 


ACID STRENGTH 


We now get 


; asa 
AH =—RT? (—) =Als 
aT\RT 


e 1 1 d1n D 
E(tat)\(i4™), @ 
2D\r+ r- oT 


There is no good theoretical expression for the 
dielectric constant of water. There are a number 
of satisfactory empirical formulae available® and 
we choose the one given by Akerléf,! 


D=A exp [—B(T—T;) ]. (8) 


A is the dielectric constant at 7; and B is a 
constant. This fits the data well in the region in 
which we are interested, 


es i 1 
a= ate+—(—+—) —BT). (9) 
2D\r+ r- 
The value of B is 0.00472 and the second term 
of (9) is found to be negative for all temperatures 
above 215°K. 

Figure 2 gives the heat change in dissociation 
corresponding to AF of Fig. 1. It is interesting 
that the heat change of the electrostatic work is 
negative although the work itself is positive. 
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Fic. 2. Schematic representation of the heat changes 
taking place during the dissociation of HA. As in Fig. 1, 
only the largest terms are shown. 


®See Darsey, Properties of Ordinary Water-Substance 
(Reinhold Publishing Corporation, New York, 1940). 
10 Akerléf and Short. J. Am. Chem. Soc. 58, 1241 (1936) 
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Fic. 3. This is a plot of Eq. (14). 

x The heavy line (parabola) gives the 

Xr first term and the dashed line gives 

\\ the sum of the two terms. The points 

ra _| for the lower aliphatic acids are 

\e plotted. The agreement is rather 

\ — rough, as is to be expected, but there 

\s is apparently a cubic trend in the 
*\\. | experimental points. 
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It follows that the entropy change is large and 
negative. For the electrostatic term 


esi 1 
AH" =— —+—)(1-B7) (10a) 
2D\r+ r—- 
e? 1 1 
AS” = —— —+—)B. (10b) 
r+ r- 


Let us examine this situation further. We 
have given formulas which imply a sort of 
continuous dielectric medium and it may not be 
clear at first what the entropy of charging an 
ion is. On the other hand we have promised to 
include part of the hydration energy in the 
electrostatic term. It is seen from the formulas 
that the numerical result depends upon the 
rapidity of the variation of D with T, i.e., on B. 
This rapid variation of D is directly a result of 
the motion of the water dipoles increasing with 
temperature ;" thermal motion overcomes orien- 
tation caused by electric fields. However the 
randomness of the thermal motion increases 
much more rapidly for the water in solution 
than for the water which forms the hydration 
on the ions. The hydrated water is essentially 
“frozen” onto the ion. Thus as the temperature 
rises the entropy becomes more. negative; the 
heat also becomes more negative, passing 


11 The general features of the dielectric constant of water 
are quite well known. See, for example, Fowler and 
Guggenheim, Statistical Thermodynamics (Cambridge Uni- 
versity Press, London, 1939); Conference on Dielectrics, 
Ann. New York Acad. Sci. 40, 289-482 (1940) ; Kincaid and 
Eyring, J. Chem. Phys. 6, 620 (1938). 





through zero and gives K, a maximum for some 
temperature 0. 

It should be emphasized that the above 
treatment merely presents a model which makes 
easy the calculation of the constants in Eq. (3). 
Later we shall discuss further the dissociation 
process, but now we proceed with the derivation 
of the constants of our equation from this model. 


IV. COMPARISON WITH EMPIRICAL RELATION OF 
HARNED AND EMBREE 
Using Eq. (9) we find 

dAH es i 1 dAHs 
A -—-=-_(— -— B+. (11) 
oT r+ r- oT 

AHg is a nonelectrostatic term and should not 
be very temperature dependent. Baughan‘ has 
shown that the temperature dependence of a 
number of acids can be represented quite well 
by Eq. (9), when (1/r+ +1/r—) and AHs¢ 
are both constants. From his values of 
(e2/2)(1/r+ + 1/r—) for the lower aliphatic 
acids, we take 4.5105 calories as an average 
value. Thus we have a consistent theory, if we 
neglect dAHs/dT. Using Dog=78.5, we obtain 


AC, = —0.124 Te?.00s72(7-28), - (12) 


Substitution into (3) gives, for values of 6 near 298°, 


0.124 
In K—In K,,= ———-(T— 6)? 
2R0 


0.067 
+——(T—6)*. (13) 
Re? 





THEORY 


Taking the coefficient = 300, which is approxi- 
mately correct for the lower aliphatic acids, and 
several others, this becomes 


In K—In Kn = —1.03X10-*(T— 6)? 


+3.72X10-7(T—6)*. (14) 
Harned and Embree! give 1.15 X 10~ as the best 
value of p in Eq. (4) for a large number of acids. 
Walde’s analysis of each of the sets of data, 
however, indicates that there is a trend such as 
one would expect from a cubic term about the 
magnitude of that in (13). It is well to point out 
that we would predict a different value of the 
constants depending upon @. In some cases it is 
obviously impossible to neglect dipole effects, 
and this would also change the constants. 
Formic, acetic and propionic acids are examples 
in which the theory outlined above should apply ; 
@ is also near room temperature and measure- 
ments have been taken on both sides of Kn. 
In Fig. 3 we have plotted in the heavy (parabolic) 
curve the first term of (13). The dashed curve 
gives the result when the cubic term is added. 
The experimental points for formic,” acetic,” 
propionic’ and butyric! acids are also plotted. 
It is immediately obvious that the points 
actually fall nearer the cubic than the quadratic 
for two of them. Formic and butyric do not fall 
so well on either, but the trend of the points 
seems to indicate a cubic. 


V. Entropy or Acipic DISSOCIATION 


The fact that there is very nearly a reduced 
equation to explain the behavior of acids in 
aqueous solution indicates that the general 
nature of the process is determined by the 
properties of water. We point out now that the 
constants in Eq. (13) depend mostly on the 
temperature variation of the dielectric constant 
of water. Variations from a following of the 
general equation may be caused by specific 
properties of the acid molecules themselves. 


ass and Embree, J. Am. Chem. Soc. 56, 1042 
a Harned and Ehlers, J. Am. Chem. Soc. 55, 652 (1933). 
Harned and Ehlers, ]. Am. Chem. Soc. 55, 2383 (1933). 

“934 and Sutherland, J. Am. Chem. Soc. 56, 2039 
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STRENGTH 
From Eq. (10b) we have 


eri 1 4.90 X 10° 
AS” = -—(—+—)s BD) nes neemnennonemetonn 
2D\r+ r- 78.5 


0.00472 = — 29.4 


when the correct values for water‘ are substituted 
for T=25°C. This is the entropy decrease due 
to the charging of the H;O* and OH ions. The 
entropy decrease for the freezing of water (at 
0°C) is 5.3 entropy units. Thus we interpret the 
electrostatic entropy of ionization of water as 
the ‘freezing’ of the equivalent of 6 moles of 
water about the two moles of ions formed. The 
water around the ions is not held as tightly nor 
precisely in the same manner as molecules of ice; 
so as yet we do not know exactly how many 
molecules are involved in the ionization process. 

Now let us calculate the total entropy of 
ionization. From the heat of ionization and the 
equilibrium constant!® at 25°C for the equation 


2H,.O—H;*0+0-H (15) 
we obtain: 


AS = — 35.30. 


The difference between this and the part due 
to the electrostatic process is —5.90 units. This 
is about the entropy of freezing a mole of water, 
and it is easily understandable when we realize 
that we have not considered that the two ions 
H,;0* and OH-, are essentially in a frozen state. 
This last entropy of —5.9 units should be 
practically independent of the dielectric constant 
and temperature and thus appear in ®. The 
electrostatic part of the entropy depends upon 
how tightly the hydrated water is bound to the 
ions, but this latter entropy depends only upon 
the fact that the ions are hydrated. 

At this point, we cannot tell how much 
entropy is lost by a freezing of a mole of water 
due to the electrostatic effect. Consequently, 
we cannot guess exactly the number of moles of 
water involved in the ionization process (15). 
However, the number estimated as follows 
should not be far from the truth. We have 
interpreted the difference between the total 
and electrostatic parts of the entropy as due to 


16 Harned and Hamer, J. Am. Chem. Soc. 55, 2198 
(1933). 
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the freezing of the two water molecules which 
become the ions. Then, it can be concluded 
that for the ionization losing 35 units 12 moles 
of water will be involved. However, the number 
12 is probably a lower estimate, since the 
entropy loss due to the electrostatic effect should 
be smaller than 5.9 units. If we assume that by 
the freezing of a mole of water the reaction 
loses about half of the entropy which attends 
the actual formation of ice, i.e. 5.3/2, then, 14 
moles of water should take part in the ionization. 
We shall later present the evidence that this 
rough estimate, 12-14, is justified. 

It might be well to mention at this point that 
since we are taking 15 as the ionization process 
for water, our K’s are smaller than the con- 
ventional ones by the factor 1/55.5."" 

It is generally recognized that Eq. (15) 
represents the reaction better than the dissocia- 
tion into a proton and hydroxyl. This latter 
process, being a dissociation, should have a 
positive entropy change; Eq. (15) contains an 
equal number of molecules on either side and 
should have no entropy change of this sort. 

In Table I we have entropy changes calculated 
from experimental values for the dissociation 
constants and heats of ionization of.a number of 
acids. From the entropies for H;BO3; and H,PO.- 
it is apparent that their ionization processes 
are similar to that for HO. These molecules are 
expected to fit very well into the structure of 
water. The differences in the heats are probably 
due to weaker O—H bonds in their molecules. 

The aliphatic acids are distinctly different. 
If we were to assume that the acid HA were 
already ‘‘frozen’”’” more or less, then we expect 
an entropy change of about —29—3=—32. 
The —29 is for the electrostatic process (see 
above) and the —3 (about half of 5.9) the 
freezing of the water molecule which becomes 
H;0+. The entropies of Table I indicate that 
the experimental values are above this, which 
means that the HA molecule is already more 
restricted than we have guessed, i.e., that water 
has already “‘frozen’’ on it. This hypothesis 


17 The pK values in Table I were taken directly from 
>! eam indicated and are ‘‘conventional’’ ones, 1.e., for 
V, 
_(H*)(OH-) 
(H,0) 
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seems to agree with the facts: Formic acid, the 
most polar of the unsubstituted acids, has the 
smallest decrease in entropy, and as the hydro- 
carbon chain becomes longer the entropy be- 
comes more negative. Furthermore the chloro- 
and hydroxy-substituted acetic acids, which 
include very polar groups, are like formic acid 
rather than acetic. Entropy of this type should 
appear in ® since the energy of binding does 
not depend very strongly on variation of the 
dielectric constant. The water is attracted to 
the unionized molecules by dipole forces. Ac- 
cording to Kirkwood,*® the dependence of this 
energy on D is (D—1)/(2D+1) which is practi- 
cally equal to $ for all values of D in which we 
are interested. 

It is to be noted that the ionization constant 
of butyric acid is larger than that for propionic 
acid. This is due to an irregularity in the heat of 
ionization rather than to the entropy. The heat 
is actually —691 calories at 25° whereas for a 
regular decrease in K one might have expected, 
say —250 calories. The reason for this larger 
heat value is somewhat obscure but may well 
be due to a ring closure.'® Evidence for this will 
be discussed in a later publication. 


VI. FurTHER EVIDENCE CONCERNING 
Acipic DIsSsoOcIATION 


The conclusion that ions have a definite 
hydration about them suggests the writing of 











TABLE I. 
AciIp REFERENCE AH293 pKoss AS298 
H:,0 1 13,481 15.741 — 35.30 
H;BO; ii 3,360 9.236 — 39.77 
H2PO,- 2 731 7.206 — 38.80 
HCOOH 3 —13 3.754 — 25.44 
CH;COOH 4 —112 4.757 — 30.41 
C.H;COOH 5 — 168 4.874 — 31.06 
n-C3;H;COOH 6 —691 4.820 —32.45 
CICH:COOH 7 —1170 2.861 — 25.19 
OHCH:COOH 8 210 3.831 — 25.02 








1 See reference 16 of text. 

¥ Owen, J. Am. Chem. Soc. 56, 1695 (1934). 
2 Nims, J. Am. Chem. Soc. 55, 1946 (1933). 
3 See reference 12 of text. 

4 See reference 13 of text. 

5 See reference 14 of text. 

6 See reference 15 of text. 

7 Wright, J. Am. Chem. Soc. 56, 314 (1934). 
8 Nims, J. Am. Chem. Soc. 58, 987 (1936). 


18 Dippy, Chem. Rev. 25, 151 (1939). 








THEORY OF ACID STRENGTH 


TABLE II. 








MOLE 
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the dissociation process for water as 


(m+n+2)H20 
H;*O(mH,0)+OH-(nH.0). (16) 


The m molecules of water are ‘‘frozen”’ to the 
H;tO ion, and m to the OH- ion. We have 
pointed out that the entropy of the ionization 
process permits us to estimate roughly the 
number of moles of water involved. Due to the 
recent work of Harned and collaborators!*—*! on 
water dioxane mixtures we now have further 
evidence to consider. 

Let us first consider an ion in a water dioxane 
mixture. If the medium surrounding the ion 
remained at the same composition as the bulk 
solution, there would be a very large increase in 
the free energy of the charging of the ion over 
the charging energy in water since the dielectric 
constant is lower. This excess free energy would 
be given by (e?/2r)(1/D—1/Do), Do being the 
dielectric constant for water and D that of the 
mixture. The ion can lower its energy by sur- 
rounding itself with water molecules which it 
unmixes from the solution. This process will con- 
tinue until the most favorable condition has been 
reached. The work of removing x (=m+n+2) 
moles of water from a water-dioxane mixture of 
mole fractions N; and Nz is —xRT In yMi, 
where ¥ is the activity coefficient of water in the 
mixture. Thus the work done in transferring 
the hydrogen and hydroxyl ions from water to 
the mixture is given by: 


Asi 1 
=~akTInyWs-—(—-—). 17 
x} Do D 


'? Harned, J. Phys. Chem. 43, 275 (1939). 
1939)" and Fallon, J. Am. Chem. Soc. 61, 2374 
(1939), 


*t Harned and Fallon, J. Am. Chem. Soc. 61, 2377 (1939). 
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The hydrated ion has a definite radius, the same 
in both solutions; it is proportional to x’ and 
A is the constant factor to give the correct 
charging energy. This model is, of course, quite 
rough. It is well known” that simple Born terms 
do not actually give the charging energy in 
molecular media. The success of the Born term 
in the first part of this paper is attributed to the 
fact that the temperature dependence of the 
dielectric constant is the most important factor 
involved; the ion radii required for the Born 
charging term to fit experiment really carry in 
themselves a correction for the error in using a 
macroscopic dielectric constant where one should 
consider the molecular behavior of water dipoles. 

Minimizing w of equation (17) with respect to 
x we find that the extreme value of w comes 
when the electrostatic term is three times as 
large as the unmixing term. For a particular 
value of N,; the extreme value is 


Wmin = —XRT In yN,—3xRT In yMi. 


Thus if In yN; and (1/D)—1/D) have the same 
variation with mole fraction there is the same 
value of x for all mixtures and Eq. (16) takes on 
significance as representing the reaction. The 
experimental values for y calculated from vapor 
pressure data* indicate that In yN, is linear 
with JN, for the range for which the experimental 
data are available. Frequently (1/Do—1/D) for 
mixtures has been approximated as a linear 
function of Ne. For water dioxane mixtures it 
varies considerably from linearity in investigated 
ranges. However, if we assume that the effective 
dielectric constant varies linearly we can write* 


Po w 
logio —= =cNo. (18) 
2.3RT 


Harned and Fallon®**! have shown experi- 
mentally that w is a linear function of N2; they 
give c=12.2. 

We have seen what assumptions are involved 
obtaining agreement between our model and 
experiment for the free energy of transfer of 


2a Webb, J. Am. Chem. Soc. 48, 2589 (1926). 

% Everett and Coulson, Trans. Faraday Soc. 36, 633 
(1940). 

23 Hovarka, Schaefer and Dreisbach, J. Am. Chem. Soc. 
59, 2753 (1937). 

* See abstract. 





426 MAGEE, RI 


ions from an aqueous solution to mixtures. Our 
model also predicts the heat and entropy changes 
individually and we shall now calculate the heat 
term in w. We write 


w=2.3RT(3.05N2+9.15N2) 
= —14RT In yNi+9.15N2X2.3RT. 


The heat for the unmixing term we shall give as 


—14RT In y= rn (19) 
Here, we assume that the heat arising from the 
better bonding of the water molecules in water 
over that in the mixtures is the total cause of 
deviation of the mixtures from ideality. The 
heat for the ionic term is of the form given 
earlier in this paper, 


—9.15N2(1—BT) X2.3RT =bAH charging. (20) 

In Table II the experimental values of AH and 
64H are given. In the fourth column the values 
of 6AH are calculated by summing Eggs. (19) and 
(20). It is seen that 6AH is first positive and 
then negative as it should be, but it does not 
become negative rapidly enough. In the fifth 
column values are calculated taking the ratio of 
the unmixing to charging term as 1 : 2.5 instead 
of 1:3.0. The improvement is immediately 
apparent. The number of water molecules frozen 
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Fic. 4. Potential curves for the dissociation of HA in 
vacuum and in a continuous medium of dielectric con- 
stant D. 
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Fic. 5. Schematic potential curves for the two lowest 
states of an acid molecule (in an aqueous solution) as a 
function of the H—A distance. 


out, however, is increased to 16 for the 1 : 2.5 
ratio. 

The initial increase in heat is due to an 
increase in the charging energy. The negative 
heat term is due to the water-water bonding 
which is obtained in the water frozen on the ion. 
Thus it seems that our model for the process is 
in agreement in a rough quantitative way with the 
experimental data. 

The fact that the charging term cannot be 
taken as large as one might have expected (i.e., 
three times the unmixing term) is probably 
related to the structure of the ions. Hydrogen 
and hydroxyl ions fit into the water structure 
as well as the water molecules themselves. Other 
ions which spoil the structure in their hydration 
layers may have different behavior. 

The measurements on acetic acid”! in the water- 
dioxane mixtures indicate that the behavior is 
very similar to water. There seems to be less 
water involved, as one would have supposed. 
This similarity seems to indicate that a carboxy] 
group forms an ion which is quite similar to the 
hydroxy] ion itself. 

The behavior of other mixtures is not neces- 
sarily expected to be as simple as that of water- 
dioxane mixtures. Whether or not there will be 
a complete sorting of the medium by the ions 
depends upon the nature of the molecules in 
the medium. Ethyl alcohol, for example, could 
not be eliminated completely from the hydration 
layer because of its large dipole. Changing of 
the composition of the layer frozen to the ions 
would add to the complication of the problem. 
It seems, however, that it is possible to write 
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Eq. (16) for the ionization of water in water- 
dioxane mixtures with (m+n+2)=16. 


VII. THe DissociaATION PROCESS 


Let us now reconsider briefly the mechanism 
of the dissociation process. The problem should 
be treated by quantum mechanics. Although 
this is much too difficult to do, we do know the 
general nature of the solution. The general 
quantum-mechanical viewpoint of the subject of 
electrolytic dissociation has been discussed by 
Ogg and Polanyi.** In Fig. 4 a schematic set of 
curves is given. Curve A represents the dissocia- 
tion of the acid molecule in the gas phase into 
radicals. Curve B gives the polar state A~H*. 
If the molecule is put into a continuous medium 


= Ogg and Polanyi, Trans. Faraday Soc. 31, 604 (1935). 
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of dielectric constant D (there is no hydration) 
the only change is that the curve B is changed 
to that given by C. The dissociation energy of 
the polar state is decreased by the factor 1/D. 
Of course there will be a resonance between the 
polar and homopolar states and the result will 
be that the two solutions follow the dotted lines 
near the crossing. From this figure it is easy to 
see that the AH value for the electrostatic part 
of the ionization process (as of Fig. 2) is really 
negative. 

Introduction of solvent molecules complicates 
the process tremendously. The hydration energy 
of H+A is greater than that of HA. The water 
also rearranges in a complicated manner when 
the ions are changed. The net result may well 
be potential curves of the type shown in Fig. 5. 
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A method for calculating rates for a class of chemical reactions in vacuum is developed and 
applied to the reaction in which tungsten reduces magnesia to give free magnesium. The calcu- 
lations were tested experimentally and the observed and computed values were found to be in 


substantial agreement. 


EVERAL chemical reactions which would 

ordinarily be considered very improbable, 
are believed to occur at appreciable rates in 
vacuum. Some of these are of importance in the 
actual operation of electronic tubes. Typical 
examples are the reduction of thoria by tungsten 
in the thoriated tungsten filament and the re- 
duction of baria and other alkaline earth oxides 
in both the oxide-coated cathode and “‘batalum’”! 
getters. One of the most readily observed reac- 
tions of this type was discovered several years 
ago in the development of the early equipotential 
cathode vacuum tubes, in which magnesia was 
tested as an insulator between the tungsten 
heater wire and the cathode sleeve. Such insu- 


'E, A. Lederer and D. H. Wamsley, RCA Rev. 1, 117 
(1937); E. A. Lederer, RCA Rev. 4, 310 (1940). 


lators proved impractical because of a reaction 
in which the tungsten was etched until it burned 
out and a film of magnesium was deposited on 
the bulb. This reduction of the ceramic by 
tungsten is rather surprising, at first thought, 
since magnesia is such a stable refractory and 
since magnesium is usually a much more powerful 
reducing agent than tungsten. Since such reac- 
tions are important commercially, it was decided 
to investigate them in some detail, beginning 
with that between tungsten and magnesia, be- 
cause of its comparatively rapid rate. 
The reaction was thought to be 


2MgO(solid) + W (solid) 


—2Meg(gas)+W0O,(solid). (1) 


For the reaction to proceed as shown at any 
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pressure and temperature, the Gibbs free energy 
of the substances on the right must be less than 
that of those on the left. To determine the free 
energy difference at any temperature, it has 
seemed most convenient to compute as a function 
of temperature the free energy of each sub- 
stance involved and to combine the individual 
values algebraically. From data analyzed and 
published either in the International Critical 
Tables or by Kelley,’ it is possible to compute 
directly the free energy of all substances con- 
cerned in the above reaction except WO, and 
this can be obtained by applying well-known 
methods to published equilibrium measurements. 

The method applied to the first three con- 
stituents in (1) will be illustrated for the case of 
Mg and the computations required for WO, will 
be indicated. The summation will first be made 
for the standard conditions and then for those 
existing in vacuum tubes. From the latter result 
it is shown why the reaction goes as found and 
its maximum rate is calculated. Computed re- 
action rates were found to be roughly equal to 
those experimentally obtained. 


FREE ENERGY OF MAGNESIUM 


The Gibbs free energy of any substance is 
given by F=H—TS where, if no phase change 
occurs between 7; and 72, H and S are obtained 
from the equations 

T2 


Hr.,=Hr,+ C,dT 


T1 


2 K. K. Kelley, Bulletins No. 371, 383, 393 and 394, U.S. 
Bureau of Mines. 


E. 





MOORE 


and 
T2 Ce 
St2,=St+ —dT, 


T1 


where 7; is a temperature at which H and S are 
known. For solid magnesium, the specific heat 


is given by Kelley as C,=6.20+1.33X10°7' 


— 67,800/7? calories per gram in the temperature 
range 273°K to 923°K. The I.C.T. give the 
following values per mole at 298.1°K: H=5170 
joules; S°=34.0 joules per degree; F°= —4980 
joules. From these data, the free energy of 
magnesium can be computed at any temperature 
up to the fusion point. For the heat of fusion, 
9040 joules as given by Kelley was used, rather 
than the earlier value of 7300 joules given by 
I.C.T. The specific heat of liquid magnesium 
was assumed constant at 30.9 joules as suggested 
by Kelley, no other data being available. Using 
these values, H, S°®, and F® were computed up to 
the boiling point (1383°K). 

The heat of vaporization given by Kelley, 
136,000 joules, was used rather than the earlier 
value 262,000 joules given in the I.C.T. In com- 
puting H, S°, and F® above 1383°K, magnesium 
was assumed to form a perfect monatomic gas 
with C,=5R/2=20.8 joules per degree. To illus- 
trate the magnitudes involved, the values of //, 
S°, and F° are plotted up to 1900°K on Fig. 1. 


CALCULATION OF FREE ENERGY OF W AND MgO 


For magnesia and tungsten, the specific heat 
equations are given by Kelley and values of 
H, S°, and F® at 298.1°K are given in I.C.T. 
The free energies were determined up to 1900°K, 
as for Mg. 


SOLID Mg 


GASEOUS Mg 


Fic. 2. 
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REDUCTION 


CALCULATION OF FREE ENERGY OF WO, 


The best data for obtaining the free energy of 
WOz seem to be Chaudron’s* measurements of 
equilibrium pressures in the reaction: 


WO2.+2H22@W+2H,0 
in the range 600-1100°C. (2) 


These measurements are quoted by Smithells* 
who gives the following relation for the equi- 
librium constant as computed from them by 
Van Liempt, which is here assumed to be valid 
up to 1900°K. 


2222 


Prot 
= -——+1.690. (3) 


logio K,=logio [= 


He 


The free energy change in the reaction is 
given by 


—AF°=RT In K,. (4) 


The free energy of W over the temperature 
range was computed above and similar data are 
available? for Heand HO. Combining Chaudron’s 
measurements of the free energy change with the 
free energies of the constituents of the reaction 
(2), the free energy of WO. was computed as a 
function of temperature. 


COMPUTATIONS OF FREE ENERGY 
CHANGE IN REACTION (1) 


The work discussed above gives the free 
energy over the temperature range for all four 
substances and the change due to the reaction can 
be computed. The result of this computation is 


shown in Fig. 2, which shows that if the reaction: 


(1) proceeds, the free energy would increase at 
all temperatures. At first thought, this might 
suggest that the reaction shown could not pro- 
ceed spontaneously. 

However, the conditions in a vacuum tube are 
such that the reaction can proceed. The above 
computations give the standard free energy 
change AF® which refers to each material in its 
standard state, i.e., below 923°K magnesium is 
assumed to be solid, from 923 to 1383°K to be 
liquid and above 1383°K to be a gas at one 


* Chaudron, Comptes rendus 170, 1056 (1920). 
*Smithells, Tungsten, A Treatise on its Metallurgy, 
HH eae and Applications (Van Nostrand, New York, 
6). 
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atmosphere pressure. These values of AF® then 
indicate merely that the reaction (1) could not 
occur in a closed vessel at uniform temperature 
to produce magnesium under these standard 
conditions. If the pressure of magnesium were 
somehow lowered below the equilibrium value, 
the reaction could proceed with the evolution of 
additional magnesium, just as the removal of 
carbon dioxide will permit the dissociation of 
calcium carbonate. In a vacuum tube, any 
magnesium formed will generally condense on 
the walls of the bulb and the action of the bulb 
in thus removing a reaction product is analogous 
to that of a vacuum pump in removing the carbon 
dioxide in the dissociation of the carbonate. Then 
in an actual vacuum tube, as distinguished from 
the closed vessel at uniform temperature, if the 
equilibrium pressure of magnesium vapor in 
reaction 1 is greater than the saturated vapor 
pressure of magnesium at the temperature of the 
bulb wall,* the free energy relationships are such 
that the reaction might proceed at a rate given 
by this equilibrium pressure. By the application 
of well-known methods’ the equation for calcu- 
lating these equilibrium pressures becomes 


In pi=In p2—AF*/RT, 


where /; is the desired equilibrium pressure, and 
p2 is the saturation vapor pressure of magnesium 
corresponding to its state at any temperature T 
for which the free energy values in Fig. 1 were 
computed. As before, AF® is the standard free 
energy difference shown in Fig. 2. Equilibrium 
pressures thus computed are shown by the solid 


* Strictly speaking, the Knudsen correction for thermal 
effusion should be applied, but it is unimportant in this 
case and is neglected. 

5 See for example Noyes and Sherrill, Chemical Principles 
(Macmillan, New York, 1923), pp. 238-242. 
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line curve in Fig. 3. In computing this curve, 
values of p2 were obtained from the vapor pres- 
sure-temperature relationships for magnesium 
given by Landolt-Bérnstein® for temperatures 
from -557°K to 788°K and by Kelley for higher 
temperatures. The data were extrapolated to 
363°K = 90°C assumed as a bulb temperature, by 
the usual method of plotting log p against 1/T. 

In the above considerations, the effect of any 
variation in the state of the tungsten oxide was 
neglected. Its properties are not well known and 
its behavior after formation cannot be accurately 
predicted. A fraction might evaporate, another 
fraction might dissolve in the tungsten and still 
another fraction might combine in some way 
with the magnesia. Experimental tests described 
below showed that comparatively little evapo- 
rated and it was assumed to remain in the coating 
in a state having approximately the standard 
free energy. Any variation in its state would be 
more apt to make the reaction proceed as shown, 
although its neglect probably introduces no 
substantial error. 

In Fig. 3, the vertical dotted line at 1145°K 
intersects the calculated curve of equilibrium 
pressure at an ordinate corresponding to the 


6 Landolt-Bérnstein, Physikalische Chemische Tabellen 


3rd Erg, p. 2434 (1936). 


saturation pressure of magnesium at 363°K. 
This means that with the bulb at 363°K, mag- 
nesium can leave the region of the heater and 
deposit on the bulb at all heater temperature: 
above 1145°K. 


EXPERIMENTAL 


Tube design 


The calculations described above were tested 
experimentally using tubes as shown in Fig. 4. 
The two filaments, each consisting of a tungsten 
ribbon 0.040 X 0.0005 in. supported at the top by 
a suitable spring, were connected in series. Each 
ribbon was coated with a suspension of C.P. 
magnesia,’ the weight of the coating being from 
4 to 5 mg per square centimeter. On one filament 
in each tube, two 0.001-in. molybdenum wires 
were welded 2.75 in. apart for measuring the 
potential drop. Some of the filaments were 
coated with a water suspension of magnesia 
directly as obtained, while others were coated 
with a mixture which had been fired in a molyb- 
denum boat in hydrogen first for 15 minutes at 
1000°C and then at 1700°C. The latter coating 
was much more dense and adhered to the 
tungsten only when suspended in amy] acetate. 
The adherence was then excellent. The unfired 
magnesia shrank to a small fraction of its former 
volume on heating the filaments so that its 
adherence was rather poor. 

The assembled tubes were pumped through a 
tubulation at the base. They were baked at 
450°C for 1 hour and the filaments were operated 
on the exhaust system at the temperature to be 
used in subsequent tests and sealed off. 


Temperature measurements 


The filament in each tube was then operated at 
a fixed current so chosen that the temperature 
range 1100°K to 1900°K was covered. The tem- 
perature of operation was measured pyrometri- 
cally correcting for the spectral emissivity at 
0.654, which was found to be 0.42 by using a 
reflectometer described by Prescott and Morri- 
son,® and the usual bulb correction. The tem- 
perature was also estimated from the resistivity 
of the filaments using the tables of Jones and 


7J. T. Baker, Light Calcined MgO for Radio. 
8 Prescott and Morrison, Rev. Sci. Inst. 10, 36 (1939). 
® Pyrometric Practice (Nat. Bur. Stand., 1921), p. 117. 
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Langmuir.!° The resistivity was determined from 
potentiometer measurements of the filament 
resistance and the weights of the bare ribbons, 
measured after all other tests. When operation 
at the higher temperature caused magnesium to 
deposit on the bulb, the pyrometer measure- 
ments were no longer reliable, and the tempera- 
ture was maintained constant by holding the 
resistance constant. 

After deposit of a suitable film, the tube was 
sealed to a gas analysis system, the break-off 
tip was broken and the free magnesium formed 
was measured by subjecting it to water vapor 
and measuring the hydrogen evolved. 


Results 


The above measurements are listed in Table I 
in such a form that the rate of evaporation of 
magnesium from a tungsten ribbon coated with 
magnesia can readily be computed. Making 
assumptions which appear reasonable, the related 
vapor pressure may be computed using the 
Knudsen equation :"! 


Z=po/(2emkT)}. 


This equation gives the number of atoms Z 
evaporating from a small orifice into a vacuum 
at a pressure p where o =area or orifice, m = mass 
of atom, and k= Boltzmann constant. The use of 
this equation is justified if the dimensions of the 
orifice are small compared to the mean free path 
of the evaporating atoms at pressure p. This 
condition was satisfied in the present work. The 
superficial area of the entire filament was taken 
as o and the reflection coefficient as zero. 
Considerably more doubt exists as to whether 
the pressures obtained by application of the 
Knudsen equation for the measured amounts of 
magnesium represent equilibrium pressures in 


TABLE I. 








MAG- 

NESIUM TUNGS- 
FOUND TEN FouND 
(MG) 


0.53 
0.34 


TuBE TEMPERATURE °K OPERATING 
No. PyroM. REsIS. TIME 


1650 1750 2.50 hours 
1543. 1577 4.25 hours 
1390 1400 118.5 hours 
1745 1836 5 min. 

1257. 1290 743 hours 
1150 1200 982 hours 


(MG) 


0.64 
0.17 
0.18 0.11 
0.15 0.13 
0.u07 + 
0.0015 + 





EC 428 
EC 430 
EC 431 
EE 644 
EE 641 
EC 432 








0 Jones and Langmuir, Gen. Elec. Rev. 30, 354 (1927). 
1 See, for example, Tolman, Statistical: Mechanics (Chem. 
Catalogue Co., New York, 1927), p. 207. 
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reaction (1). Because all magnesium found must 
have formed and diffused away from the inter- 
face between the magnesia and tungsten,. the 
pressure observed could be less than the equi- 
librium pressure either because the area of con- 
tact is less than the superficial area, or because 
of the existence of pressure gradients across the 
coating, or occupation of reaction sites by WOz, 
a product of the reaction. On the other hand, 
the formation of compounds such as magnesium 
tungstate or the ability of the tungsten oxide to 
diffuse away or dissolve might increase the 
equilibrium pressure. No attempt was made to 
resolve these questions and the experimentally 
determined amounts of magnesium converted to 
pressures were plotted directly as points on 
Fig. 3, together with the equilibrium pressures 
calculated from thermodynamic data. The points 
agree with the curve within a factor of less than 
10 for all but one case and might be taken to 
indicate that equilibrium pressures exist over the 
reacting system. 

The total amount of tungsten both as metal 
and oxide found on the bulb walls was also de- 
termined. The bulbs were washed with 0.5N 
KOH and the tungsten contents of the resulting 
solutions were measured spectrochemically by 
comparison with synthetic standards, using a.c. 
arc excitation. The results are given in Table I. 
The total quantity of tungsten found is much 
smaller than the amount to be expected if all 
the tungsten oxide evaporated. Thus the major 
part remains in the reaction zone. 

The agreement between the experimental re- 
sults and the calculations indicates that the 
reaction actually proceeds at approximately the 
rate corresponding to the computed equilibrium 
pressure. This agreement also suggests that the 
free energy of the tungsten oxide, in whatever 
form it actually exists, is not very different from 
that of the assumed form. 

It is a pleasure to acknowledge the help of 
members of the Chemical Research Department 
of these Laboratories, especially for determina- 
tions of free magnesium by Mr. W. G. Guldner 
using methods developed by Mr. L. A. Wooten, 
and spectrochemical determinations of tungsten 
by Mr. A. E. Ruehle. Consultation with Mr. C. 
H. Prescott, Jr., in formulating the theory was 
also very valuable. 
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The Thermal Reaction Between Hydrogen and Oxygen 


III. The Temperature Coefficient of the Steady Thermal Reaction 


O. OLDENBERG AND H. S. Sommers, JR. 
Research Laboratory of Physics, Harvard University, Cambridge, Massachusetts 


(Received February 1, 1941) 


The thermal combination of hydrogen and oxygen, investigated in various vessels (quartz, clean 
Pyrex, Pyrex covered with KCl), follows the Arrhenius law over a considerable temperature 
range in which, according to Hinshelwood, a chain reaction in the gas phase, broken at the wall, 
must be assumed. It is concluded that in this temperature range the chains are not branched. 





I. PROBLEM 


N preceding papers! we discussed certain ob- 

servations on the thermal reaction between 
hydrogen and oxygen, some of which seemed 
incompatible with the present theory of the 
reaction. These observations were largely con- 
cerned with the behavior of the reaction in the 
neighborhood of the second and third explosion 
limits. The present paper deals with the steady 
reaction between these two limits, with main 
emphasis on its temperature coefficient. 

Knowledge of this coefficient is important from 
two points of view. First, it should distinguish a 
nonbranched from a branched chain reaction. 
Second, comparison with the photochemical re- 
action, which has been investigated under the 
same conditions, should yield more detailed evi- 
dence of the mechanism of the thermal reaction. 
This second point might prove very helpful, for, 
though the starting of the thermal reaction is not 
known, the first step in the photochemical reaction 
is definitely the creation of H or O atoms by 
suitable incident light. 

The first definite conclusions as to the charac- 
ter of the thermal reaction were reached by 
Hinshelwood and his co-workers.’ They observed 
that the reaction rate increased with vessel size or 
with addition of an inert gas and so concluded 
that the reaction proceeded in chains broken at 
the vessel wall. Previous to this, chain reactions 
had been identified among photochemical re- 
actions by difficult measurements on reaction 


1Q. Oldenberg and H. S. Sommers, Jr., J. Chem. Phys. 
8, 468 (1940); 9, 114 (1941). 

2C. N. Hinshelwood and A. T. Williamson, The Reaction 
between Hydrogen and Oxygen (Clarendon Press, Oxford, 
1934), p. 36. 
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yields, but this was the first proved case of a 
thermal chain reaction. 

Their work, however, left it uncertain whether 
these thermal chains were branched or non- 
branched. Later, Lewis and von Elbe’ offered an 
interpretation, concerned primarily with the 
various explosion limits, based on Semenoff’s' 
concept of branched chain explosions. Such 
chains are supposed to branch even in the region 
of the steady reaction, but in this range powerful 
chain breaking prevents the explosion. At 
the explosion limits, branching overcomes the 
breaking. 

The rate of a reaction can be expressed in the 
form 


Rate = Roe~#/*7, 


For a reaction which obeys the Arrhenius law, 
the constant E is independent of the temperature 
and is interpreted as an activation energy. 
Experimentally a constant value of £ is recog- 
nized by a linear plot of log rate versus 1/T. The 
activation energy E is constant for a nonbranched 
chain reaction, though here E may receive contri- 
butions from several elementary steps. According 
to Semenoff, however, a different situation occurs 
for branched chains; his mathematical treatment 
leads to an activation energy increasing with the 
temperature, which would be recognized by a 
nonlinear plot of log rate versus 1/T. This is 
plausible, as the probability of a branching, which 
is equivalent to the creation of a new chain, 
increases with the temperature. 


3B. Lewis and G. von Elbe, Combustion, Flames, and 
— of Gases (University Press, Cambridge, 1938), 


4 N. Semenoff, Chemical Kinetics and Chain Reactions 
(Clarendon Press, Oxford, 1935), p. 41. 
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Consequently, if the thermal reaction has a 
constant factor £, true branching cannot be 
present as a branched chain can only be postu- 
lated if’the reaction has a factor E which in- 
creases with the temperature. (Conversely, such 
a factor increasing with the temperature does not 
uniquely define a branched chain, since other 
processes could be responsible (see Section 


III, C).) 
II. METHOD 


We used the method of quiescent mixture in 
which the gases are mixed in the reaction vessel 
in the stoichiometric ratio and the progress of the 
reaction followed by measuring the falling pres- 
sure. A typical run is shown in Fig. 1. We found 
the inhibition period (the time necessary for auto- 
acceleration to bring the rate up to a maximum) 
to be extremely erratic, while the maximum rate, 
averaged over several centimeters of pressure, 
was far more consistent. So it is to be understood 
that whenever we speak of rate, we mean the 
constant maximum rate. (For the packed vessel, 
where the inhibition period was consistently 
zero, we recorded the equivalent of this rate, the 
time for the first 10 percent of the reaction.) 

Our “‘standard”’ size of vessel, in which most of 
our work has been done, was a cylindrical vessel 
of about 5 cm X12 cm. For comparison, we also 
used a few vessels with larger and with smaller 
ratios of surface to volume. We used quartz, 
Pyrex, and KCl-coated ones. 

The oven was provided with a thermostat to 
permit rapid changes in temperature to prede- 
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Fic. 1. Typical measurement of the reaction rate, pressure 
vs. time. 
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Fic. 2. Log rate vs. 1/7. Quartz vessel, uncoated. 


termined values. Temperatures were measured 
by a platinum and platinum-rhodium thermo- 
couple calibrated against a zinc sample furnished 
by the National Bureau of Standards. Pressures 
were measured by a mercury manometer which, 
though heated to prevent condensation of water 
vapor, furnished no mercury vapor to the vessel 
except the 0.001 mm present at room tempera- 
tures. We also checked that additional mercury 
vapor did not affect the reaction. Tank O2 and 
purified tank He were used. Further details of the 
arrangement were given in the preceding paper.' 

Rate measurements were carried out over as 
extended a range as feasible. On the low tempera- 
ture side they were limited by the rate being 
exceedingly slow (one run requiring several 
hours), and on the high temperature side by the 
high pressure explosion. The highest observed 
rate was in excess of 50 percent combination per 
minute. 


III. REsULTs 


A. Quartz vessel 


Hinshelwood and Thompson® were the first 
who measured the temperature coefficient in a 
quartz vessel, but because of the erratic results 
they soon changed to other materials. Neverthe- 
less, since its transparency to ultraviolet light has 
restricted all photochemical measurements to 
quartz, we felt it worth while to spend con- 
siderable time in measuring this coefficient in 
quartz as carefully as possible. 


5C. N. Hinshelwood and H. W. Thompson, Proc. Roy. 
Soc. A118, 180 (1928). 
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Fic. 3. Log rate vs. 1/7. Pyrex vessel, uncoated, 
standard size. 


After a large number of experiments we found 
out that the rates, although highly erratic in a 
new quartz vessel, after many successive runs 
became more consistent. This gradual subsiding 
of the inconsistency led us to the following pro- 
cedure. With a new vessel many runs were made 
in rapid succession. As soon as a better consist- 
ency was obtained, the rate was measured as a 
function of the temperature in another set of 
runs, with but ten minutes of pumping with a 
Megavac pump between successive runs. It is 
plausible that such a procedure provides more 
uniform surface conditions than extensive pump- 
ing and baking between runs. 

The results for this vessel (spherical, 7 cm 
diameter, of newly blown transparent quartz) are 
given in Fig. 2. The first set of nine runs (numbers 
9-17) represents a straight line over a range of 
rates as wide as 1 : 20. (After leaving the vessel 
heated and with 40 cm H.O inside for 15 hours, 
the rates were surprisingly decreased to 18-19). 
The same quartz vessel, which was removed for a 
couple of days for repairs to the furnace, settled 
down after a group of eleven erratic runs to 
another consistent sequence of runs (runs 32-41, 
Fig. 2; 36 was spoiled by a leaky stopcock). In 
this second group of runs we avoided taking 
them simply in the order of increasing tempera- 
ture. Instead, we increased and decreased the 
temperature in an irregular manner. Again a 
good straight line resulted with the same slope as 


before, but it was materially shifted, although the 
experimental conditions, as far as we could 
reproduce them, were unchanged. 

This result, it is true, is open to the Criticism 
that a row of successive runs was selected 
according to its fitting into a simple law, but we 
do not believe that this objection makes the 
experimental result of Fig. 2 invalid. Our reason 
is that we did not discard results within these 
long rows of consecutive runs but only results 
obtained during the initial seasoning period. 
Furthermore, it would be a very strange coinci- 
dence if a straight line representing so long a row 
of successive measurements is only accidental ; in 
order to interpret the straight line by such a 
coincidence one must assume that some property 
of the reaction (like chain branching) leads to a 
deviation from the straight line and that, by 
mere chance, this deviation is just compensated 
by a gradual change of surface conditions. This is 
particularly improbable as the runs were not 
taken in the order of increasing temperature. 

The change of the rate which occurs when we 
leave the heated vessel overnight filled with 
water vapor, or when we remove and remount it, 
but which does not occur within a shorter period 
of time is presumably due to slow activated 
adsorption of gases. This adsorption seems to 
change only the numerical rate factor Ro but not 
the slope of the curve or, in terms of the theory, 


Rate (em/min) 
= 


Ol 
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Fic. 4. Log rate vs. 1/T. Pyrex ‘‘shell,” uncoated. 
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it changes the area covered by ‘“‘active centers’’® 
but not the energy of activation. 

The energies of activation of the reaction in a 
quartz vessel, derived from the two straight lines 
of Fig. 2, are 130,000 and 170,000 cal./mole 
respectively, accurate to perhaps +15 percent. 


B. Pyrex vessel 


In a variety of Pyrex vessels, the results were 
far less ambiguous than in quartz. Also for 
Pyrex the aging of the vessel by many pre- 
liminary runs is far less important. Therefore, we 
extensively investigated the reaction in Pyrex, 
measuring the temperature coefficient at various 
pressures and with several ratios of volume to 
surface. 

In our “‘standard size’’ vessel the characteristic 
curve (log rate against 1/7) is given in Fig. 3. 
The order of the runs is indicated by numbers. 
(Nos. 1-4 were not concerned with the measure- 
ment of the rate; 13 was discarded because of an 
accident with a stopcock.) The pressure was 60 
cm, the temperature range from 493°-543°C, the 
ratio of extreme rates 1 : 60. The points form a 
straight line, with only the first run showing 
much deviation. Considering the difficulty of 
these measurements, in which there is ample 
chance for accident, this cannot be regarded as 
invalidating the linear character of the line. 
There is no deviation from the Arrhenius law up 
to the explosion. The activation energy is 95,000 
cal./mole +4 percent. 

The validity of the Arrhenius law was further 
checked by a set of runs in a vessel of smaller 
ratio of volume to surface. We avoided the 
conventional method of filling the vessel with 
small glass tubes or spheres.,.Such a vessel cannot 
be expected to give the rate that would belong 
to the temperature of the oven because the heat 
developed by the reaction in the middle of the 
vessel is not easily conducted to the walls. 
Therefore, such a vessel would be expected to 
give too high rates because of internal heating. 
Instead, we built a vessel with small ratio of 
volume to surface in the shape of a shell or a 
Dewar flask (0.d. 5.5 cm, distance between walls 
0.52 cm, total volume 160 cm*). 

Figure 4 represents thirty-six consecutive runs 


_°H. S. Taylor, Twelfth Report of the Committee on 
Catalysis (National Research Council, 1940), p. 41. 
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Rate ¢ em/min.) 
Pe 


550 325 500 


1.220 1260 1300 


$210? (°K) 


Fic. 5. Log rate vs. 1/T. Pyrex vessel, uncoated, 
large diameter. 


made with this ‘‘shell’’; no run is omitted, nor 
were there preliminary runs aging the surface. 
Runs 1—13 were made at a pressure of 62 cm, runs 
14-23 at a pressure of 44 cm, and finally runs 
24-36 alternating between these two pressures in 
irregular fashion. Again, for each pressure a 
straight line results, indicating that the Arrhenius 
law governs the reaction. Only one run, number 
34, deviates badly. In particular, there is no 
deviation from the straight line for high tempera- 
ture. The figure shows the very wide range over 
which the Arrhenius law applies (temperatures 
495°-565°C; ratio of extreme rates 1 : 45). The 
energies of activation are the same for the two 
pressures; the diagram gives a value of 82,000 
cal./mole accurate to about +8 percent. 

In a larger Pyrex vessel (diameter 5.5 cm, 
length 12 cm) a series of twenty-six runs was 
made, most of which alternated between 60 and 
42 cm pressure. Some runs with 42 cm of reacting 
gases +20 cm helium, which were inserted at 
random in the same row, will be discussed in 
Section III, E. There were no preliminary runs 
and no run was omitted from the diagram (Fig. 
5). Again the characteristic curves are straight 
lines over a very wide range of temperature and 
rate; for the pressure of 42 cm the rates varied 
over a range of 1 : 40. Only two early runs (3 and 
5) deviate materially from the straight line. 
Again, there is no deviation from the Arrhenius 
law for high temperature. The energy of activation 
results as 74,000 (+4 percent) and 82,000 (+10 
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percent) cal./mole for the pressures of 42 and 60 
cm, respectively ; the difference between the two 
values is within the limits of error. 

Some more measurements were made in a 
packed vessel which originally was built for the 
investigation of the second explosion limit. A 
Pyrex vessel of 5.2 cm i.d. and 16 cm length was 
packed with Pyrex tubing of 6-mm o.d. In this 
vessel the reaction behaved quite differently. 
First of all there was no inhibition period. There- 
fore, we easily measured the rates by the time 
required for the first tenth of the reaction. A 
certain aging during the first runs was evident. 
The rate was proportional to the pressure, while 
in large vessels it was proportional to a higher 
power of the pressure. All these features indicate 
that here we are dealing with the surface reaction 
as would be expected in a packed vessel. These 
results, therefore, have no bearing on the gas 
phase reaction the mechanism of which is our 
main problem. The energy of activation was 
70,000 cal./mole to about 10 percent, perhaps 
smaller than the values obtained in vessels of 
larger ratio volume to surface. Moreover, it is 
evident from the above discussion that in the 
packed vessel the rates observed at high tempera- 
tures may be observed too high because of the 
internal heating of the vessel. If this source of 
error had an effect, the energy of activation of 
the surface reaction would be smaller than the 
value of 70,000 just mentioned. 

The principal result is that in Pyrex vessels, in 
which the rates are more consistent than in 
quartz, the gas phase reaction and, presumably, 
the surface reaction both follow the Arrhenius 
law up to temperatures as close to the third 
explosion as can be investigated. 

The measurements reported (rate versus tem- 
perature measured at different pressures in 
different size vessels) permit a comprehensive 
determination of the order of the reaction in 
Pyrex. The result is an order of 3.8 in the large 
vessel (d=5.5 cm), 2 in the vessel for smaller 
ratio volume/surface (Dewar flask shape) and, 
finally, 1 in the packed vessel. This agrees well 
with the results obtained by Williamson, Pickles 
and Hinshelwood’ in quartz vessels. 


C. Porcelain vessel 


In a porcelain vessel the temperature coeffi- 


7 Computed from the table given in Hinshelwood and 
Williamson, reference 2, p. 36. 
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cient of the thermal reaction was measured by 
Gibson and Hinshelwood.® Their curves (log rate 
against 1/7) are linear for the temperature range 
up to 530°C. Since a straight line with the same 
slope is obtained in a packed vessel up to far 
higher temperatures (580°C), one should at- 
tribute this slope to the surface reaction. In the 
unpacked vessel the curves lose their linear be- 
havior above 530°C; but they reach a steeper 
straight line and continue as such between 560 
and 578°C, the limit of the observations. So they 
indicate that, in this case, the nonlinear part of 
the curve is to be interpreted not by branched 
chains, which would cause an ever increasing 
steepness of the curve, but by a new type of 
reaction coming into play at 530°C and becoming 
predominant above 560°C. This presumably is 
the gas phase reaction. (These temperature limits 
depend, to a certain extent, on the size of the 
vessels.) From the slope of this straight line, 
published by Gibson and Hinshelwood, an energy 
of activation as high as 200,000 cal./mole is 
derived. But one cannot judge how uncertain 
this value is since the authors did not publish the 
individual measurements; the straight line pre- 
sumably indicates only the general trend of the 
curve. 


D. Vessel covered with KCl 


The temperature coefficient in vessels covered 
with KCI has been measured by Prettre.® He 
emphasized that introducing a KCI coat gave far ° 
more reproducible results. By the same tests 
discussed above—effects of pressure and inert 
gases*—he made sure that the reaction observed 
in KCl must be attributed to a gas phase chain. 
The slope of his curves (log rate against 1/7) 
which are linear between 540 and 565°C (not 
observed for lower temperatures because of too 
low a rate) indicates an activation energy of 
95,000 cal. /mole. 

We extended the measurements in the KCl- 
covered vessel in order to check the linear charac- 
ter of the curves and make sure of the tempera- 
ture range over which it applies. However, a 
difficulty arose from gradual changes of the 
surface. Although consecutive runs are rather 
well reproducible, a gradual change of KCI on 
quartz has been noted by Prettre. Even on 

8 C. H. Gibson and C. N. Hinshelwood, Proc. Roy. Soc. 


London, A119, 594 (1928). 
9M. Prettre, J. de chim. phys. 33, 196 (1936). 
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Pyrex some change takes place; after many runs 
the KCl. crystals, although originally easily 
visible, disappear. Another objection to the KCI 
cover of the walls is that at high temperature the 
vapor pressure of KCI, small as it is, may affect 
the chain reaction or, instead, such an effect may 
be due to the Cl, gas if K is built into the glass 
wall. Because of such uncertainties, not very 
many measurements were made with KCI. 

The runs we took gave reasonably straight lines 
in agreement with the results of Prettre, at 
temperatures between 540 and 570°C, leading to 
an activation energy of 115,000 cal./mole-about 
15 percent. 


E. Effect of added helium 


A theoretical discussion (to be published) 
prompted us to investigate the effect of added 
helium. We applied tank helium since its im- 
purities are known to be inert. For further 
purification we passed it through a liquid-air trap 
containing activated charcoal. 

In a long series of runs at random tempera- 
tures, we were completely unable to distinguish 
between a run with 40 cm of 2H2+O2 and 20 cm 
of added helium, and a similar run without the 
helium. Both sets fell on the same linear log plot. 
This result is not very different from those of 
Hinshelwood,' who found that, among the vari- 
ous gases added, helium had the smallest, though 
a still noticeable, effect of increasing the rate. 
Anyway, it is certain that the rate is not de- 
creased by the addition of helium. 

In Fig. 5 we entered only the runs without 
helium in order to avoid confusion. Six runs with 
helium, distributed over the whole temperature 
range, fit very well into the curve for 40 cm 
pressure. 


IV. CONCLUSIONS 


Before discussing the theoretical implications 
of our measurements, we must make sure whether 
we are dealing with the gas phase reaction or the 
surface reaction, first recognized as different 


processes by Hinshelwood and Thompson." 
There seems to be no doubt that the straight lines 
of Figs. 2-5 refer to the gas phase reaction. First 

- N. Hinshelwood and A. T. Williamson, reference 2, 
p 


on C. N. Hinshelwood and H. W. Thompson, Proc. Roy. 
Soc. A118, 170 (1928). 
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of all, this follows from the strong effect of the 
vessel size on the rate which serves as the main 
criterion for a gas phase reaction whose chains 
proceed in the volume and are broken at the wall 
(comparison of Figs. 4 and 5). Also the other 
features of our reactions agree with Hinshelwood’s 
description of the gas phase reaction: the high 
order of the reaction in large vessels and the 
temperature range which in our experiments 
reached even up to 570°C (small ratio volume/ 
surface, pressure 44 cm). In some cases the 
highest temperature for which we measured the 
rate was only a few degrees below the third 
explosion limit. 


A. The starting of the chains 


Hinshelwood proved that the thermal reaction 
consists of chains, continuing in the gas phase, 
broken at the wall ; it remained uncertain whether 
they start at the surface or in the gas phase. 

The material of the vessel has a strong effect on 
the activation energy (up to 170,000 in a quartz 
vessel, 82,000 in a Pyrex vessel of similar di- 
mensions). We must find out whether thiseffect of 
the material is to be attributed to the starting, the 
continuation, or the breaking of the chains. 
Presumably it is not due to the continuation, 
because, according to Hinshelwood, the continu- 
ation takes place in the gas phase (rate increasing 
with vessel size, etc.) and would not be affected 
by a change of surface. Furthermore, it should 
not be ascribed to the breaking because the 
removal of free atoms or radicals—in contrast to 
their creation—is supposed to require no activa- 
tion energy. Hence we should conclude that the 
observed effect of the material on the activation 
energy is due /o the starting process. 

As surface changes are responsible for the 
erratic behavior, it follows that the starting takes 
place at the surface. In particular, it follows that 
the emission of carriers from a Pyrex surface 
requires less energy than the emission from a 
quartz or porcelain surface. 

B. Argument for nonbranched chains 

According to Semenoff,” a gas phase reaction 
which proceeds by branched chains is recognized 
by a deviation from the Arrhenius law. The 
straight lines, however, given above (Figs. 2-5) 
manifest so good an agreement of the thermal 


12 N. Semenoff, reference 4, p. 49. 
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reaction with the Arrhenius law that, within the 
temperature ranges of these straight lines, no 
pronounced branching of chains should be as- 
sumed. This argument seems to us incompatible 
with the present theory of the thermal reaction. 
In this theory the chains are branching, even in 
two consecutive steps, I and II.’ 


If one wanted to adapt this scheme to the Arrhenius 
law, one might exclude continuation of the chain by the OH 
radicals (step I) and so take away the quality of branching 
from reactions II and III. But, presumably, this assump- 
tion would entail a considerable accumulation of OH 
radicals. This seems hardly compatible with our experi- 
ments! on the absorption spectrum of OH radicals in the 
most rapid stages of the thermal reaction; in spite of a 
very sensitive arrangement we completely failed to observe 
any. Moreover, excluding a chemical activity from OH 
radicals would be inconsistent with certain arguments of 
Taylor and Salley which will be discussed in a sub- 
sequent paper. 

13 QO, Oldenberg, E. G. Schneider and H. S. Sommers, Jr., 
Phys. Rev. 58, 1121 (1940). 


14H. S. Taylor and D. J. Salley, J. Am. Chem. Soc. 55, 96 
(1933). 
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Another conclusion from the linear character 
of Figs. 2-5 is that, in the temperature range 
covered, one mechanism is preponderant, or else, 
if several mechanisms are involved, that their 
activation energies do not materially differ. 

Although this argument excludes branched 
chains from the steady reaction, that is, the range 
where the straight lines were observed, it has no 
bearing on the explosion itself which occurs above 
this temperature range (third explosion limit). 
The highest rates we could measure are rather 
close to the third explosion limit; for the origin 
of this explosion it seems entirely sufficient to 
assume internal heating, but branching processes 
are not excluded. 

The Bonhoeffer-Haber scheme seems more 
compatible with the results reported here, as it 
involves no branching of chains. But by the 
comparison with the photochemical reaction, this 
scheme leads to other difficulties which will be 
discussed in a subsequent paper. 





MAY, 1941 


JOURNAL OF CHEMICAL PHYSICS 


VOLUME 9 
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UCHS! reported that considerable expansion 

resulted from the mixing of gases. He worked 
with all possible pairs from the group Os, No, 
CO», N.O, and found expansion as great as 3.5 
percent of the total volume. He did not measure 
the volume change directly, but calculated it 
from the light refraction in the pure gases and 
in mixtures. These experiments were carried out 
under ordinary conditions of temperature and 
pressure. Such large changes in volume on mixing 
do not seem consistent with the behavior of these 
gases when pure. Little work is reported by 
others that would indicate the validity of his 
conclusions. Leduc? found that air had the same 
density as calculated from its composition and 
the density of nitrogen and oxygen. He found 
contradictory results in two experiments on the 
mixing of carbon dioxide and nitrous oxide. 


1 Fuchs, Zeits. f. physik. Chemie 92, 641 (1918). 
2 Leduc, Ann, Chim, Phys. 15, 5 (1898). 


Richardson and Woodhouse’ reported that this 
mixture expanded substantially as reported by 
Fuchs. Others have reported volume changes at 
high pressures, for instance Amagat? at pressures 
around a hundred atmospheres, and Gillespie® at 
pressures ranging from thirty atmospheres up, 
the data being extrapolated back to zero pressure. 

The present work was undertaken to see if by 
direct means the conclusions of Fuchs could be 
confirmed. Such conditions are of practical sig- 
nificance in numerous ways. Thus, for instance, 
the analysis of gas mixtures must involve correc- 
tions for the effect, as pointed out by Kobe.® The 
same gas mixtures that Fuchs used have been 
investigated over the complete concentration 
range, all at 25° and at atmospheric pressure. 

3 Richardson and Woodhouse, J. Am. Chem. Soc. 45, 
2368 (1923). 

4 Amagat, Comptes rendus 127, 88 (1898). 


5 Gillespie, Phys. Rev. 34, 352 (1929). 
® Kobe, Ind. Eng. Chem., Anal. Ed. 3, 262 (1931). 
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The apparatus is shown in part in Fig. 1. 
It consists of two bulbs, A and B, of about 70 cc 
capacity each. Bulb A is calibrated by filling 
with mercury from cock 1 to the mark ¢ and 
weighing. Connection is made at 6 to a mercury 
leveling bulb. The end a is sealed to the buret 
and manometer previously described.’ The bulbs 
are immersed in a thermostat, and the buret and 
manometer are water-jacketed and maintained 
at 25°. 

To carry out an experiment, the bulbs and 
buret are first filled with mercury by raising the 
leveling bulbs, and opening the cock at the top 
of the buret to the atmosphere. Connection is 
then made between the gas source and the top 
of the buret, and gas drawn in by lowering the 
leveling bulbs. When the bulb A is filled with gas, 
the mercury meniscus is set exactly at the mark 
c and cock 2 closed. When the gas has had time 
to reach thermostat temperature, the pressure is 
accurately adjusted to 760 mm by manipulation 
of the mercury level in the buret. Then cock 1 is 
closed, the buret filled with mercury, and the 
other gas of the mixture drawn into the buret, 
to any desired height. After temperature equi- 
librium is reached, the pressure is adjusted as 
before, cocks 1 and 2 are opened, and by manipu- 
lation of the leveling bulbs, the gas is passed 
back and forth several times between the buret 
and the bulbs. Mixing is completed in this way 
very quickly. The mercury level is again set at c, 
the pressure adjusted to 760 mm, and the volume 
change noted by the change in buret level. 

The gases used came directly from commercial 
cylinders. The CO. was 99.8 percent and the 
NO was 99.7 percent pure. In view of the 
results, the impurities can have little effect. 








Fic. 1. Apparatus used (in part). 
7 Markham and Kobe, J. Am. Chem. Soc. 63, 449 (1941). 
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Fic. 2. Volume changes on mixing gas pairs. 


The error caused by temperature fluctuations, 
and the error in reading the pressure and volume, 
have a probable maximum of about 0.02 cc of 
gas. This amounts to an error in AV of about 
0.02 percent. 

The results are presented in Fig. 2. The AV 
found, instead of several percent, as reported by 
Fuchs, was only several hundredths of a percent 
of the total volume. Hence the total effect was 
not far outside the maximum possible error. 
The plots show, however, that in most cases the 
results are consistent with each other within less 
than 0.01 percent. 

It is significant that with all gas pairs, except 
oxygen and nitrogen, a small but definite in- 
crease in volume was consistently observed. This 
increase.is about one percent of that reported by 
Fuchs, and indicates that the change that he 
observed in the refraction had a cause other than 
the change in volume. However, Fuchs reported 
that AV was large or small, depending on the 
difference in the critical temperatures of the 
gases. The same statement applies qualitatively 
to the results presented here, in which the 
O2—Ne and the CO.—N,O mixtures show the 
smallest AV. 

Calculation by the van der Waals equation of 
the increase in volume when carbon dioxide 
expands from one to one-half atmosphere shows 
an increase of 0.26 percent. This should be the 
maximum, as the presence of other gas molecules 
should retard expansion. The volume changes 
determined are so small that the additive volume 
law of Amagat is as well satisfied at ordinary 
temperatures and pressures as the other gas 
laws. At higher pressures, where others have 
reported substantial volume changes on mixing,* 5 
this law may not hold. 
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Solutions of Long Chain Compounds 


MAvRIcE L. HuGcGIns 
Kodak Research Laboratories, Rochester, New York 
April 4, 1941 


ppt and Rushbrooke! have derived statistically 
equations for the activities of the constituents of a 
dilute binary solution, in which the solute molecules are 
elongated, each occupying twice the volume occupied by 
a single solvent molecule, with the molecules of the two 
kinds otherwise so similar that there is no heat of mixing. 

Extending their methods to the idealized case of a 
dilute solution of long chain molecules, each composed of 
randomly oriented ‘‘submolecules,” similar in nature and 
size to one of the solvent molecules, I have derived the 
following relations: 


a, = N,*!!" exp(—k, N2*) 


N.* N.* 
Pe “ exp| 2hin(INs* ) -)} 


Here, a; and az are the activities of solvent and solute, 
Ni* and N,* are related to the mole fractions N; and N2 
by the equations 


a:= 


N, nNe 
Ni+nN2 ~ Mi+nN,’ 


n is the number of “‘submolecules” in the solute molecule 
chain, and k; and ke are constants, equal to 0.45 and 11, 
respectively, for a close-packed distribution of solvent 
molecules and solute submolecules. 

These activity equations can, of course, be used to 
derive expressions for osmotic pressure, freezing point 
lowering, solubilities, etc. Application to osmotic pressures, 


Ni*= and N,*= 
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for example, leads to the relation 
n=“F(- InN,*!/"+k, N2*?) 


AEM t+ hintN), 


where R is the molal gas constant, T the absolute tempera- 
ture and V, the molal volume of the solvent. Except for 
the constant ;, this last equation is equivalent (for my 
hypothetical model) to 


n= A701 +(Q-1)9s 1 


just published by Powell, Clark and Eyring,? in which 
Q is the number of segments into which the solute molecule 
is divided and ¢2 is the volume fraction of the solute in 
the solution. Their equation was derived by assuming 
“that the osmotic pressure is determined by the effective 
mole fraction” and “that the probability of a polymer 
molecule’s moving in segments is a linear function of its 
environment as expressed by the volume fraction of 
polymer.” 

According to the above equations, II/c increases ap- 
proximately in proportion to c (the concentration in g/cc 
or similar units), a relationship which has been known 
empirically for some time.’ On the theoretical side (besides 
the papers already mentioned), Haller* has presented a 
kinetic interpretation and Meyer® has recently discussed 
the subject from a statistical point of view, but in a 
qualitative manner. 

Details of the present treatment and application to 
other properties and to real solutions will be given at the 
Wilder D. Bancroft Colloid Symposium in June. 

1R. H. Fowler and G. S. Rushbrooke, Trans. Faraday Soc. 33, 
ar | (1937). 

. E. Powell, C. R. Clark and H. Eyring, J. Chem. Phys. 9, 268 
aa. Mark, Physical Chemistry of High Polymeric Systems (Inter- 
oust Publishers, New York, 1940). 

W. Haller, Kolloid Zeits. 56, 257 (193 


1). 
ok. H. Meyer, Zeits. f. physik. Chemie B44, 383 (1940); Helv. Chim. 
Acta 23, 1063 (1940). 





Erratum: Spectroscopic Evidence for Hydrogen Bonds: 
Comparison of Proton-Attracting Properties of Liquids. 
IV. 

(J. Chem. Phys. 9, 215 (1941)) 


WALTER GorpDy 
Mendenhall Laboratory of Physics, Ohio State University, Columbus, Ohio 
April 10, 1941 


The captions of Figs. 2 and 3 are interchanged. 
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